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t-h ! Abstract 

m ■ 

Color transparency is the proposal that under certain circumstances the strong inter- 
^ actions can be controlled and in some cases reduced in magnitude. We give a com- 

prehensive review of the physics, which hinges on the interface of perturbative QCD 
with non-perturbative strong interactions. Color transparency is expected to occur 
^ . in many kinds of quasi-exclusive reactions with either electron or hadrons beams. 

j> \ We review the interplay of color transparency with nuclear filtering, which is the 

conversion of quark wave functions in hadrons to small transverse space dimensions 
5_i ■ by interaction with a nuclear medium. A complete description of the phenomena 

as a multi-scale QCD process is given in terms of light-cone matrix elements. We 
also review a number of other approaches, including pictures based on modeling 
the time evolution of hadronic wave-packets and the use of hadronic interpolating 
states. Spin plays an intrinsic role in hard exclusive reactions and the role of spin 
in testing and understanding color transparency and nuclear filtering is discussed at 
length. We emphasize the use of data analysis procedures which have minimal model 
dependence, while comparing and contrasting various models. We review existing ex- 
perimental data and the experimental program planned at various facilities. We also 
discuss making use of the subject's strong scientific complementarity and potential 
for progress in exploring hadron physics at current and future facilities, arguing that 
the combination of different experiments can lead to a broad new program to study 
the strong interactions and hadron structure. 
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We will begin our survey with an overview of the basic idea of color transparency 
and how the concepts have evolved. Following an introduction, we will present the 
basic idea of "survial of the smallest" , then "back up" to add explanation, and answer 
some common questions. The Section is closed with a brief discussion of a program 
in strong interaction physics at current and future facilities. 

1 Introduction 

Color transparency is the predicted phenomenon of reduced strong interactions under 
certain conditions. The prediction hinges upon what is known about hard scatter- 
ing processes, the coherent cancellation of perturbative interactions, and a body 
of experience in strong interaction physics. Color transparency presents many new 
experimental opportunities as well as many interesting theoretical paradoxes. It chal- 
lenges our understanding of the quark and gluon-based picture of strong interactions 
as it confronts an established picture based on nucleon and meson degrees of freedom. 
While the new and old pictures are complementary, it seems that color transparency 
would never have been anticipated in conventional strong interaction physics, because 
it represents an entirely new approach to how hadrons interact. 

Recall that the proton-nucleon inelastic cross section is about 30 — 35 mb in the 
energy region of 2 — 30 GeV. The mean free path before having an inelastic collision 
inside nuclear matter should be about 2 fm, so that a proton would almost never 
escape intact from a nuclear target in the conventional theory. However, QCD says 
that under certain circumstances we can prepare a state to make it pass through many 
Fermi of nuclear matter. For all practical purposes, sometimes the strong interactions 
can be turned off. This causes a drastic conceptual shift in the interpretation of the 
physically realized strong interactions. It is a direct consequence of hadrons being 
composite objects made of internal degrees of freedom, the quarks. 

Color transparency theory is based on the rather well established methodology of 
perturbative QCD. However, the experimental situation is not settled; there is much 
discussion at this time whether color transparency has actually been observed, or 
will be observed, at laboratory energies and in realistic experimental circumstances. 
The subject is quite lively and controversial; no matter which way it is presented 
it sometimes seems to contain fantastic and unbelievable ingredients. If research in 
this area will probably test the limits of perturbative QCD, it will also explore more 
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about the "real" strong interactions using QCD as a tool. 

We hope to learn fundamentals of the strong interactions by having a new uni- 
verse, namely reactions inside nuclear targets, to compare with the usual universe of 
experiments made in "free space" . The new methodology, of using the nucleus as a 
test medium to compare and contrast with the "control" reactions of hadrons in free 
space, may eventually be one of the most significant advances ever in learning about 
strong interaction physics. At the present time, color transparency plays a central 
role in the growing, broad subject of the interface of perturbative (p)QCD with the 
non-perturbative structure and interactions of hadrons. The ultimate goal of this 
and other current study at the interface is: 

to solve the physics of the strong interactions. 



1.1 The basic idea: survival of the smallest 



Hadrons are made of quarks and gluons. While hadrons are asymptotic energy eigen- 
states, the quarks and gluons are not. It follows that hadrons fluctuate between 
different quark and gluon configurations when they interact; the process of measure- 
ment involves time dependent dynamics which we wish to understand. To observe 
color transparency one needs to measure a hadron when it has fluctuated to a con- 
figuration which will have small interactions, by definition. 

The key to having unexpectedly small interactions in a gauge theory is coherence, 
in which emissions from different parts of a composite system can cancel with de- 
structive interference. Because coherence is ubiquitous in electrodynamics, there are 
several instructive precedents for color transparency. An early theoretical prediction 
by King|80| was studied by Perkins [ |lld|1 . He observed tracks in emulsion following 
e + e _ pair production in cosmic ray experiments, and found reduced ionization in 
the spatial region close to production of the pair, where the opposite electric charges 
were close together in space. This is a case where the continuum spectrum of the 
e + e~ system is important, and the bound state contribution irrelevant. Anomalous 
transmission of electrically neutral systems such as positronium through thin foils 
has also been observed. In this case properties of the bound state have been used 
to make order of magnitude estimates. Comparing these two examples, it should be 
clear that both bound state and continuum states can be important in coherence- 
the main thing is that an important region of the wave functions produce canceling 
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results. 

In elementary scattering theory, coherence also has a major effect on the calcula- 
tion of cross sections. The total cross section for scattering of an electron by a bare 
Coulomb field is infinite. However, the elastic scattering of low energy electrons on 
a neutral Hydrogen atom has a finite, geometrical cross section, given approximately 
by 

d " - <M? ; /to)- eV ' (2/0)< 



dn ' <f V [(2/a) 2 + g 2 ] 2 

a tot = 47ra 2 , 



where a is the Bohr radius. This comes from the Born approximation to the screened 
Coulomb potential of a central proton surrounded by an electron cloud. [|117| This 
illustrates the rule that if a neutral system in a gauge theory is spatially small, 
then its reaction cross section will be small. Obviously this is not the only way 
for coherence to reduce interactions, so the argument cannot be reversed. A small 
cross section does not necessarily imply a spatially small object. We will see the 
consequences of this later. 



In QCD useful calculations have been done [93, |108| , |59| for the analog of neutral 
"atom-atom" scattering. As in QED the cross section is found to be small if associ- 
ated dipole moments are small, cutting off the amount of radiation exchanged. The 
scattering cross section a(b) for the region where two fermions in neutral "atoms" are 
separated by a transverse position separation b scales like b 2 . This suggests that re- 
duced interactions, and color transparency, should occur for configurations of quarks 
which have a transverse separation small on the typical hadron size scale of about 
1 Fm. As we will see, defining these regions is somewhat subtle, and there are dif- 
ferent renditions in the literature of what this means. Observing color transparency 
may require tuning of the experiment, so that particular regions of the quark wave 
function can make sufficiently large contributions to an observable process. 

It follows that the search for color transparency is a search for short distance 
regions where simple color coherence will be found. According to current theoretical 
work the short distance regions might dominate in certain hard exclusive reactions. 
("Hard" means more than high energy; it means that every invariant momentum 
transfer Q 2 , the Mandlestam variables s, t and u for 2 — > 2 reactions is big compared 
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to a few GeV 2 . "Exclusive" means that all the momenta of all participating hadrons 
are measured.) The relevance of short distance properties to hard exclusive reactions 
is a very important question, controversial for something like twenty years. By now, 
the theoretical statement of short distance dominance at large Q 2 is controversial 
mainly in detail: how big does Q 2 have to be for an asymptotic argument to apply? 
In what experiments does it work? A second ingredient from hard scattering pQCD is 
that at large Q 2 the most important contributions come from the minimal Fock space 
number projection of the hadrons: 3 quarks in a proton, 2 in a meson. |l8|, Again, 
how large does Q 2 have to be to enforce this? How reliable are the perturbative 
estimates? 

The concepts of coherence and the "dominant region" of a wave function are in- 
trinsically quantum mechanical, and depend on the measurement conditions. "Short 
distance" has an intuitive meaning, but also has a precise definition in QCD in the 
impulse approximation, as we discuss below. Because of all these conditions, the 
requirements for color transparency are far from guaranteed in general, and may or 
may not be obtained under conditions such as diffractive scattering. As for semiclas- 
sical or classical arguments, we will find that a correct view of color transparency 
cannot be obtained from their consideration. 

In the classic hard scattering conditions, the typical quarks selected are separated 
by a transverse distance b that gets smaller as Q 2 increases, with the important 
integration region being b 2 = 1/Q 2 - To test whether this is the nuclear target 

provides the ideal detector. Let an experiment arrange for a hard exclusive reaction 
right inside a nucleus, which will be standing there to attenuate the participating 
particles. An example experiment would be knocking a proton out of a nucleus with 
a fast incoming proton, the reaction pA — > p'p"(A — 1), where in this case Q 2 ~| t \. 
A signal of color transparency is a reduced attenuation cross section in the nuclear 
target as Q 2 increases. Ideally, we would measure the cross section da/dt at fixed cm 
angle as a function of Q 2 , extract the nuclear attenuation rate from the data, and see 
the typical "size" of the quark separations decreasing rapidly with Q 2 . This is "part 
one" of the basic idea, as originally expressed by Brodsky and Mueller p0|, |100| . 



The nucleus is "passive" if we think of it as measuring the event; it is "active" if it 
plays an essential role in giving us the configurations sought. "Part two" of the color 
transparency story exploits this aspect of quantum mechanics; this is the concept 



we have denoted "nuclear filtering" . [TO, |119|, |120|| If a spatially large region would 
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contribute much, it would also proceed to have large inelastic interactions, and the 
experiment will simply reject it with a trigger on quasi-elastic data. This is quantum 
mechanics, which depends crucially on the experimental conditions. Nuclear filtering 
again is utterly quantum mechanical and cannot be replaced by a classical picture. 

There is a precedent to nuclear filtering in Sudakov effects, reviewed in Section 
4. In the past, Sudakov effects have been considered in "free space". But free space 
is a physical medium, much like the nucleus because most of the probability ends up 
in highly inelastic events. In the free-space Sudakov analysis, the rare and special 
amplitudes which remain exclusive, (through inability to radiate away gluons), must 
come from limited transverse separations between the charged particles (or small 
dipole moments in QED language). This tends to select short distance even when 
kinematics might not suffice. Similarly, in a nuclear medium, only the amplitudes 
which have small interactions will survive to be measured. The regions involved must 
strongly select short distance between quarks, because the nuclear medium is dense 
(compared to free space.) There are therefore strong parallels between traditional 
Sudakov effects and the pQCD treatment of color transparency which we will explore. 
An extra parameter, and an extra experimental handle, is the nuclear target length 
Ra ~ l.2A l ^Fm, which regulates the transverse smallness for partons not to have 
an interaction. 

With large Q 2 and large A, we have two excellent tools that tend to select the 
short distance regions of quark wave functions, and study them. 

1.2 Back Up! 

In presenting the basic idea, we have skipped over many essential details, which will 
be clarified below. We have also presented the idea from the picture in the quantum 
mechanical basis of quarks, as originally conceived and formulated. Needless to say, 
theory can employ different bases, leading to seeming disagreements in the physical 
picture. The basis of hadrons as a complete set has been extensively used, and 
deserves a discussion on an equal footing. 

In perturbing the system to measure color transparency, it is quite sensible to 
consider the resulting system as a superposition over asymptotic hadronic states 
rather than the fundamental Fock space. In this basis, pQCD has to be reinterpreted. 
The most common interpretation is that traditional strong interaction physics applies. 
This is a dynamical assumption, implying that under conditions of color transparency 
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the decreased strong interactions will be due to destructive interference among the 
different hadronic modes in a wave packet. This is far from trivial: there is no 
systematic way to go from predictions in the quark basis to those made in a basis of 
asymptotic states. Thus, models in the hadronic basis represent new hypotheses to 
be tested separately from the quark picture. Underlying this is a broader conceptual 
question: to what extent can asymptotic hadrons actually be used to imitate the 
predictions of QCD? The answer is unknown. 

Current calculations in the hadronic basis are intrinsically neither more nor less 
rigorous than those using the quark and gluon Fock space. In the hadronic basis 
one has to rely on a postulate of completeness, and a host of unknown, uncalculable 
matrix elements. Completeness is a pre-QCD postulate from the days of the S-matrix. 
However, it has never been quite clear how the collective coordinates, the hadrons, 
are to reproduce the internal coordinates, the quarks, in the system with confinement. 
Setting aside this question for now, the hadronic viewpoint's advantage is that one 
can use some measured information about the interactions of hadrons to model the 
system. This is powerful, and many theorists have exploited it. A corresponding 
disadvantage is that any quantum mechanical superposition being proposed has to be 
"tuned" to reproduce the perturbative results, where they apply, and as much as they 
are understood. It has often been assumed that a spatially small wave packet made 
by superposing hadronic states must reproduce the quark picture. But as we shall 
see in Section (2.2) below, tuning the hadronic wave packet is not straightforward: 
contrary to some expectations, the pQCD dominant region of b 2 < 1/Q 2 is not the 
same as a wave packet with < b 2 >< 1/Q 2 . 

1.3 Common Questions 

Why are exclusive processes emphasized? Couldn't one observe color transparency in 
all kinds of reactions ? 

The answer to this is maybe so, but probably not. Inclusive reactions sum over 
channels and disparate kinds of configurations of the quarks inside the hadrons (or 
hadronic final states). In the sums there are many color flows, and different degrees of 
inelasticity. Perhaps some of the configurations are "transparent" , but it is reasonable 
that these should be a small part of most inclusive cross sections. An active theoretical 
occupation is to seek "solid gold" reactions which are color transparent while relaxing 
the difficult experimental conditions of exclusivity. 
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On the positive side, exclusive processes are a wonderful frontier. Their study 
is well justified as a field in itself. Exclusive processes are difficult to measure, but 
this is because there are many, many channels, and much work to be done. At 
large momentum transfer, they involve projections of hadrons onto just a few quarks. 
These are exactly the projections we are interested in learning about on the interface 
of non-perturbative QCD physics. 

Even if the perturbative treatment becomes self consistent, isn't there still room 
for doubt in its predictions? 

Absolutely yes. Many interesting questions - such as what happens to a hadron's 
"cloud" of soft pions and/or soft gluons, the way it sits in the chirally broken vacuum, 
and the way a hadron reforms after being struck hard - are all questions beyond 
perturbation theory. It is conceivable that perturbation theory will fail. But if color 
transparency does not work, then we will have observed where perturbation theory 
fails in a region where it might have worked. This is a way to learn the fundamentals 
of the strong interactions better. 

Why are small transverse separations between quarks involved in hard scattering 
reactions ? 

This is a straightforward but detailed technical question. The true answer is: 
the association of large momentum transfer with short distance is not always true. 
Understanding when short distance between quarks is involved or not is a current 
research problem. One must check for it on a reaction-by reaction basis. By studying 
it, we hope to learn more about the way the quarks themselves are arranged. It is 
important to note, however, that we should not, and do not, always insist on having 
short distance. If different distances are involved in different reactions, we can take 
a snapshot of the quark wave functions to see how the hadrons are made. 

1.4 Exclusive Reactions: a Program for Solving Strong In- 
teraction Physics 

QCD is now more than 20 years old, and maturely established as the fundamental 
origin of the strong interactions. The time is ripe to begin using that theory, and 
understanding it fully, to solve the many outstanding puzzles of hadron structure 
and nuclear interactions. Central to these puzzles is confinement, not understood 
theoretically, a topic in which experiment must play a central role. Physics needs a 
bridge between the non-relativistic picture of nucleons making up a nucleus on the 
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one hand, and the ultra-relativistic perturbative approach used to establish QCD on 
the other. Color transparency provides a means to make this transition. The subject 
interpolates between quarks interacting weakly at short distances to the "real" strong 
interactions of hadrons at long distances. 

To study the interface of perturbative and non-perturbative QCD, and move 
toward the goal of solving the physics of the strong interactions, a re-focusing of 
experimental emphasis will be needed. The key element in understanding hadron 
structure will be exclusive reactions. By systematically studying exclusives, impor- 
tant questions of how quarks are confined in hadrons will become answerable. Many 
questions are extraordinarily interesting: how "big" in space are quark and gluon 
wave functions spread out, hadron by hadron, Fock-space channel by channel? How 
good is a few-quark approximation, and at what value of Q 2 1 How do gluons dress 
quarks in detail, and how does their color flow as mesons are exchanged? What is 
the role of quark orbital angular momentum, in making up the spin of hadrons and 
in exchanges between them? A singular thing about these questions is that they 
cannot possibly be answered in the bulk, inclusive-averaged experimental framework. 
Instead, theory and experiment must study channel by channel, spin by spin, in as 
much quantum mechanical detail as possible. 

The framework of hard exclusive reactions, along with nuclear filtering, promises 
to make such a study systematic and theoretically well defined (Sections 3 and 4). 
A bonus from the "parton revolution" of the past 20 years is a strong theoretical 
structure in defining what matrix elements are measured, how, and when. Entirely 
new matrix elements are measured with exclusive interactions, and can be used con- 
structively to compare and make cross-links between diverse reactions. This is future 
science we hope to see. But at the moment color transparency is still an emerging 
field, with active discussions and widely different interpretations. A calculation now 
using meson exchange might be better than perturbative QCD for some reactions, 
or maybe worse. Diversity of ideas and interpretations is part of the fun: the sub- 
ject is open and interdisciplinary, without any picture yet commanding the field. In 
this report we have done our best to unite the valid quantum mechanics used in 
nuclear physics approaches to the much different Green function methods of pertur- 
bative QCD. From these interdisciplinary roots a new discipline of strong interaction 
physics is ready to spring. 

Investigating the strong interactions in these new ways will require several exper- 
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imental technologies. What is needed now for a program in quark-hadron structure 
is development of experimental facilities. Every different beam and target has its ad- 
vantages. Hadron-initiated reactions are of extraordinary interest in identifying the 
color, flavor and angular momentum flow between particles. For fixed-angle reactions, 
a beam energy in the 5-30 GeV region is ideal, with high flux, high acceptance, and 
state-of-the-art polarization and spectrometers the most important requirements. As 
will become clear (Section 3), the comparatively large data taking rates of hadron ini- 
tited reactions, and the fact that several particles interact with the nuclear medium, 
also gives hadron initiated reactions a special attractiveness. The Indiana LISS con- 



cept |89j includes many of these features and seems ideal. We refer the reader to a 
very thorough study preliminary to a formal proposal [jB9||. 

At the same time, facilities at BNL are expected to continue providing pioneering 
measurements, which so far are the foundation of our belief that color transparency 



has been observed ||28|| . A second round of hadron beam experiments has been ap- 
proved and a new detector, named EVA, with much higher acceptance has been taking 



data with beam energies in the 6-20 GeV range for about one year|24|]. Due to the 
revolutionary detector design and increased beamtime, this is hoped to expand the 
amount of data taken by a factor of perhaps 200, allowing a larger energy range and 
an analysis at different scattering angles. It is also planned to study meson-nucleus 
scattering. In particular, the reaction it A — > p(A — l)p at 90° is very interesting as a 
test of short distance physics expected from nuclear filtering (see Section 3.5). 

Many interesting results can come from the diffractive regime. Already there is 
valuable data [TTL |5{J from Fermilab experiment E665. This also provides evidence for 
observation of color transparency, and will be reviewed in Section 3.3. The Hermes 
detector |6^] at HERA is also beginning operation. It is a nuclear gas target installed 
on the 25 — 30 GeV electron (or positron) beam. The use of light, possibly polarized 
nuclei is foreseen in the next few years. The experiment may enable extension of 
the E665 data (see Section 3.3) on p meson diffractive production at moderate Q 2 
values to smaller values of energies 10 < v < 22 Gev. Hermes' small luminosity 
will only allow integrated measurements. Due to limitations of energy resolution, the 
experiment cannot assure that diffractive events are not mixed with inelastic events, 
much like FNAL. In the near future one does not expect detection of the recoiling 
proton, but nevertheless there will be much to learn from the results. 

Electron-initiated reactions also have extraordinary appeal: the electromagnetic 
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probe is the cleanest tool to investigate the structure of hadrons at short distance. 
Accelerator technology has progressed dramatically. Despite the relatively small sig- 
nal - due to the size of a em - it is now practical to conduct precise experimental 
measurements of exclusive reactions with electron beams. The breakthrough is the 
availability of high duty factor and high luminosity electron machines at Mainz, Bonn 
and now CEBAF. CEBAF has accepted several proposals to study color transparency 
and hopefully will continue to add more in the future. There have been many inge- 
nious ideas; elsewhere Section 3) we will point out that the high precision facilities 
of CEBAF have special advantages which have not been exploited. It is simply not 
established whether color transparency will be observed or not at CEBAF energies. 
For this reason, and also due to the intrinsic interest of the subject, it would seem 
logical to answer experimental questions answerable at CEBAF by experiments. 

The ideal choice of the energy range for an electron accelerator suitable for color 
transparency studies is fixed by three constraints. (1) One must have sufficiently 
high energy and momentum transfer to describe the reaction in terms of electron- 
quark scattering. High energy creates a very fast process where the struck quark 
can propagate quasi-free. High momentum transfers are necessary to probe short 
distances. (2) The energy of the incident electron beam should be tunable to match a 
characteristic interaction time r to the diameter of the nucleus. Starting from the rest 
frame time r G ~ 1 Fm/c and taking into account a Lorentz dilation factor 7 = E/M 
this means a time r of several fm/c's in the laboratory. If the energy transfer is 
too large, the building-up of hadrons may occur outside the nucleus which can then 
no longer be used as a microscopic detector. Theory cannot reliably determine the 
interaction time yet, which remains an experimental question to be investigated. (3) 
Charm production, if desired, requires a minimum electron beam energy of 15 GeV 
to have reasonable counting rates. 

To get a large enough rate and avoid a prohibitively large number of accidental 
coincident events a high duty cycle is imperative. Finally, good energy resolution 
is necessary to identify specific reaction channels. A typical experiment at 15 GeV 
(quasielastic scattering for instance) needs a beam energy resolution of about 5 MeV. 
At 30 GeV the proposed experiments only require one to be able to resolve pion 
emission. These requirements are satisfied by the ELFE proposal, an exciting and 
ambitious prospect to which we refer the reader for more detail. j|, [5| 

There have been a large number of studies of accelerator and detector require- 
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ments for a future program in hadron physics. In general, these studies point at 
measurement capabilities in quality of tracking, particle identification, and high ac- 
ceptance which have previously been used in high energy physics experiments, and 
are comparatively new to the nuclear /particle physics interface. We refer the reader 
to several recent specialized meetings. pSI Sufficiently high momentum resolution 
(5 x 10~ 4 ) and high luminosity (10 38 nucleons/cm 2 /s) can be achieved by magnetic 
focusing spectrometers. For semi-exclusive or exclusive experiments with more than 
two particles in the final state, the largest possible angular acceptance (~ 4tt) is 
highly desirable. Quality and reliability of large acceptance detectors have improved 
substantially in the last two decades. Research and development studies, however, 



are still needed, as the literature will show. [95]. 

We now turn to our main topics, beginning with a deeper look at the basic quan- 
tum mechanics in Section 2. Section 3 is concerned with defining the process exper- 
imentally, as well as summarizing existing experimental results, and different ways 
to extract a signal. Many suggestions for future experiments also appear here. Sec- 
tion 4 is a small "handbook" for the beginner to understand better the methods of 
pertubative QCD. There are also some new results that may interest the experts. 
Hopefully the reader will see the link between this physics and the (mostly non- 
relativistic) quantum mechanics of Section 2. Section 5 is concerned with reviewing 
diverse models of propagation. 
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2 Quantum Mechanical Basics 



Color transparency is an intrinsically quantum mechanical phenomenon with some 



semi-classical interpretations. The correct picture has been controversial. ||122| In this 
section we will examine these issues, mostly using elementary quantum mechanical 
arguments. Important points include the motivation for short distance in hard exclu- 
sive processes, the different roles of elastic and inelastic channels, the expansion time 
scale problem, and use of the hadronic basis as a substitute for the quark basis. To 
make this section more accessible, we have separated detailed discussion of QCD into 
Section 4, and used non-relativistic concepts when practical. As we will see, some 
basic discrepancies between different theoretical approaches can be traced to different 
underlying quantum mechanical assumptions. The way in which these assumptions 
are implemented in specific calculations is also a separate topic, to appear in Section 
5 on Models of Propagation. 

It turns out that there are plenty of basic quantum mechanical issues to resolve 
before turning to detailed calculations. We hope that this review of elementary ideas 
can be understood by the beginner, while refreshing the understanding of the experts. 



2.1 What is the right picture? 

It is often asserted that large momentum transfer in a reaction will select short 
distance by the uncertainty principle, AQAX > 1. A widespread notion also exists 
that simply absorbing a virtual photon with large Q guarantees some kind of short 
distance. But certainly one must specify what is going to be measured! There is 
no dispute that a hard inclusive process in QCD such as deeply inelastic scattering 
probes short distance aspects of hadron physics. However, it worth reviewing the 
role of short distance in exclusive reactions, and in some cases of current interest in 
QCD, the bland assertion of generic short distance can be utterly wrong. 

Exclusive versus inclusive processes 

Consider a typical exclusive process such as 7* + p — ► p'. The matrix element 
involved can be represented by two form factors: 

(p + Q I J? m I p) = u(p + Q, /)[F 1 (g 2 ) 7 ^ + T^F 2 (Q 2 )a^QMP, s) (2.1) 
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The process is illustrated by Fig. (2.1). 

Suppose we do not interpret this process yet in terms of "quarks" or QCD. On its 
own, objectively, there is no kinematic information in the matrix element to indicate 
short distance. Without a model for the substructure, the matrix element simply 
stands alone (as it did in the 1960's). The standard assertions of short distance are 
motivated by fairly detailed properties of the proton as a composite object. These 
are specific to QCD and the proton as made of quarks. 

Depending on the kinematics, the proton can also absorb a virtual photon and 
"explode" into a multitude of final states in a hard inclusive process, namely deeply 
inelastic scattering. The matrix element measured in that case does not consist of 
a single amplitude, but rather a sum over many amplitudes to produce many final 
states: 

W v» = f d A ye iQ - y <p,s\ J"{y) \ N >< N | J"(0) | p, s > 

N J 

= J d 4 ye iQ - y <p,s\ J tM (y)r(0) \p,s> (2.2) 



The appearance of a Fourier transform with respect to the large momentum trans- 
fer Q is enough to begin motivating a small spatial separation as Q — > oo. However, 
care is still required; in the parton model, the matrix element above is reduced to 
a correlation function of quark fields. This occurs through the use of the impulse 
approximation and the "hand-bag" diagrams (Fig. 2.2) 

The expression then for W^ v is 

= J d 4 kTi[Y^ + @)Y$(k, Q)]5{k + - xp + ) 



where $ is a correlation function to measure quark fields ip(y), given by| |121 

d 4 y 

W) 



®ap(k, Q) = I <p,s\ MVMv) \p> s > elk ' y ; ( 2 - 3 ) 



I d 4 k5(k + -xp + )$ a(3 (k,p) = J^-<p,s\ Tp a {0)Mv± = V + = 0, y~)e iv ~^ \p,s> 

"Short distance" means that the quark corellations are measured practically at zero 
transverse and null plane time separation, up to scaling violations. Nevertheless, not 
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all the distances involved decrease as Q — ► oo: the remaining y~ coordinate is finite 
even at large Q. This shows that "short distance" is never automatic, and can be 
misunderstood. Even for Q 2 >> GeV 2 , short distance can fail in some kinematic 
regions, such as Feynman momentum fraction x — > 0. 

Color Transparency Involves Both Exclusives and Inclusives 

In color transparency we will be concerned with both microscopically exclusive 
and inclusive processes. Consider quasi- exclusive scattering in a nuclear target. The 
process will consist of a "hard" collision, and subsequent propagation through the 
nuclear medium. There are two distinct channels: 

1) the purely exclusive process of knocking a proton into a proton which exits the 
nucleus with only elastic interactions. 

2) the contribution of microscopic inelastic scatterings, which by further interac- 
tion, can feed back into the exclusive process. 

These two are separate theoretical subprocesses. While relating them separately 
to observables is very intricate, a proton propagating through the medium and a 
proton making two transitions are distinctly different. Many of the disputes in the 
interpretation of color transparency originate because of basic disagreements on the 
importance of these different channels. Needless to say, we will not settle those 
disputes here, but instead we hope to illustrate the merits of different viewpoints in 
a growing field. 

Let us discuss the separate channels separately, first illustrating the physics of 
elastic collisions in free space, and then bring in the complexity added by the atten- 
uating nuclear medium. 

2.2 Scattering in Free Space 

We again want to see step by step which assumptions are behind the assertions of 
short distance. We will present a fresh analogy, temporarily removing the "virtual 
photon" bias entirely from the discussion. Consider instead a pedagogical example of 
a non-relativistic atom scattering off a plane surface we can call a "partially silvered 
mirror". The setup is shown in Fig. (2.3): an atom enters carrying momentum P, 
and leaves carrying momentum P' = P + Q. 

The reason for introducing the mirror is that it is precisely not a short-distance, 
microscopic probe, but instead a crude macroscopic object. When the atom reflects 
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from the mirror its momentum component parallel to the mirror is unchanged. The 
atom's perpendicular component of momentum Pt is reversed, so we have Q = 2Pt- 
Let the mirror be arbitrarily large, in fact translationally invariant in the direction 
along its surface. Let the mirror's reflectance be very low, so that at most one of the 
constituents is reflected: the probability of more than one scattering from the mirror 
as an "external field" will be considered negligible. The mirror interacts so weakly 
with the constituents that it can be treated perturbatively, and its distortion of bound 
state wave functions (except for the scattering itself) is also totally negligible. Then 
the picture of scattering for the atom's constituents looks like Fig. (2.4). 

Observing the initial and final momentum of the system, one can not tell where 
the collision with the mirror occurred, which could have been anywhere on its sur- 
face... an obvious point due to translational invariance. Since we do not know the 
location of scattering, we cannot sensibly compare the relative position of the re- 
flected constituent with the unscattered one. The mirror therefore shows us that 
we could have an arbitrarily large momentum transfer and not necessarily probe the 
short distance separation of anything inside the atom. 

In fact, the set-up simply gives us a concrete physical picture of the interaction 
of a virtual photon carrying space-like momentum in a Q° = "Breit Frame" in- 
teracting with one of the atom's internal constituents or "quarks". Although the 
photon interacts only at a single point, the photon itself is not well localized; like 
the mirror, it has translational symmetry properties and serves only to deliver a def- 
inite momentum Q. The photon is just like a mirror for this reason. (One proviso 
should be made: in an ultrarelativistic system, only two transverse "invariant" coor- 
dinates behave non- relativistically, as in this treatment. The third, and the time-like 
"longitudinal" coordinates have a relativistic dependence that scales with the overall 
momentum of the particle and which should be treated separately.) 

Let us analyze the mirror scattering in more detail. The momentum flow of the 
internal atomic constituents is shown; a constituent's momentum is half the center 
of mass momentum (P/2) plus or minus a relative momentum (called k/2). With at 
most one constituent striking the mirror, the other constituent must fit itself into the 
outgoing atom through the fluctuations in relative momentum in the wave function. 
Let the incoming atomic wave function be 4> P (P/2+k/2, P/2— k/2), as shown. Before 
the atom strikes the mirror, this is just the ground state wave function. Because the 
collision occurs very fast, the wave function has no time to evolve to a different state. 
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This is the impulse approximation, which is a very good approximation in QCD. (It 
is so standard in QCD that it is often formalized into "hard scattering kernels" and 
"distribution amplitudes" , a notation carrying concepts equivalent to the ones being 
used here, and which will be introduced in Section 4. 

After striking the mirror, the scattered outgoing constituent has momentum P/2+ 
k/2 + Q. The full amplitude M(P, Q) for the two particles to look like an outgoing 
atom is given by the overlap back onto the ground state wave function, given by the 
Feynman rules of perturbation theory as a loop integral: 

M(P, Q) = (2tt) 3 5(P' -P-Q) j d 3 k(f)* P ,(P/2 + k/2 + Q, P/2 - k/2) 

x (p P (P/2 + k/2,P/2-k/2) (2.4) 

= J d 3 kd 3 k'd 3 k"5 3 {P' -k- k")<p* P ,(k, k")5 3 (P — k — k') 

x <f)p{k,k')5 3 {k" - k'-Q) 



In the second line we have simply inserted some delta functions for later convenience. 
We have included the overall momentum conserving delta function, often suppressed 
in the set-up of amplitudes, in order to see everything in the calculation. 

The "wave functions" are actually certain Green functions. By translational in- 
variance, for the atom carrying total momentum P, then the wave function of the 
quarks at positions x and x' should be given by exp(iP • R)(f>(r) where r = x — x',R = 
(x + x')/2 are relative and center of mass coordinates. This wave function is related 
to the momentum space wave functions above by: 

±^_ e - lk . x - lk '. x ' etP . {x+x ' )/ 2~ i ^_^ ) 

The relativistic case is exactly analogous, replacing the 3-dimensional variables by 
4-dimensional ones. To see the spatial regions of the wave function that are involved 
we insert these coordinate space representations into Eq. (2.4): 

M(P,Q) = [ ^ x< ^ x ' d 3 yd 3 y' r jz^zy^' e -ik- x -ik>- x '+ik- y +ik" y e iP-(x+x')/2-iP'-(y+y')/2 
x fa - x')$(y - y')5 3 (k" -k'-Q) 



Doing the momentum integrals we have 

k" = k'+Q; J d 3 ke ik < x -^ = (2n) 3 5 3 (x-y) ; J d 3 k'e ik '< x '-^ = (2n) 3 5 3 (x' -y') 
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This leaves some coordinate space integrals to do, which can be written 

Finally we obtain an expression we can call a "mirror form factor" F(Q): 

M(P,Q) = 8 5 3 (P' - P-Q)F{Q) ■ 

F(Q) = J d 3 r e tQ - r <f>*{r)4){r) (2.5) 

The form factor which tells us how the external, observable process of rebounding 
from the mirror depends on the internal wave function of the atom. This quantity 
would be fairly informative if the atom's constituents were hidden from us, as is the 
case with quarks in hadrons. The reader can recognize the mathematical expression 
for the mirror form factor as a familiar result, the same as the electromagnetic form 
factor for absorbing a virtual photon, because the fundamental kinematic ingredients 
are the same. 

Short Distance 

Let us can examine the expression of Eq. (2.5) for the dominance of "short 
distance". Take the momentum transfer Q to be very large. From what integration 
region in r does the form factor get its main contribution? It is not quite safe to say, 
U Q — > oo, then r — * 0, by the uncertainty principle." This is incorrect because there 
are many mathematical functions whose most rapidly varying sections are not at the 
origin! Consider, for example, the Fourier transform of a Gaussian with a node, as 
one finds for the first excited state of a harmonic oscillator (See Figure 2.5). In a 
saddle point approximation, the dominant contribution as Q — > oo to the Fourier 
transform for the form factor can be shown to come from a region near the "bumps" . 

oo 

J dx X *e- x2 l 2a2+ ^ x ~ v^QVe"^ 2 / 2 (2.6) 

— oo 

In this case the dominant integration region for the large Q behavior is r = a and 
not the r — ► region. This is an elementary result, but good to keep in mind. Note 
that the large Q dependence is exponentially damped because the wave function is 
"smooth" in this region. 
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To make a serious argument about short distance in QCD, one evidently must 
supply information about the wave function to find the quarks. We can look for 
inspiration in a more familiar gauge theory, ordinary electrodynamics. In atoms the 
wave function has a "kink" near the origin which dominates the large Q Fourier 
transform. The hydrogen atom, for example has a ground state wave function going 
like exp(— r/a). The corresponding form factor using such a wave function is given 
by 

J d 3 rexp(—iQ ■ r)exp(—r/a) = 8ira/(Q 2 + 1/a 2 ) 

~ 87ra/Q 2 , Q 2 » 1/a 2 . (2.7) 

It is an elementary exercise to verify that here the dominant region actually is 
r — > 0: the power-law decrease with Q 2 can be traced directly to the "kink" at the 
origin. The same occurs in QCD but in that case only a perturbative model of the 
wave function (valid at short distance) is known. The result in transverse momentum 
space kx is 

Mkr) ~ Xlog^(k 2 T /A 2 QCD ) , (2.8) 

where N is an x-space projection and 7at is an anomalous dimension. (See Section 
4). It is no coincidence that Eq. (2.7) and (2.8) show the same power-law decrease, 
because both QED and QCD have dimensionless couplings and vector exchange; 
perturbation theory reconstructs the short distance part of the bound state wave 
function in both cases. 

So here is the "short distance"; bin by bin in integrating over the relative coor- 
dinates r, the important region is r — ► 0, due to both large Q and the special short 
distance properties of the wave functions in the gauge theory. 

Remarks 

The mirror-photon example illustrates several points: 

* In free space it is not necessary for the system to "compress to a small size" 
to undergo a hard scattering. While self-compression is an attractive and often 
quoted semi-classical picture, to hold to such an interpretation would be a basic 
misunderstanding of the quantum mechanical impulse approximation. [] 

*The impulse approximation is sometimes mistakenly called the "frozen approximation". As we 
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In the example, the atomic wave functions are full sized, and ignorant of the fact 
that they are going to smack into the mirror. They do not time evolve significantly 
during the short time of the collision. Nor is it necessary to model the interaction 
of the system with the "mirror" in terms of complicated initial state superpositions. 
One state, the ground state, is all that is relevant in the scattering of a ground state, 
in the free space case. The same conclusion follows from inserting a complete set in 
matrix elements: 

< P + Q\ J^\P >=J2<P + Q\N >< N\J"\P >=< P + Q\J^\P > 

N 

This may seem to be a trivial point, but it is worth mentioning before considering 
the interactions with the medium. 

* The uncertainty principle is inadequate to predict much here, because it does 
not have enough detail. If we used the uncertainty principle to estimate a dominant 
region, it would generally fail. The precise statement is that an integration region as 
small as the one indicated by the uncertainty principle can be a dominant region - in 
the Hydrogen atom case, the region Ar m 1/Q turns out to dominate, as much due 
to dimensional analysis as anything else. 

* The observations made so far apply to more processes than just the scattering 
of a single photon. The scattering of two atom on each other can go through "two 
mirrors". In the center of mass frame, the picture looks like Fig. 2.6. Needless to 
say, the naive inference of "short distance" in this case is utterly wrong, because it 
ignores the need for information on the relative real space positions of the "mirrors" 
or scattering centers. It is obvious that one can have large momentum transfer 
and not have short distance in this case. We will discuss this interesting process of 
"independent scattering" in hadronic reactions in Section 4. 

* Finally, the form factor discussion is simply an example of "classic" color trans- 
parency, in which (as we will show) all of the ingredients of conventional QCD come 
together. This does not exclude other possibilities, and other motivations for small 
color singlet systems even if the short distance motivations from the impulse approx- 
imation are not evident. Certainly the validity of any proposed "small" hadronic 

have seen, the impulse approximation applies when the time scale of the scattering event is much 
faster than the time scale of the internal constituents. The "frozen approximation" is a separate 
concept applied to the relativistic diffractive limit - an issue to be discussed below. 
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system, a point like configuration, and so on, depends on the model and the con- 
text. It may be eminently reasonable and give correct predictions under the proper 
circumstances. 

Inelastics, Again 

Now we return to discuss the inelastic processes, which form the bulk of all ampli- 
tudes possible. What usually happens when an atom moving fast strikes the mirror? 
The atom explodes (Fig. 2.7). 

We seek the amplitudes to form various final states \ N > , N ^ the proton. 
Suppose we propagate into the future the set of amplitudes | s(0) > that are defined 
to strike the mirror at time t = 0. The future amplitudes are given by: 

| s (t) >= e~ iAt | N >< N(0) | s(0) > 

where H is the Hamiltonian operator. The time evolution of each energy eigenstate is 
given by | N(t) >—\ N(0) > exp(-iEt). The outgoing wave packet has a complicated 
time evolution. It expands as time progresses. We can calculate the "size" of the 
outgoing wave packet in the following way: 

< r 2 (t) > = < s(t) | r 2 | s(t) >= E < N (°) I r ' I N '(°) > e-^'-^ 

N,N' 

x < N'(0) | s(0) >< s(0) | N(0) > 

If, for example, we truncate the sum over states to two states for simplicity, then we 
find a time dependence given by 

(r 2 (t)} = A cos(E' N - E N )t + B, (2.9) 

where A and B are certain matrix elements. This indicates a characteristic time for 
the packet's expansion is set by energy differences. The uncertainty principle would 
give an acceptable estimate of this "coarse grained" flow of probability. 

However interesting, this line of reasoning is unfortunately not relevant to the 
elastic scattering already calculated. The size of the inelastically scattered wave 
packet has almost nothing to do with the elastic process. This can be illustrated 
in Figs (2.8, 2.9), showing the wave packet (in transverse separation coordinates) 
"before" and "after" the scattering from the mirror (or virtual photon). The figure 
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shows that the final state wave packet is quite broad in transverse size: it is not 
particularly small! What has been changed is relative phases among the real space 
components, so that the wave packet is broad but has become wildly oscillating. 
Because of these oscillations, what matters is the value of the overlap from the region 
near the origin. This region is an entirely different thing from the size of the full 
wave packet or its configuration, which is destined for inelasticity. If, for example, 
we mixed a symmetric ground state with a typical excited state, the excited state 
would be quite important for the bulk expansion of the inelastic wave packet. Yet, 
the node at the origin of the excited states would contribute almost no overlap to 
the dominant region at the origin for the elastic process. Finally, neither of these 
is the same concept as a point-like configuration in which the system is quantum 
mechanically prepared to exist entirely within a small region! 

We have reviewed an elementary yet significant point: the value of the ground 
state wave function at the origin is not at all the same thing as the average size of 
the inelastically scattered wave packet size. A time scale for the expansion of the 
inelastically scattered wave packet simply has little to do with the elastic scattering 
in free space. 

More Remarks 

There are a few additional remarks: 

* In many experiments the direction of the momentum transfer Q is well mea- 
sured. Suppose this is the x direction, so Q = xQ x . The "x" coordinate space 
variable conjugate to Q x is the naive region where "small" separations are measured 
in exp(— iQ x x). Nevertheless, all transverse regions that generically contribute have 
"small" dimensions. To see this, project exp(— iQ x x) into SO (2) spherical harmonics 
and find the coefficients in terms of br, the magnitude of the transverse space radial 
coordinate: 

J dcpe- tm v e - bTlQTlcosip = 27rr n J n (b T \Qt\) 

Oscillations in the Bessel functions in | Qt | &t generally wipe out long distance con- 
tributions in the magnitude br- As we will see in Section 4, the same conclusion is 
reached in the pQCD hard scattering formalism, in that case because the gluon prop- 
agator is isotropic, and selects participants from regions nearby in space regardless 
of the direction of Q. 

* Relativity allows one to repeat this type of analysis in a large set of reference 
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frames. Each frame has its own interpretation. For example, for an electromagnetic 
form factor it is convenient to use null plane coordinates (reviewed in Section 4) for the 
fast proton with momentum P^ = (P + ,P T ,P) = (P + , 0, m 2 /2P + ). In this frame, 
one can set Q» = (0, Q T = y/(-Q2),Q~ = Q 2 /2P + ). The point of this reference 
frame is that Q M is practically "all transverse" in the limit of P + — > infinity. This 
focuses attention on the transverse space coordinates, b T . The transverse coordinates 
also have the advantage of being invariant with respect to boosts along the "z" axis. 
The kinematic analysis of the previous results then applies. As a good estimate, the 
region of transverse separation b T < 1/Q is the important one to estimate overlaps 
of wave functions in hard scattering calculations. 

Because of the scaling of the "+" coordinates under a boost, it is also useful to 
express momentum fraction of quarks using the Feynman xp variables, which are 
also invariant. In Fig (2.9) we have plotted the wave function of a typical "pion" 
in a mixed (transverse space, xp) representation. This coordinate system turns out 
to be very useful and will be discussed further in Section 4. As one sees, the result 
of absorbing a virtual photon is creation of oscillations in the transverse direction 
without any decrease in the average size of a state. Compressed states postulated 
in semi-classical argument are not the result of kinematics nor calculations in QCD, 
but instead are dynamical proposals made on a case-by case basis. 

Next, one might ask, what happens in an attenuating medium? 

2.3 Attenuating Medium: Diffractive Mixing of Inelastic Pro- 
duction 

The above discussion might give the impression that inelastic processes are irrele- 
vant. This is far from obvious; in fact, many current models of color transparency 
quite reasonably emphasize the possible importance of inelastic contributions. The 
"probe" of color transparency, which is the nucleus, also interacts significantly with 
the wave packet in the process of making a measurement. In so doing, the intermedi- 
ate channels which are irrelevant in a free space exclusive experiment may feed into 
the elastic process finally measured. A further complication is that at high energies 
many states of quite different masses can contribute through diffractive excitation. 
This realization has led to the application of a pre-existing body of theoretical ex- 
pertise, developed on both sides of the iron curtain before QCD. Since the theory 
of diffractive processes has not evolved significantly in the last twenty years, the po- 
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tential for conceptual clash between the "pre-QCD" diffractive picture and the "post 
QCD" hard scattering picture is considerable. Neither picture is a substitute for the 
other; each is useful when appropriate. 

The kinematics of the diffractive limit is one where the momentum transfer to a 
spectator is small, or at least fixed, as the energy of a beam becomes large. Contribu- 
tions from each momentum transfer channel become focused kinematically into the 
forward small angle region, creating a highly coherent quantum mechanical problem. 
It is worth noting that the diffractive problem has nothing to say about the "hard 
scattering" part of color transparency experiments, if such is present. Consequently 
the impulse approximation and the concept of "short distance" , which have become 
part of QCD, are not contained in the diffractive discussion. Diffractive scattering is 
relevant to the "propagation" part of interaction with the nuclear medium before or 
after a hard scattering occurs. 



In a classic paper in 1956, Feinburg and Pomer anchuk |43|1 discussed the quan- 
tum mechanics underlying high energy perturbative calculations in electrodynamics. 
This paper is particularly useful for understanding the concept of "formation time" 



in inelastic phenomena. Some time later, Good and Walker [55] discussed the physics 
of strongly interacting diffractive production with deep quantum mechanical under- 
standing. 

Good and Walker's Polarizer 

Feinburg and Pomeranchuk (F&P) and Good and Walker (G&W) observed that 
the relevant quantity for time evolution of single particle states with mass m is not 
the energy E, but the differences in energies E for states with the same momentum 



P z . The relation for single particles is E = J (p 2 + m 2 ). In the ultrarelativistic limit, 
the energy differences go like the inverse of the big momentum: 



E x - E 2 = AE = (ml - m\) /2P Z (2.10) 

This strange, literally "asymptotic" result of relativity can be easily visualized 
by an energy-momentum diagram (Fig (2.10)). (Note that this energy difference is 
proportional to that obtained from the null plane Hamiltonian P = (m 2 + p 2 ^)/2P + .) 
It is slightly more precise to think in terms of the null plane formalism. For a given 
momentum transfer q~ to the beam, we can make an arbitrarily large invariant mass- 
squared m 2 = 2q~P + simply by having large enough P + . From this we can construct 
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a time scale, the so-called formation time, from the inverse energy. 

This concept generalizes to the case of the nuclear medium, but the free space 
formula (2.9) for it certainly does not. The problem is that the dispersion relation 
for hadrons propagating in the nucleus is unknown. The problem has the potential 
to be serious in the high energy limit, because one is asking for a small difference 
between two big, relatively uncertain numbers. We will postpone discussion of this 
question and possibly important dependence on transverse momentum, set to zero 
for the discussion, to Section 5. 

The next step in the argument is to observe that, due to relatively small effective 
energy differences, many different particle states become effectively degenerate and 
can be considered to strongly mix under an interaction. Due the relativistic effects, 
the effects of strong interactions become even stronger. 

The quantum mechanics of such a strongly mixed system is quite interesting. 
G&W make an instructive analogy with diffraction of light from a circular disk of 
polarizer. (Fig. 2.11) The degenerate states which occur for light are the two orthog- 
onal polarization states. Suppose first that an unpolarized electromagnetic wave is 
moving along the z-axis and impinges on the polarizer. Then, one component of the 
electric field is transmitted with no change, while the other component is absorbed, 
leading to a diffraction pattern with the corresponding polarization. 

Suppose, next, that the incoming wave is already polarized in the x direction. If 
the polarizer is set in the y direction, then it is a black disk for the incoming wave 
and leads to a diffraction pattern. This also is quite simple. 

Imagine finally the case where the polarizer is set at 45° to the x axis. Then the 
outgoing wave with this projection of polarization is undisturbed. The outgoing wave 
in the perpendicular polarization direction shows a black disk absorption pattern. It 
is elementary, but somewhat disturbing, to express the outgoing amplitudes again in 
the parallel and perpendicular coordinates of the original wave. In the final state, 
both x and y polarizations occur. Yet there was no y component in the initial state; 
the polarizer has produced it, by filtering away waves from the incoming x polarized 
ones! 

All this is in accord with elementary quantum mechanics. If the initial state, 
polarized at 0° to the x axis, is called | E >, and a state polarized at 45° to the x 
axis is called | E 45 >,then the polarizer corresponds to the operator | E 45 >< \ 
that covers a certain spatial region. The polarizer is diagonal in this basis. It could 
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also be written as 

| E 45 >< E 45 |= 1/2(| E > + | £90 >)(< E I + < E 90 |) 

The transition from to 90 degrees polarization is clear as the off- diagonal cross 
terms | E >< E 90 \ + \ E 9Q >< E | in the operator. Thus, conclude G&W, the 
interaction with the polarizer has, through the process of measurement, produced 
a component of the polarization that did not exist before. In the same way, they 
reason, a number of particles which are practically degenerate with the initial state 
are quantum mechanically resolved and appear to be produced by the process of 
scattering on a strongly interacting target. 

The G&W argument reminds us that quantum mechanics is needed to make a 
description of amplitudes. The strange results of quantum mechanics cannot be de- 
scribed by semi-classical or probabilistic arguments. A "classical filter" , for example, 
could exponentially attenuate the power carried by a wave. Yet a classical filter could 
not cause the creation of the new amplitudes as described above. A "quantum filter" , 
by resolving different components in different ways at the amplitude level, will appear 
to create observable phenomena merely through the process of measurement. 

The Quantum Filter in Color Transparency 

We have seen that a picture invoking bulk attenuation of hadron numbers crossing 
a nucleus would be inadequate to describe the relative propagation of different ampli- 
tudes. Instead, color transparency involves a quantum filter in the form of the nuclear 
medium. To visualize this, let us continue with the analogy of the atom scattering 
from a mirror. This time we immerse the "mirror experiment" in an atom-absorbing 
medium (Fig. 2.12). 

In the quantum medium, we will specify the dependence of the absorption on 
the internal configurations of the atom's constituents at the amplitude level. The 
theoretical object carrying this information will be called the survival propagator 
G a- Ga is the amplitude for two (or three or more) particles entering the medium at 
x and x' to be found again at positions y and y' (Fig. 2.13). This 4 (or 2 N) point 
Green function consists of 

G A (x, x'...y, y') =< A \ T[^(y')^(y) . . . ^(x')^(x)] | A > 
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where | A > is the target state and (ip*) ip are the field operators to create (annihi- 
late) the quarks. Ga is defined to be 2PI (2 particle irreducible) or 3PI depending 
on the number of quarks being propagated. We will be concerned with local color 
singlet parts of the 2N-point function suitable for hadron propagation. Flavor in- 
dices are suppressed, and we can assume we are examining the propagation of two 
(or more) unlike quarks, such as a u and d, which propagate together to finally form 
a hadron such as a n + meson. As remarked earlier, the transverse spatial separation 

— * — * 

coordinates (br = %t — x' T , b' T = yx — jf?) are the important ones in a relativistic 
treatment, while the longitudinal coordinates are considered separately. 

The survival propagator Ga is the object that contains information about inter- 
actions with the medium. Let GU(A;, k'\ £, £'; w) be the Fourier transform with respect 
to the quark momenta, with w the total amount of momentum transferred from the 
nucleus (Fig. 2. 13). The effects of interaction can be expressed in two steps. First, 
the outgoing wave function is converted to an interacting object we call 4>A{k, k'\ w): 

<P A {k, k'; w) = Jdl G A {k, k' ; I, I', w)<j>{l, I') (2.11) 

The integration over the variable 1 is indicated by the loop in the diagrams and in a 
relativistic treatment (Section 4) is 4-dimensional. The object 0a, which can be called 
a "filtered wave function," acts inside the momentum convolutions just as the original 
wave function does. In terms of its Fourier transform in the relative coordinates k — k' 
conjugate to quark separations r, we again have for the full amplitude 

M(P, Q; w) = 5{P' -P-Q-w) j dr <^(r)0(r) e -* Q - r (2.12) 

Short distance in the separations of 4>a are again selected by the oscillations in 
exp(— iQ ■ r). 

The Survival Propagator in QCD 

Obviously, the relation between <pA and the original wave functions <fi depends 
on the model for Ga- Current ideas based on QCD hold that Ga can be predicted 
perturbatively in the region where all pairwise transverse separations between the 
quarks are small compared to a Fermi. The diffractive approach, and hadronic basis 
calculations (Sections 4, 5) are different models with the same goal. For example, if 
one proposes that many intermediate states are created by the fast momentum P + Q 
plus the nuclear momentum transfers totalling w, then the model for Ga and 4>a will 
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reflect this. The outputs of these models can usually be re-expressed as models for 
4>a, so the decomposition above is rather general. 

Due to the integrations over the internal coordinates in Eq. (2.5), one sees that 
the hard scattering is tied into the filtered wave function in a non-trivial way. In 
particular, it would be an approximation of unknown validity to propose that the 
hard scattering simply multiplied the nuclear effects. 

Color transparency is obtained using the perturbative picture for Ga- As a first 
remark, the survival amplitude is "unity" minus the scattering amplitude. For the 
scattering amplitude, we practically have an almost entirely imaginary (absorptive) 
object, closely related by the optical theorem to previous calculations of the total 



cross section between "atoms". We can adapt the eikonal method of Ref. |59|] to 
calculate the scattering amplitude in pQCD. 

Consider, then, the scattering of two quarks (system 1) on a spectator atom in 
the medium (system 2) made of two quarks. Consider the 2-gluon exchange diagrams 
shown in Fig. (2.14). Flavor is ignored; since helicity is conserved by a vector interac- 
tion, spin indices can also be suppressed. Each gluon attachment is accompanied by 
a color matrix (X a )ij. Taking the singlet projection indicated, one finds conveniently 
that all matrix orderings in the same system are the same, and simply produce an 
overall factor Tr [X a \b] = CV<W 

The amplitude itself can be set up directly in coordinate space. Let B be the 
impact parameter between the two systems, with net momentum transfer w. Let the 
quarks be located at positions Xi,x' i} given by 

x x = B/2-r x 2 = -B/2-s 

x\ = B/2 + r x 2 = -B/2 + s 



In the relativistic case these are transverse vectors. Let the coordinate space kernels 
of the gluon exchange between points X{ and Xj be denoted V(xi — Xj). There are 16 
possible combinations of the pairwise interactions, of the form V(x\ — X2)V(x' 1 — x' 2 ), 
etc. These combinations, denoted K(xi,Xj), will be listed momentarily. Then the 
amplitude is 



G A (x,x; y,y f ; w) = 5(x-y)5(x'-y') 



1- / dBe 



—iwB 



ds | cp(s) | 2 K(xi, x'j) 



(2.13) 
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Now we consider the kernel K. Each exchange between opposite color charges comes 
with a relative "minus" compared to the exchanges between like charges. As the 
reader can check, the total of all the combinations of the form 

Kix^x'j) = V(x 1 -x 2 )[V{x 1 -x 2 )-V{x 1 -x' 2 ) + V(x' 1 -x' 2 )-V(x' 1 -x 2 )} 
- V(x\ - x 2 )[V(x 1 - x 2 ) - V(xi - x' 2 ) + V(x[ - x' 2 ) - V(x[ - x 2 )] 
+ ... (2.14) 



can also be written 

K( Xi , zj) = [V(xi - x 2 ) - V{x x - x' 2 ) + V(x[ - x' 2 ) - V(x[ - x 2 )} 2 (2.15) 

One sees the results of coherence in the "minus" signs, which represent destructive 
interference between emissions from opposite color charges. We obtain the perturba- 
tive expression for Ga scattering from one atom spectator: 

G A (x, x', y, y', w) = 5{x - y)S(x' -y') I- J dBe~ iwB J ds\ <p(s) | 2 

x [V{ Xl - x 2 ) - V{xi - x' 2 ) + V{x[ - x' 2 )V(x[ - x 2 )f 

(2.16) 

We are interested in this in the short distance region selected by exp(— iQr). 
A power series approximation to the interaction part of Ga in this variable, with 
r =| xi — x'i |, yields 

K( Xi , 4) = [r • V(V( Xl - x 2 ) - V( Xl - x' 2 ))} 2 (2.17) 

This goes like r 2 , showing the characteristic "geometrical" behavior of dipole-dipole 
scattering (and the associated dipole moments in the expansion of V). 

An exactly similar result follows in the relativistic treatment, if one expresses 
space variables using B s , the transverse separation between the center of P + of the 
right-moving quarks, and center of P~ of left moving spectator, in the center of 



momentum frame. One recovers in that case the results of Ref . [|124j| . Although the 
perturbative kernels seem to suffer from infrared divergences, the calculation is finite. 
For example, the transverse Fourier transform of a massive gluon is given by 

e ik T (x T -x T ) 



V(x -x) = / d 2 h 



k 2 + m 2 
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- I x — x' I m Ki{\ x — x' | m) , (2.18) 



where K\ is a Bessel function of imaginary argument, and where m is a gluon mass. 
This is badly behaved as m — > 0. Yet, the region of transverse position sensitive to 
this is erased by coherence, and replaced by scales characteristic of the bound state, 
so the massless gluon limit is "safe" . K: 



The main result is that the effects of interaction in the perturbative treatment 
are a factor of the spatial variable "r 2 ", at least as r — > 0. Because this is the region 
selected, and it also vanishes in the region of large Q, we have color transparency as 
a rather general result in this limit: 



lim / dre- l ^ r cp* A (r)(f(r) = / dr e~ lQr V * {r) V (r) (2.19) 

Q^oo J J 

Color transparency in pQCD is thus a feature of kinematics, coherence, and the 
overlaps of wave functions, rather that of superposition over states; a "one hadron" 
universe is clearly sufficient for it to occur. The fact that a completely perturbative 
treatment of color transparency is self consistent has not been widely appreciated. 
Because the interaction probability is small, even the presence of many targets (a 
complicated nucleus) does not disturb the result asymptotically. Detailed models 
of the nucleus as a target could, in principle, lead to straightforwardly calculable 
pertubative results. 
A Model Independent Framework 

The framework here can accomodate many models; we have devoted the entire 
Section 5 to "models of propagation" . Perhaps the most controversial item at finite Q 
is the effects of the nuclear momentum transfer w. In the hadronic basis treatment, 
many physical states are used to create a superposition which moves through the 
medium, created by the large P + and moderate w- values from the nuclear target. 
This can be expressed as a superposition of states model for Ga- In addition, inte- 
gration over the energy-like variable of w is thought to create a time-scale (or many 
time scales) which might give complicated dependence on the target length. These 
ideas will be discussed in detail in Section 5, while general remarks on the time scale 
issue will be given soon. There seems to be no disagreement that the formalism of 
Ref. ||124|| is the correct description at very large values of Q. There is, of course, 



plenty of controversy whether this can apply at current laboratory values. f42|, |125 | 
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Before closing this subsection we point out an amusing point. Among all bound 
states, the value of the wave function at the origin is the largest in the lowest lying 
states, and generally decreases monotonically as the states number increases, at least 
for atoms. 

Since the quantum filter of G a has caused the region near the origin to propagate 
most favorably, it is bound to dominate the measurements. This again is called 
"nuclear filtering" . If, for example, the value of Q were too small to make the Fourier 
motivated approximation of short distance very good, then the region near the origin 



might anyway be selected by sheer "survival of the smallest". [124] Transitions to 
excited states involving less wave function at the origin and more wave function in 
the peripheral regions getting "filtered away" would be suppressed. Thus, large A and 
large Q provide two complementary but independent motivations for short distance. 

The Time Scale Problem 

So far we have not discussed whether the path taken in the attenuating medium 
is "long" or "short". If the attenuating medium is shallow, then it hardly matters 
what the attenuating interaction is, and color transparency may require cleverness 
to measure. (As we will see, electroproduction has to contend with this problem.) 
The path length is this context is long or short compared to some characteristic 
scale. What scale? Are we to think of the size of the nucleus in units of some hard 
scattering, short distance scale, an interaction scale, a mass scale, a scale modified 
by a Lorentz boost, or what? 

Closely associated with this problem is the "time scale" problem, addressing 
with non-relativistic concepts the dynamics of the time evolution of a hadronic wave 
packet. In most models, if the wave packet expands rapidly, it is assumed to lose the 
capacity for exhibiting color transparency. In the same models, the observable energy 
dependence of the processes, which of course includes the effects of the Lorentz boost, 
is often interpreted rather directly as a measure of the time evolving dynamics. 

The details of the initial state become crucial in this discussion. In many treat- 
ments, the initial state is set up as a boundary condition. There are two extremes: 
a geometrically small wave packet might be formed at time t = 0. This proposal is 
ad-hoc, but worth considering. On the other hand, in the impulse approximation, a 
hadron-sized wave packet is disturbed by the addition of large momentum transfer 
oscillations at time t = 0. This is what happens in the theory of hard scattering 
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QCD. 

Because color transparency may exist in many circumstances, and the clear-cut 
kinematic conditions of hard scattering theory may often not be fulfilled experimen- 
tally, we will consider each of the above extreme cases here. For example, in charmed 
quark pair production the initial state is unknown and a small wave packet can be 
invoked as an idea to be explored. To complicate things further, the time evolution in 
an attenuating nuclear medium is actually unknown, because the energy eigenstates 
in interaction with the medium are unknown. Thus, many disparate pictures are 
possible. 

We now turn to more specific remarks. Many authors have proposed the initial 
state is "compressed" or quantum mechanically prepared in some way to exist in the 
region r < l/Q. From the uncertainty principle, the dominant momentum space 
region of such a wave packet has k = Q, leading to a rapid expansion in size. On 
the other hand, in the impulse approximation we have emphasized here, the final 
state wave function is exp(i(P + Q))<f>(k + Q). The relative coordinate momentum 
space wavefunction is centered at P + Q, a "big" number. Because this is just a 
translation, the fluctuation in momentum is still set by the hadronic size. To see 
this, let Ak 2 =< k 2 > - < k > 2 . Then 

Ak 2 = J dk^*{k + Q)k 2 v {k + Q) - J dk v *{k + Q)k v {k + Q) 2 

dk(p*(k)k 2 ip(k) (2.20) 



It follows that the case of the impulse approximation has a different time evolu- 
tion, because the fluctuations are not large. (Fig. 2.15). We will present a simple 
calculation to verify this. 

Of interest are the short time dynamics, for which we can use the power series 
solution of the operators in the Heisenburg picture as a function of the time t: 

f(t) = r(0) + it[H, f (0)] + (it) 2 [H, [H, f (0)]]/2! + . . . (2.21) 

The Hamiltonian in such a non-relativistic framework (which is adequate to illustrate 
our point) is written as H = p 2 /2m + V(r). A short calculation shows that the time 
dependence of the relative coordinate-squared r 2 is given by: 

f- 2 (t) = r 2 (0) + (t/m)[p(0)f(0) + r(0)p(0)] + t 2 [p 2 (0) - (l/2m)dV/dr] + ... (2.22) 
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Now consider the average r 2 as the struck wave function <f> s (r) time-evolves to- 
wards the final state wave function <f)f(r), measured by < r 2 j(t) >, where 

< r 2 sf (t) > = X - f rf 3 r0}(r)e +l ^f 2 e-^V s (r)/F + c.c. 

= -jd*r<t>* f {r)r 2 {t)<l> s {r)/F + c.c.- (2.23) 
F = j d 3 r0}(r)0 s (r) , 

where c.c. stands for complex conjugate. (We keep the real part; the reader may 
verify that the imaginary part works the same way.) The reason for using < r^(t) > 
rather than a diagonal expectation value in the struck state is to probe the regions 
involved in the "particular probability" to go from the struck state and reach the 
final state, and not the "coarse grained probability" for the whole wave packet to 
expand to any old inelastic state whatsoever. Of course, < r 2 s ^{t) > is also an average 
r 2 measured in the form factor, which we have already seen is the relevant thing for 
interaction with the nucleus via the survival propagator Ga- 

Because we want representative short distance behavior, the final state wave func- 
tion will be been chosen to be Coulomb-like, as opposed to other possibilities such as 
a Gaussian: 

(f) f (r) = e~ r/a (2.24) 

where a is a bound state distance scale (the Bohr radius for an atom). The two 
pictures of the wave function <p s just after being struck are: 

4>sij) = e~ r ' Q ! compressed case ; (2.25) 

(f) s (r) = e-r/o-tQ-r impulse case ; (2.26) 

We have used the fact that in the pQCD kinematic analysis the wave function — > 
0exp(— iQ-r) after being struck, while in the compressed case it is assumed to reside 
in a region r < 1/Q . 

The details of the calculation depend on the choice of the Hamiltonian, but the 
main features are easy to anticipate. Suppose we use a harmonic oscillator Hamilto- 
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nian, with V(r) = mw r /2. Upon doing some integrals, we find for the two cases: 

2 . . 12 1, Q 6«7q . 

< r s fW >= 77m — 7TT7 + * ( 9 ~ n ^ i m / ^ + • • • compressed case 

s/w (| Q | +l/a) 2 v am 2 (| Q | +l/a) 2 ' 

(2.27) 



2/N 8Q 2 -24/a 2 2 / l 4Q 2 u> 2 -12«; 2 /a 2 \ 

< r s f W >= 7777- 71 ono + * — t~7 7777- 77 9N9 — + • • • impulse case 

; (Q 2 + 4/a 2 ) 2 \m 2 a 2 (Q 2 + 4/a 2 ) 2 / 

(2.28) 

In both cases short distance is being selected, in one (impulse case) via the Fourier 
oscillations, and the other (compressed case) by simply imposing it. However, the 
relevant rates of expansion of the "compressed" wave packet and the "hard struck" 
wave packet are strikingly different. The expansion of the compressed wave packet 
is set by the initial size and goes like Q. The constant expansion rate of the wave 
packet in the impulse case is due to the hadronic time scale.Q 

Relativity makes the situation even more interesting. In a typical high energy 
hard scattering observation, the momentum transfer Q 2 goes like the beam energy E 
(for example, Q 2 = const s = const 2mE, in fixed angle hadron scattering; Q 2 = 2mE 
where E = the photon energy v in electroproduction where m is a mass scale). The 
Lorentz boost to the lab frame expands every time scale in the struck hadron rest 
frame by a factor of E/m. In the lab frame we will observe the above formulas with 
Q 2 -> 2mE, t -> tm/E: 

< r 2 sf (t) > LA B= Q/(mE) + V2mVEt 2 /(E 2 a) + ... compressed case (2.29) 



< r 2 f (t) > LA B= 4:/(mE) + t 2 /{a 2 E 2 ) + ... impulse case (2.30) 

In the compressed case one does not gain as much by going to higher energy, because 

the system simply expands more rapidly to compensate. This has caused many 

difficulties for this proposal, and led some authors to conclude that large Q 2 leads to 

"rapid expansion." But, in the hard scattering case from the impulse approximation 

kinematics, there is an advantage to high energy, because every effective time scale - 

^The reader may wish to repeat a similar calculation for the time evolution of f 2 (t) in the diagonal 
case / d 3 rip*r 2 (p s , which measures the coarse-grained probability. 
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even including the increased hardness of the scattering - increases in the relativistic 
high energy limit. 

Perhaps more effort has gone into studying the time scale issue in color trans- 
parency than any other aspect. The example shows not only that no unique "time 
scale" exists; even the functional dependence on energy depends on the assumptions. 
The resolution is subtle and depends on the experiment. It is reasonable to suppose 
that the true situation might be intermediate between the two pictures with quantum 
mechanical superposition of the two pictures playing a role. Nevertheless, the rapid 
expansion dynamics of a "compressed state" force it to be the most difficult in which 
to observe color transparency. By the same token, the dynamics favors selection of 
hard scattering impulse amplitudes because these expand at a slower rate. This is 
another case where the quantum mechanics of nuclear filtering are important, be- 
cause even if hard scattering conditions seem to be fulfilled only marginally, there is 
every reason to believe that selection by triggering on elastic survivors should make 
this contribution relatively more important. 

Color Transparency in the Hadronic Basis 

Much of the discussion of color transparency in the hadronic basis focuses on the 
role of inelastic intermediate states. Since there are no "quark relative coordinates" 
in this basis, it is not clear how to interpret the perturbative prediction for Ga- Most 
workers have tried to understand how color transparency might happen by creating 
various assumptions or parameterizations of Ga- 

The role of the hadronic basis and inelastic channels can be explored [pTj] by in- 
serting a complete set of states in the s-channel cut of Ga- Note that a complete set 
means all physical processes, namely all numbers of particles of all types, including 
resonances, continuum, all possible momenta, the works. One can write an s- channel 
dispersion relation (Fig. 2.16): 



where disc is the s-channel discontinuity (aka "imaginary part"). Such expressions 
should be model independent and exact. State by state in s', we can imagine doing 
the integrals over the internal coordinates of the constituents for each on-shell final 
state. The resulting expression makes no references to quark coordinates, and should 
permit one to survey the contributions of the various hadronic state channels. 
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Recall the discussion of the wave function after the hard scattering (Eq. (2.5)). 
The wave function of the struck state is just the original one multiplied by exp(— iQ-r). 
The contribution of this state onto each final state < N\ is given by the overlap 
< N\exp(—iQ ■ r)\P >. This, of course, is the transition form factor to the particular 
state: 

< P + Q\J' l \P >=J2<P + Q\N >< NIJ^P > 

N 

where the field theory operator J M was just replaced by e~ l ®' r in non-relativistic one 
photon exchange. The set of channels which can be inserted is potentially very large: 
for a typical momentum transfer of 300 MeV, and with a beam momentum of 10 
GeV, one can involve all states with the proper quantum numbers up to an invariant 
mass of 6 GeV 2 ! From deeply inelastic scattering we know that the continuum states 
make a large contribution. (Conversely, saturation by a few resonances is a poor 
approximation.) With so many states the possibility of difractive rescattering back 
into the original state may be considerable, and the hoped-for understanding of color 
transparency in terms of hadrons might occur. 

Unfortunately, there is no universally agreed proceedure for making the calcula- 
tion. Successful models exist, but almost all of them make the "compressed state" 
boundary condition as an initial assumption, leading to a rapid (non-pQCD-like) ex- 
pansion rate. Because this is an initial condition in the hadronic basis treatment, 
subsequent calculations do not determine whether it is "right" or "wrong" . An excep- 
tion to the compressed initial state assumption is the recent work of Markovoz and 
Miller We will discuss various models in more detail in the section on "models 
of propagation" . 

Few Particle Unitarity 

We finish the section with an observation. Ordering the set of intermediate states 
by increasing energy, the first state to contribute is the proton itself. In the quark 
picture this "diagonal" contribution is sufficient for predicting a significant part of 
the color transparency. No other states are needed, because the wave function for 
small transverse separation propagates almost freely, and has a comparatively strong 
overlap with the final state. It is a very detailed question whether the single state 
would be a "good" approximation. Nevertheless, the quarks say that there is no strict 
need for mixing among hadronic states to have color transparency. We can drop all 
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the other states, truncating the universe to the diagonal contribution, and still see 
the effect! We need only to refer back to the impulse approximation in the quark 
basis Eq. (2.5) to verify this. 

This is a very peculiar thing, because if we had independently measured the atten- 
uation of Hydrogen atoms in a medium without a hard scattering mirror experiment, 
we might not anticipate it. But it is actually simple. The total attenuation cross 
section of an atom not undoing hard scattering would be set by its overall size, of 
order one angstrom. For the total cross section experiment the sum over all regions 
of tranverse separation is needed; the region of the largest size in the wave function 
is the most important. The total cross section for attenuation is a square Angstrom. 
This is a measurement based on one, "center of mass" coordinate as determining the 
dynamics. Compare this to the hard scattering experiment followed by attenuation. 
This depends sensitively on the transverse region with the smallest spatial size, which 
is entirely different, and which probes the internal coordinates. Should we attribute 
this to the role of inelastic processes? No. The relation is simply unitarity in the 
single atom sector. (Fig. 2.17). 

Because an atom is made of only two constituents (a very good approximation, 
for atoms), then we will get the total cross section by integrating over all possible 2 
particle final states of the constituents. This relation is 



where G\ is the one-state contribution. Since disc G\{r = r') is a positive definite 
quantity, the integral over it is certainly larger than the contribution from the short 
distance region of r' < 1/Q. Thus it is rigorous that a process selecting short distance 
will make a comparatively small contribution to a tot . The one particle contribution 
will kinematically give color transparency (Fig. 2.17). 

The same conclusion is very difficult to obtain by simply considering the "atomic 
center of mass" (hadronic) basis alone. This is because the matrix elements needed are 
new diagonal contributions not readily accessible from bulk "center of mass" hadron 
cross section measurements. Moreover, one can easily see that the mixing of the 
different channels might be entirely irrelevant. Imagine a hypothetical case where 
the spacing between excited states of an atom (the different "hadrons") were made 
arbitrarily large in energy. For the Hydrogen aton this can be done by taking the 
Rydberg energy me 4 /2ft 2 — > oo. With an arbitrarily large separation of levels, 
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mixing of levels will not occur. At the same time, we can arrange not to have the 
wave function at the origin change at all. We simply take the limit in a way so as to 
keep fixed the Bohr radius a = h 2 /me 2 . The joint limit is arranged by taking e 2 — > 
and m — > oo with m 2 e 2 fixed. According to our earlier calculation, the scattering 
of the atom via the elastic "ground state" channel not change in this limit, while 
the mixing of off-diagonal channels would decrease substantially. We can even take 
the opposite limit: let the separation between levels be reduced as much as desired, 
while still keeping the wave function at the origin fixed. In this case a host of levels 
become available for mixing. They do not matter; same story, same result, for the 
diagonal contribution. Yet, by rescattering, these might add probability separately 
in this case. 

Some arguments on color transparency miss this point completely. The "quantum 
filter" of the nuclear medium can operate this way. It occurs because the different the- 
ories based on quarks and hadrons are not the same; quarks act as if they have more 
degrees of freedom than hadrons. But as a separate issue, the mixing of states might 
nevertheless be important when we do the experiment in an attenuating medium. The 
diagonal contribution does not exhaust all possible contributions to the "s-channel 
cut" of Ga by any means - there are many other channels. It might be possible to 
create an observable which can separate diagonal from off-diagonal transitions. In 
the absence of any clear separation, the relative strengths of the different processes 
constitute a research problem. In Section 5, we review the various models. 



37 



4 Figure Captions 

.1 The elastic form factor. Assumptions are required to show that short distance 
is probed by large Q. 

.2 The handbag diagram for deeply inelastic scattering. The sum over the final 
state (dashed line) runs over all number of all kinds of constituents, including 
both the struck quark and undetected final states of all kinds. 

.3 An atom scattering off a weakly interacting mirror (left) receives momentum 
in the same way it absorbs a virtual photon with spacelike momentum transfer 
(right). 

.4 Kinematics for scattering with momentum transfer Q from a mirror, or a virtual 
photon. The wave functions are indicated by blobs. 

.5 A wave function's large Q Fourier transform may receive its dominant con- 
tributions (arrows) from regions not representing small separation (left) or the 
region of small separations (right) depending on where the most rapidly varying 
region occurs. 

.6 Independent scattering from two "mirrors". Assuming large momentum trans- 
fer, short distance is nevertheless not guaranteed; in particular, the separation 
of the two mirrors in coordinate space can be large. The same phenomenon can 
occur when quarks exchange gluons. 

.7 An inelastic event in the quark basis (left) and in the hadronic basis (right). 
Sums over all final states in the quark basis are done using unitarity. In the 
hadronic basis, one can attempt to model all the terms individually. 

.8 (a) The relative coordinate wave function of a full-sized Hydrogen atom, an 
exponential, (b) The real part of the same function multiplied by exp(— iQx), 
which is the result of absobing a virtual photon. Plot is to scale with Q = 
10. The overlap of this wave function onto the ground state gives the form 
factor; the result is dominated by Ax — 1/Q. (c) A hypothetical wave function 
(arbitrary normalization) described by exp(— 10a; 2 ) from the "compressed intial 
state" hypothesis, in which a wave function is quantum mechanically prepared 
to exist in a region Ax = 1/Q. 
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2.9 Plots of a typical wave function for quarks in a pion in a mixed representation, as 
a function of the longitudinal momentum fraction xf and the transverse space 
coordinate br- (a) The wave function is exp(— — x), where br is the 
transverse coordinate in units of Fermis. (b) Real part of same wave function 
after absorbing a photon in the bx direction with momentum Q = 1 GeV. 
The width of the wave function in bx is unchanged, but the photon produces 
oscillations which select the region bx < 1/5 Fm from the region near the origin. 

2.10 Even large mass differences become small energy differences for states moving 
close to the speed of light. To predict the small energy differences used in 
the hadronic basis approach, one must know the dispersion relation with great 
accuracy. 

2.11 An incoming electromagnetic wave polarized in the x direction hits a polarizer 
oriented at 45 degrees. The outgoing wave is measured to have both x and y 
polarizations, illustrating the creation of new degrees of freedom by a quantum 
filter. 

2.12 The mirror experiment is placed in an absorbing medium, which creates a 
quantum filter. Regions of the wave function with small transverse separation 
propagate well, while large separation regions are more attenuated. The sur- 
vival propagator Ga represents this interaction. Momentum Q, received from 
the hard interaction, is distinguished from momentum w transferred from the 
medium. 

2.13 (a) the survival propagator Ga for a quark anti-quark, (b) the survival propa- 
gator Ga for 3 quarks. 

2.14 Diagramatic expansion of Ga- The spectator atoms (system 2) have wave 
functions indicated by blobs; a propagating color singlet (system 1) is to be 
attached to a hard scattering on one side, and a wave function for the final 
state on the other. Other spectators are not shown. Space coordinates used in 
the calculation are indicated; the separations between quarks in system 1 and 
2 are r and s, respectively. 

2.15 Momentum space models for the struck wave function, (a) A full-sized wave 
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function <p(k) before being struck. It is centered at k — 0, with typical hadronic 
fluctuations of relative momentum of order 300 MeV. (b) The same as (a) after 
receiving large momentum Q, which causes a translation. The relative mo- 
mentum is order Q, but fluctuations remain of order 300 MeV. (c) Momentum 
space picture of a "compressed" state, which is centered at large Q and also 
has a large fluctuation of order Q. Large fluctuations cause rapid expansion 
not occuring in the impulse approximation. 

2.16 The survival propagator Ga can be written using a sum over on shell inter- 
mediate states N, an s-channel dispersion relation, for each fixed momentum 
transfer w to the nucleus. If the intermediate states are very massive, a few 
resonances will not be a good approximation. 

2.17 Unitarity and color transparency, (a) A single hadron contribution to the imagi- 
nary part (s-channel discontinuity) of the survival propagator (G$) for a quark- 
antiquark. The momentum transfers and relative separation coordinates r of 
the quark antiquark are held fixed; the region r — > is of interest, (b) Unstitch- 
ing diagram (a) and reassembling it backwards, the integral over all r will give 
the scattering amplitude a to t in this channel. Since the integrand is positive 
definite, the integral over any finite region of r of disc (G^\r)) is always less 
than atot, final state hadron by hadron. 
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3 Experimental Definitions and Criteria 



3.1 Overview 

In this section we review some of the ideas which have been put forward for finding 
color transparency and review them critically. This subject is bubbling with ideas, 
and it is important at present to have a liberal view of what might be an interesting 
experiment. 

A basic requirement for color transparency is that the color flow should be con- 
trolled. For example, in exclusive processes we believe we can account for the flow of 
color charge by studying Feynman diagrams of the microscopic processes. When a 
color singlet flows coherently through a reaction, as occurs when we trace the quark 
lines in an elastic scattering, then one can argue that all that is needed is short 
distance dominance and the reduced attenuation signal of color transparency should 
be observed. Insisting on large Q 2 may not be necessary because attenuation in 
the 'quantum filter' may be able to create short distance (or small color separation) 
conditions on its own. 

In the earliest work, the conditions of color transparency were fairly tight: to 
enforce exclusivity, one desired sufficiently restrictive cuts so that no pions would be 
produced beyond those discussed in the reaction. As the field has evolved, many 
experiments (e.g. at FNAL energies) are being examined in which this restriction 
cannot be enforced. This seems to be healthy development to the extent that some- 
thing can be discovered in kinematic regions where theory may not be well developed. 
The search for color transparency is primarily an experimental task; like gold, it is 
going to be 'where you find it'. We will begin by organizing the discussion under 
common questions, and then turn to reviewing suggestions for experiments. 

What cuts are required to define an acceptable exclusive reaction? 

In the first BNL experiment on color transparency in pp scattering [2S|, only 



processes which satisfied the nominal conditions of elastic scattering within an energy 
resolution of about 20 MeV were included. This cut, of course, implied that there 
would be no pions in the final state on a target at rest. Since the targets in the 
nucleus are not at rest due to Fermi motion, the experiment also included a veto 
detector designed to reject cases where pions were emitted by the final state nucleus. 
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Although such restrictive definition were necessary for a pioneering experiment, we 
argue that the experimental cuts can be made less restrictive and still represent 
coherent color flow. At Fermilab energies, for example, such restrictive cuts would 
require unrealistically precise determination of the energies of all the particles. A 
small error in the measurement of energy Ap of the incident particle would lead to a 
large uncertainty in the squared invariant mass, AM 2 rs Ap^, where s is the center 
of mass energy squared. Therefore at Fermilab energies it would be very difficult to 
require that no additional particles are produced unless specific detectors are added. 

The upgraded BNL experiment will measure all the momenta of incoming and 
outgoing particles, allowing for an overdetermined kinematic situation. This is a 
significant experimental accomplishment, and may also make possible many unusual 
experimental cuts (such as direct observation of the system "center of mass"). The 
interpretation of the data can be somewhat ambiguous, however, because the sys- 
tem being measured interacts with the nucleus. Typical momentum transfers in a 
strong interaction scattering of a few hundred MeV are comparable to typical Fermi 
momenta, so some modeling of data which 'measures' Fermi momenta directly is 
necessarily involved. Depending on the observable- we will discuss several in the fol- 
lowing sections- a combination of modeling ordinary nuclear effects and clever choices 
of what to measure certainly extends the interesting allowable energy resolution into 
the 100's of MeV range. The opposite extreme, of requiring very high resolution 
so that individual nuclear shell structure can be picked out, may be interesting but 
seems at the present time to be "overkill" . Finally, if there is some information avail- 
able from production of pions, then such processes need not necessarily be dismissed; 
theory may get to the level of understanding the color flows. 

In high energy diffractive p production at Fermilab |39], f|(], [l[ there is an indication 



that color transparency might be observed. The kinematics of this case are such that 
final states with huge invariant masses may likely sneak through the cuts imposed 
on the detection process. This again is not necessarily catastrophic, because it may 
be perfectly OK for a nucleus to "explode" after the event has occurred. One should 
be cautious, however. In some cases, such as hard scattering, the criteria of no 
pion production seems feasible and straightforward. In other cases, and until various 
theories reach a high enough level of accountability, it is rather difficult to claim that 
one is making a definitive test. 

There are two characteristic limits: hard scattering and diffractive. Generally, 
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in hard scattering all scales are large; in the diffractive limit, s >> other scales. 
Perturbative QCD generally applies in the hard scattering limit, its application in 
the diffractive limit is a research problem. 



3.2 Hard Scattering Experiments 

There currently exists two different hard-scattering experiments that may show the 
phenomena of color transparency. The first type of experiments involve wide an- 
gle hard hadron-nucleus scattering, the second lepton-nucleus scattering. The BNL 
pA experiment and the SLAC eA experiment fall under these categories. The phe- 
nomenon of color transparency apparently has been already observed in at least one 
of these types of experiments but requires further confirmation. 

The hard scattering experiments require s » 1 GeV 2 and Q 2 >> 1 GeV 2 with 
s/Q 2 of order 1, where yfs is the center of mass energy and Q 2 is the momentum 
transfer. The impulse approximation applies to these experiments and therefore it 
should be possible to factorize the hard and soft scattering kernels for such processes. 

BNL: a brief review 

The pioneering color transparency experiment involving hard scattering was done 



at BNL by the group of Carroll et al. [28] . A proton beam was used to study 



pA — > p'p" (A — 1), (Fig. (3.1)). Data was taken simultaneously on 6 targets, namely 
1 if(in the form of polyethylene), 7 Li, 12 C, 27 Al, 65 Cu and 198 Pb. The beam energies 
were 6, 10, and 12 GeV, although the data was reported at 12 GeV only for Al and 
Cu. The data was measured at a quasi-elastic point of 90° cm scattering. (At this 
point the data depend only on one scale Q 2 = —t = —u ~ s/2.) The spectrometers 
were set at the correct angle and energy for an elastic collision, within the resolution 
of the instruments, and a veto was placed on the nuclear targets to exclude events 
where a large amount of inelastic energy was deposited. In practice any events in 
which pions might have been produced were substantially reduced. Nevertheless one 
could not determine what excited nuclear state was left behind. The four-momentum 
of one of the final state protons and the track direction of the second proton were 
measured. The Fermi momentum of the target proton was then deduced from the 
kinematic information available by assuming that the event was a two body elastic 
collsion. Fig. (3.2) shows the projection of the nucleon momentum distribution along 
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the direction perpendicular to the two body nucleon-nucleon scattering plane 0, |6^] . 
The incident momentum is 6 GeV for Li, C, Al, Cu and Pb and 10 GeV for Al. 

Fermi motion of the target nucleon inside the nucleus makes the experimental 
extraction of transparency rather complicated. It introduces a separate kinematic 
variable which has to be deduced from the available experimental information and 
assumptions involving elastic kinematics. To give an idea of the magnitude of the 
effect we compute the shift in cm energy-squared s due to Fermi motion of the target 
nucleon in BNL kinematics. Let p' be the four momentum of a proton at rest. 
Including Fermi momentum kp, the initial momentum of the struck proton is p' + kp. 
The cm energy-squared for a collision with a beam proton carrying momentum p b is 
given by, 

s = (p b + p' + k F ) 2 = 2m 2 + 2E b m - 2p b k F . 

where m is the mass of the proton. Taking p b equal to 10 GeV and kp equal to 0.3 
GeV the change in s from its value if the target nucleus were at rest is about 3 GeV 2 . 
A small change in the value of s can make a big difference in the cross section. An 
experimental cure for this problem is to measure enough variables to determine the 
Fermi momentum of each struck proton, which was done in the BNL experiment. 
Thus Figs. (3.3) and (3.5) show data for values of the effective beam energy which 
range around the values actually used, as the Fermi motion effects were incorporated 
by the experimentalists. A weakness in the procedure is that one must know the 
Fermi momentum distribution: the BNL experiment used its own measured out-of- 
plane momentum distributions as well as standard models for this. Needless to say, 
the role of the Fermi motion may be more subtle than this and it has received much 
theoretical attention. 

The experiment introduced and reported a ratio T(Q 2 , A) called the "transparency 
ratio" , defined by 

T (da/dt)(pA^p'p"(A-l) 
Z(da /dt)(pp — > p'p") 

The experimentally measured values of T as a function of incident momentum are 
shown in fig. (3.3). The reason for reporting the transparency ratio is that certain 
systematic experimental uncertainties in the data cancel out. 

It is often loosely said that the free space pp — > pp scattering in the denominator 
scales like s -10 . The transparency ratio, then, conveniently takes out a very rapidly 
varying function of the energy. The factor of 1/Z represents the expectation that the 
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hard scattering is incoherent: the only interference that should occur at large momen- 
tum transfer comes from one proton being struck at a time. The transparency ratio 
is an observable quantity which takes into account a basic idea of factorization, which 
roughly means that the hard scattering effects multiply the soft nuclear interaction 
effects. However, the transparency ratio is tricky to interpret and other observables 
can be constructed. This is discussed in more detail later. 

SLA C NE1 8: a brief review 

The second experiment involving hard scattering was performed at SLAC and 
involved quasi-elastic electron scattering on nuclear targets as shown in fig. (3.4) [J94"[ 
109|| . In this case five nuclear targets 1 H, 2 D, 12 C, 56 Fe and 197 Au were used. The 
momentum transfers reported were Q 2 = 1,3,5 and 6.8 (GeV/c) 2 . The experimen- 
talists reported a transparency ratio which is defined slightly differently than in Eq 
(3.1). We denote this transparency ratio as T'. It was defined as 

rpl _ ^TZ Ndata 2^ 

J2k Npwia 

where N refers to the number of events in the kinematic region TZ and PWIA refers to 
the theoretical calculation using the plane wave impulse approximation. The region 
TZ for the SLAC experiment was chosen to be TZ : [—25 < E m < lOOMeV; < 
p m < 250MeV/c], where E m and p m are the missing energy and momentum due 
to the Fermi motion of the target nucleon inside the nucleus. In the plane wave 
impulse approximation the proton is ejected from the nucleus without final state 
interactions by exchange of a virtual photon with the incident electron. It therefore 
corresponds to treating the target protons as being free. The experimental results 
for the transparency ratio T are shown in figs. (3.5) and (3.6). Fig. (3.6) shows the 
comparison of the results for Carbon nuclei with different theoretical predictions. The 
transparency ratio was found to be roughly flat with a very slow increase with energy 
for all the nuclei studied by NE18. The Fermi momentum distribution extracted by 
the NE18 group is shown in fig. (3.7). 

phenomenological remarks 

There have been many theoretical studies to predict or explain the results of 
color transparency experiments involving hard scattering. The rise and fall of the 
transparency ratio with increasing energy, seen by the BNL experiment, apparently 
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contradicted the expectations of color transparency as well as of the Glauber model 



2^, IT, |9Cj. After some inital confusion it became clear that some of the observed 
structure was caused by oscillations with energy in the denominater of Eq. (3.1). 
This behavior of the transparency ratio is understandable [E3, [TT5] if there are long 
distance contributions to the exclusive wide angle pp elastic scattering in free space. 
In the Brodsky-de Teramond [J23| (BdT) model the long distance contributions are 
postulated to arise due to several charmed dibaryon resonances. These resonant 
contributions are filtered by the nuclear medium leading to the observed bump in the 
transparency ratio. Ralston and Pire [ 119 1 instead argued that exclusive processes 
generically get considerable contributions from long distance components, an example 



being given by the independent scattering process [88]. Interference of these long 
distance amplitudes with short distance contributions can explain the oscillations 
seen in the free space exclusive pp wide angle elastic scattering cross section. The long 
distance processes will be absent in the nuclear medium due to "nuclear filtering". 
With oscillations absent in the nuclear process, the theory predicts oscillations in the 
transparency ratio 180° out of phase with the free space oscillations [ 119|| . Future 
experiments should be able to clearly distinguish between these two explanations of 
the BNL data. 



Farrar et al. [41] have argued that at the hard scattering a wave packet of very 
small spatial dimensions is formed. The origin of this wave packet is ascribed to the 
hard scattering. Such a wavepacket is not an eigenstate of the strong interaction 
Hamiltonian and should expand during its propagation through the nucleus. The 
authors use two semiclassical models of this expansion namely i) the naive model 
in which the transverse separation increases linearly with time and ii) a model in 
which the expansion goes like the square root of the time. The results of Farrar et 
al are shown in Fig. (3.8) for the BNL experiment. In this paper and subsequent 
work no mention was made of the oscillations in the denominator of the transparency 
ratio. Nevertheless, data for the transparency ratio was fit to a model of the nuclear 
effects in the numerator of the ratio, thereby leading to a monotonically increasing 
transparency ratio. The agreement with experimental results shown in Fig. (3.3), in 
which the theory predicts a dramatic decrease beyond E beam = 12 GeV, is clearly very 
poor. The same is also true for the predictions for the SLAC experiment as shown in 
Fig. (3.6). Other consequences of the idea is presented in several reviews [^7, 48| [49 
of this approach. 
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Exploring the effects of the denominator in the transparency ratio, Hepplemann 
presented a plot of the numerator, that is the cross section for pA — > p'p (A — 1) 
(up to a normalization), a result shown in Fig. (3.9). Due to the experimental 
procedure this was obtained by multiplying s 10 T by the denominator in Eq. (3.1). 
The results are a reasonably flat s-dependence, supporting the idea that whatever 
causes the oscillations is filtered away in the nuclear target. 

Jennings and Miller J73L [7J], [75], [7J], [77j, [73] have given detailed treatments of 
modelling the expansion of the wavepacket, assuming it begins from a compressed 
initial state. They use non-relativistic quantum mechanical models and also attempt 
a less model dependent treatment in which their point-like wavepacket is expressed 
in terms of experimentally observed hadronic states. They calculate the transparency 
ratio using several assumptions such as the BdT|fj3] and the RP [119] mechanism. 
They find that their calculations do not agree with the results of BNL experiment 



with [[751 , 76 1 and without [74] the inclusion of BdT or RP mechanism. In reference 



77p, however, by inclusion of a more careful treatment of Fermi motion and of the 
longitudinal component of the momentum of the detected proton, the authors found 
agreement with the BNL data. 



Similar conclusion was also obtained by Anisovich et al who used the RP [ 119 



mechanism within a quark-diquark model of proton to calculate color transparency. 
They find that the BNL results are consistent with color transparency. 

Kopeliovich, Nikolaev, Zakharov and collaborators [EJ], BJ^, BJ)|, 107 | have also used 



their quantum mechanical model to calculate color transparency, with and without 
the inclusion of the interference of long distance components. Using theory technology 
adapted from the diffractive limit, Kopeliovich and Zakharov [|85| argue that their 
calculations are unable to reproduce the BNL results. In a recent review of their 
approach, Nikolaev [|105|| favors the RP explanation but argues that the experiment is 
kinematically flawed in coupling Q 2 and energy at fixed-angle kinematics. (From the 
hard scattering point of view, this flaw is a virtue.) Moreover, considerable latitude in 
varying these independently exists in the new BNL experiment due to more complete 



kinematic measurements. Nikolaev et al [|106|| also calculate the transparency ratio 



for the SLAC NE18 electroproduction experiment within their quantum mechanical 
model. They find that the onset of color transparency is very slow in their model in 
agreement with the NE18 results. The authors argue that this slow onset of color 
transparency is obtained because of the large number of excited states needed to 
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produce wave packets of small transverse size. 

There appears to be a general agreement among the theorists |23|, |119| , 0, [77], ^3 



that a proper understanding of the BNL results involves some interfering long wave- 
length contributions to the free space process, which are filtered in the nuclear 
medium. If this is the case then the hard scattering in nuclear medium cannot safely 
be assumed to be the same as the hard scattering in free space. It is reasonable to 
introduce a new method for analysing color transparency experiments in which the 
hard scattering rate can be treated as an adjustable parameter along with the atten- 
uation cross section ||68 |. This procedure has been applied to the BNL experimental 
data; details will be given in Section 3.4 below. The best fit requires that the hard 
scattering rate be significantly different from the corresponding rate in free space and 
also that the attenuation cross section is considerably smaller than 40 mb expected 



due to a Glauber model [68, 3S]. The nuclear medium hard scattering rate extracted 
by this analysis is in much better agreement with the rate expected from short dis- 



tance pQCD predictions of Brodsky and Lepage PH , compared to the corresponding 
rate in free space scattering. The attenuation cross section decreases as l/Q 2 again 
in agreement with QCD predictions. Based on this analysis of the BNL data, the 
transparency ratio for electroproduction in the SLAC experiment was also studied. 
In this case use of the pQCD short distance model for the hard scattering in the 
nuclear target is again in agreement with the experiment. 

There exist several studies which attempt to better understand the role of Fermi 
motion in color transparency experiments. It has been claimed by Jennings and 
Kopeliovich |7I| and by Bianconi, Boffi and Kharzeev |T3 that in electroproduction 
at the elastic point, the phenomenon of color transparency is possible only due to 



Fermi motion. However Nikolaev et al [|104|1 do not agree with this and show that the 
effects of Fermi motion are much weaker. We find no theoretical support for stating 
that color transparency requires Fermi motion, and will come back to this point in 
Section 5. 

Since the Fermi motion effects are controversial, it is clearly useful to reduce 
the dependence of the observed transparency signal on the Fermi motion. Such an 
attempt has been pursued by Frankel, Frati and Walet [(SJ who suggest an experiment 
in which color transparency effects are claimed to be separated from nuclear structure 
effects. The authors give a definition of the transparency ratio T which is argued to 
be independent of the spectral function. 
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3.3 Color Transparency in Diffractive Scattering 

We next turn to another very interesting set of color transparency experiments which 
involve diffractive scattering. Within this class there are both low Q 2 experiments 
and those with transverse momenta of final state hadrons larger than a GeV. These 
experiments are expected to yield important information about the interface of per- 
turbative and non-perturbative QCD. A basic process of interest involves the produc- 
tion of quark antiquark pair by a high energy virtual photon. The qq pair eventually 
converts into a vector meson which could be either p, p', ip, if)' etc. The kinematics 
of this process is shown in Fig. (3.10). The incident muon with four momentum k 
emits a virtual photon with momentum q which eventually forms a vector meson with 
momentum r with exchange of a momentum t with the nucleus. The momentum of 
the scattered muon is denoted by k'. The relevant kinematic variables for this process 
are the invariant mass squared of the photon — Q 2 — (k — k') 2 , the energy lost by 
the muon in the laboratory frame v = p ■ q/M, the four momentum transfer squared 
between the vector meson and the target proton t = (q — r) 2 and the fraction of the 
energy lost by the incident muon that is carried by the vector meson z = E p j v. Here 
E p is the energy of the vector meson and M is the proton mass. The experimentalists 
also define M t' — t — t min where \t min \ is the minimum value of \t\ allowed by the 
kinematics, which corresponds to the limit of zero angle between the momentum of 
the virtual photon q and the vector meson r. 

Suppose Q 2 > a few GeV 2 . The initial state is domnated by a qq pair produced 
with a transverse separation of about 1/Q, which is much smaller than a typical 
hadronic scale. The resulting meson is expected to exhibit color transparency in 



the nuclear medium [21, |25 |. An important consideration in such experiment is the 
production time for the qq pair and the time of formation of the vector meson. In light 
cone coordinates the production time t p is given by x + ~ 1/Q~~ ~ E p /(Q 2 + my), 
where my is the mass of the vector meson. We note that at the kinematic conditions 
of E665 experiment where E p was of the order of 100 GeV and Q 2 of the order of 1 
GeV 2 , t p is extremely large. For a wide range of kinematic variables it is vastly larger 
than the nuclear diameter. The time of formation Tp of the vector meson also turns 
out to be much larger in comparison to nuclear dimensions. It is roughly equal to 
the time that the qq pair takes to reach a transverse separation equal to radius of the 
hadron r#. This implies that Tp ~ th/vt, where vt is the transverse velocity of the 
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quark. The transverse velocity is of the order of kr/E p , where kx is the transverse 
momentum and E is the energy of the hadron. We then find that the formation 
time is of the order of rnEp/kx, which for large E p is considerably larger than even 
the largest nuclei. Due to these long time scales, the impulse approximation does 
not apply. However, the "frozen approximation" is reasonably invoked because time 
evolution is slowed by the Lorentz boost. Because of the large formation time the qq 
pair is expected to escape the nuclear medium well before it forms the hadron. This 
implies that the state that traverses the nucleus is best thought of as a well localized 
qq pair and not an ordinary hadron. 

FNAL E665: A Brief Review 

Depending on the value of t' the qq pair can be produced coherently or incoherently 
from a nuclear target. In fig. (3.11), we show the graph from the E665 experiment, 
Ref. [0], which shows the if distribution of events. The coherent and incoherent 
production regions are clearly separable in the curve for Calcium. For the events 
that correspond to coherent scattering the transparency ratio is expected to go like 
A 1 / 3 in the limit of infinite Q 2 when the nucleus is transparent to a point like color 
singlet object, since the scattering is now taking place along the entire length of the 
nucleus. Again, since the scattering is coherent over the entire nucleus the impulse 
approximation is not applicable and the process is not factorizable in terms of hard 
and soft scattering kernels. The interpretation of these experiments is considerably 
different from the hard scattering case. For the case of incoherent scattering the 
transparency ratio is expected to go to 1 in the high energy limit. 

The results of the E665 experiment for the case of p mesons are shown in Figs. 
(3.12-3.14). The typical values of Q 2 reported by the E665 group are of the order 
of 0.1 GeV 2 to about 10 GeV 2 . The conditions for color transparency are generally 
expected to hold only for Q 2 greater than about 1 GeV 2 but some authors argue 
otherwise. The beam energy is extremely large, 470 GeV. The momentum transfer 
to the nucleus, —t, ranges from about 0.1 GeV 2 to about 1 GeV 2 . The experiment 
observed a rise in the transparency ratio with Q 2 as shown in fig. (3.13), in the case of 
coherent scattering and in fig. (3.12) for the case of incoherent scattering, providing 
evidence for reduced attenuation with increasing Q 2 . 

The Q 2 dependence of diffract ive p production was predicted in advance of the 
experiment by Kopeliovich et al |37|. The authors calculate the transparency ratio 
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using the wave function \j* > of the qq fluctuation of virtual photons. This wave 



function was derived in Ref. [|103|| in the mixed (6, z) representation, where b is the 
transverse size of the hadron and z is a fraction of the photon's light-cone momentum 
carried by the quark, < z < 1. The authors argue that the important quantity that 
enters the calculation of transparency ratio is the product a(b)\j* > which peaks 



roughly at 26q, where bq ps 2/yQ 2 + m 2 is the size of the relevant qq fluctuation and 
a(b) is its interaction cross section. This implies that the process measures the wave 
function at 26q and therefore can be used to scan the wave function. In the calculation 
final state diffractive scattering leads to considerable off diagonal contributions. 

The authors give predictions for the transparency ratio for p, p', J/ip, ip', T 
and T' states. In contrast to their predictions for p', no significant difference in the 



transparency ratio was observed between p and p' in results presented so far [40 j. 

In diffractive scattering one probes an end point region with Q 2 of the order of 1 
GeV 2 and s >> Q 2 , at the interface of perturbative and non-perturbative regimes. 
For a wide range of kinematic variables, factorization of hard and soft scattering 
amplitude breaks down, and perturbative QCD is not applicable in a straightforward 
way. Experimental cuts were set to suppress events in which the photon produced 
other fast particles in the final state besides the p meson. Moreover the experiment 
did not exclude the possibility of peripheral pion production. Strictly speaking, 
there is no guarantee of the exclusive kinematics typically necessary to observe color 
transparency. A further complication is that in this end point region the process is 
not necessarily dominated by any short distance part of the hadron wave function. 
Nevertheless, the process showed a very interesting rise with Q 2 of the transparency 
ratio which can be interpreted as evidence for color transparency. This may be 
an indication that the extremely tight kinematic cuts sometimes felt necessary for 
observing color transparency are overly restrictive. 



3.4 Experimental Analysis Procedures 



Ratio method 

There exists currently two methods to test for color transparency. The ratio 
method, which is still widely used, was the original method introduced in the BNL 
experiment. The signal for color transparency was expected to be a monotonic rise 
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in T with energy such that it asymptotically approaches 1. This expectation is 
based on the assumption that in free space the hard scattering is dominated by the 
short distance components of the hadron wave function. If this is true, and if the 
hard scattering in the nuclear medium is the same as the free space scattering and 
further if the nucleus is indeed transparent to incoming and outgoing protons, the 
nuclear scattering should be identical to free space scattering. The data from the 
BNL experiment, shown in fig. (3.3), disagree with this expectation and showed a 
rise and fall of the transparency ratio with increasing energy rather than showing a 
monotonic increase. 

This anomalous behavior of the transparency ratio at first caused a great stir. 
By now it is clear that the structure is due to the process of making a division by 
the free-space pp — > pp data. On a log-log plot, that data falls roughly like s~ 10 , 
Fig. (3.15), but also are modulated by oscillations with the logarithm of the energy 
(Fig. (3.16a)). The oscillations are rather clear evidence for quantum mechanical 
interference of amplitudes. Apparently the numerator of the transparency ratio - 
which is the nuclear scattering - has little, if any, oscillations, so that the ratio must 
oscillate. There is good evidence for this: consulting Fig. (3.16b), the rise in the 
transparency ratio is exactly correlated with a falling oscillations of the free space 
data, indicating that the nuclear target is showing reduced oscillations. Making a 
plot of the nuclear differential cross section by itself [Q, with no division by the free 
space denominator but taking out a factor of s~ 10 , one sees that the oscillations, if 
present at all, are much reduced in the nuclear target Fig. (3.9). 

Evidence for Filtering 

We have strong reasons for believing the oscillations in the free space data are 
understood as part of a QCD phenomenon studied for more than 10 years. But 
temporarily setting aside the explanations proposed for the oscillations, if there is 
interference then the free-space pp — > pp amplitude has not settled into its asymptotic 
form. It must be different from the processes occurring inside the nucleus, which do 
not show the same oscillations. Model independently, something interfering in the 
free space pp — > pp process has apparently been killed in the nuclear target. There 
are two explanations proposed for this: 

Brodsky and de Teramond |23| propose that there are several charmed dibaryon 



resonances in the region of interest, which are eliminated in the nuclear case. The 
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states have not been observed elsewhere, but the value of the threshold kinematics 
do coincide with calculations based on charmed quark masses. By adjusting several 
parameters the shape of the differential cross section and the spin analyzing power A 
were reproduced. Presumably the charmed states could be observed to test the idea. 
It could also be applied to other reactions. 

Ralston and Pire ||119| proposed that the free space pp — > pp scattering consists 
of roughly two perturbative QCD regions, which have an energy dependent "chromo- 
coulomb" phase difference ||112| . One of the regions corresponds to "large" quark 
configurations which would be filtered away in the nuclear medium, leading to dis- 
appearance of the oscillations in the nucleus. By proposing that the oscillations are 
reduced for A>>1 the energy dependence of the BNL transparency ratio is repro- 
duced with no free parameters. The fact that such energy dependent phases occur in 



QCD has been confirmed by Sen [ 129 1 and Botts and Sterman [17, 16]. Spin observ- 
ables in this model of the free-space scattering have been fit by Ramsay and Sivers 



[126], and Carlson, Myrher and Chashkuhnashvilli [gfj. However, the fundamental 
parameters needed to reproduce the data have not yet been obtained from the (quite 
complicated) pQCD calculations. For example, the calculation of Ref. [|I7]] obtained 
an oscillation with too slow an energy dependence to fit the data. This is not too 
surprising, since the calculation was not set up to search for more rapidly varying 
phases, leaving open the possibility that the oscillations might still be calculable from 
first principles. 

The fact that there are "large" configurations in the hard scattering description 
of pQCD was not widely appreciated until recently. Such processes, which involve 
independent scattering of quarks at separated scattering centers, simply do not obey 
the assumptions of the "basic idea" presented earlier. Nevertheless, independent scat- 
tering occurs in the perturbation theory and must be dealt with to understand color 
transparency. Both of the explanations above agree on the proposal that the struc- 
ture in the transparency ratio occurs due to dividing a smoothly varying, somehow 
filtered numerator by an oscillating denominator. 
A Crossroad 

Conceptually the appearance of an oscillating transparency ratio brings us to a 
crossroad. One learned that the experimentally attractive transparency ratio was not 
a reliable observable to measure color transparency on its own. What else might be 
used? 
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We know that the experiment on the nuclear target measures a combination 
of the hard scattering rate and the nuclear attenuation. If one honestly admits 
that the nuclear hard scattering rates are unknown, how can one extract a signal 
of color transparency from the nuclear attenuation? Suppose, for example, that one 
overestimated the nuclear hard scattering rate- could this simply be compensated 
by a corresponding underestimate of the nuclear transparency? This ambiguity is 
not just a problem with the proton initiated processes but also, although less well 
appreciated, a difficulty with electroproduction experiments. 

Attenuation Method 



This method [gg, |69| is designed to test color transparency without making the 
assumption that the hard scattering in the nuclear medium is the same as the corre- 
sponding free space scattering. Instead, the nuclear hard scattering rate is extracted 
directly from data. 

The attenuation method is based on the observation that the slope of the curve of 
transparency ratio T vs. nuclear A contains information about the nuclear attenua- 
tion and can be used to extract an effective attenuation cross section. To understand 
how this is done, consider measuring experimentally an effective absorption cross 
section cr e // in some material. If the test beam flux is known, this is easy and can 
be done by ratios. But if the test beam flux is unknown - the case of undetermined 
hard scattering rates in the nuclear medium - then a ratio is not definitive. One 
can, however, study the absorption as thicker and thicker targets are used, and fit 
a e ff from the target thickness dependence. This is the basic idea of the attenuation 
method. 

Since the target thickness is set by the nuclear number A in color transparency, 
one should extract <r e // from the A dependence at fixed kinematic conditions (Q 2 and 
energy). This can be done rather model independently. In the fit there is an overall 
normalization, proportional to the hard scattering rate. There is a separate variation 
in the shape of the curve with A (e.g. the slope on a log-log plot). The shape of the 
A dependence is fit by varying cr e ff, while the hard scattering rate determines the 
overall normalization. These two independent observables are easily separated as a 
function of two "probe" variables Q 2 and A. 

To simplify the analysis, we make a few limited assumptions, which can be relaxed 
if needed. There must be an assumption that the observed process can be factored in 
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to hard and soft interactions. For A » 1 any variation of the hard scattering rate 
with A after filtering should be negligible compared to the rapid A dependence of 
the soft attenuation process. Finally, Fermi motion effects are modeled as necessary. 

Let us illustrate this in more detail. The impulse approximation can be used to 
separate mathematically the hard scattering and the process of propagation through 
the nuclear medium: 

measured rate(Q 2 , A) = [scattering rate (Q 2 )} [survival probability (cr e ff, A)] (3.3) 

Here by "measured rate" we mean data for either a cross section or a transparency 
ratio. If one studies the transparency ratio, then "scattering rate" is actually the 
ratio of the hard scattering rates in the nuclear target to the analogous free-space 
scattering. On the other hand, if (3.3) is applied to a cross section on a nucleus then 
"scattering rate" actually equals the hard scattering rate in the nuclear target. 

We next observe that the A dependence can be experimentally isolated by taking 
the logarithm of Eq. (3.3) at fixed Q 2 = Q^: 

log[measured rate(Q 2 , A)] = log[survival probability] + const(Q ) . (3.4) 

Here const(Ql) is a term independent of A: it comes from the hard scattering rate. 
Whatever the const, the attenuation cross section a cS is the parameter which deter- 
mines the shape of the survival probability as A is varied. Thus to determine a c s we 
want a fit to the shape of the A dependence of the right-hand side, but a fit only up to 
an additive constant. In practise this is quite easy. For example, plot the measured 
rate as a function of A on log-scale paper: adding an overall constant is just a rigid 
translation of the fit up or down on the paper. 

Relating the shape of the survival probability to cr cS requires a theoretical model. 
If one uses a simple exponential attenuation model this dependence is quite explicit 
and straightforward. If this fits the data set well then one effective cross section 
<7 e ff at each Ql seems to be all that is observable with the data set. Other models 
are possible: for a huge class of models of the survival probability the attenuation 
procedure applies just the same. 

After determining the survival probability, the const (Q 2 ) values contain the 
information about the hard scattering rate. The interpretation depends on whether 
(3.3) is applied to the transparency ratio or a cross section on the nuclear target. If 
"measured data" is a cross section then one finds the scattering rate in the nuclear 
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target directly. No normalization needs to be set using isolated hadron events. If 
there is nuclear filtering we expect the empirically determined hard scattering rate 
generally will disagree with the scattering rate of free space data. 

This method was applied to the BNL data. The results provide evidence for 
observation of color transparency and nuclear filtering. pB"| , |6T?[ The A dependence of 
the data at different fixed Q 2 is shown in Fig. (3.17). The data was analysed by 
extracting the cross section in the nuclear target by multiplying the transparency 
ratio by the known free space cross section. This yields the nuclear cross section 
divided by the number of protons in the nucleus Z. By assuming some factorization, 
this quantity can be written as 

s w ^^ = N(Q 2 )P ppp (a cS (Q 2 ), A) . (3.5) 

The factor s 10 and 1/Z are inserted to take out the typical order of magnitude of the 
cross section; their usage introduces no bias. The factor P ppp , which is the survival 
probability, contains information about the attenuation of the three protons that 
cross the nucleus in pA — > p'p" (A — 1). The factor N(Q 2 ) contains information about 
the hard scattering. In equation (3.3) we have made the reasonable assumption that 
once A >> 1 then N(Q 2 ) is independent of A. In other words, the hard scattering 
rate saturates once A is large enough. We return to this point below. 

To proceed further one must model the survival probability. The simplest model 
assumes that protons are exponentially attenuated in their passage through the nu- 
cleus. Thus the survival probability of a proton after propagating a distance z through 
nucleus is equal to exp(— pa e ffz). The survival probability for the three protons can 
then be easily modelled by a Monte Carlo event generator by treating the a c s as a 
free parameter. Using the Monte Carlo calculation the data were fit to determine a e s 
and N(Q 2 ) for each value of Q 2 at which the data is available. The resulting fits for 
the BNL data were found to be: 

E a e ff(E) N(E) x 2 

6 GeV 17±2 mb (5.4 ± 0.4)C 0.28 
10 GeV 12±2 mb (3.3 ± 0.4)C 0.53 

where \ 2 = ^[(Vi ~ di) / Adi] 2 / N , di are the N data points, Adi is the error in di 
and yi are the theoretical values calculated using the Monte-Carlo, and Q = 5.2 x 10 7 
mb GeV 18 is a constant containing the overall normalization of the cross section. 
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The Q 2 values for the two beam energies 6 and 10 GeV are 4.8 GeV 2 and 8.5 GeV 2 
respectively. The fit to the data (Fig. (3.18a,b)) clearly shows that <7 e // decreases 
with energy and is considerably lower than a typical inelastic cross section of 36 mb. 

Let us return to the issue of saturation of the hard scattering rate for A » 1. It 
should be clear that the method is very general, and can easily be modified to model 
the A-dependence of hard scattering due to filtering. We will not do that here, but 
instead look to the data to estimate small- A effects. As remarked in Ref. [BBl, the 



effects on the fit of discarding the Li data point are negligible. We regard the issue 
of whether the data point for C might be discarded as an empirical question. It 
should be noted that even with the choice to discard A < 26, the last 3 data points 
in Figs. (3.18a,b) are sufficient to eliminate a 30-36 mb absorption cross section 
as a contender. The Attenuation method is thus considerably more definitive than 
the Ratio method, and a standard "Glauber Model" has been ruled out by the A 
dependence at each fixed Q 2 . 

Having fit the data's A dependence, the fit can be examined at fixed A to find 
the survival factor as a function of Q 2 . The BNL data was taken only at two energies 
for a good range of A, so N(Q 2 ) and a e ff(Q 2 ) were found at two points. However, 
for the Aluminum target the data at intermediate energies was also given. One can 
make a function for N(Q 2 ) which interpolates with Q 2 smoothly between the two 
endpoints where Q 2 was reported. It is given by: 

»(Q>) = (QV J^ 2) o,a • (3-6) 

Then inverting (3.5) one can determine the survival factor P ppp as well as cr e // 
as a function of Q 2 by using the data for Aluminum at intermediate energies. The 
resulting best consistent fit for P ppp is shown in Fig. (3.19). The corresponding 
a e ff(Q 2 ) is shown in Fig. (3.20). The survival factor would be flat with Q 2 if one 
used a traditional Glauber model. The fact that the survival probability rises with 
Q 2 is clear evidence for color transparency. 

A final consistency check involves looking at the the normalization factor. After 
taking out the nuclear attenuation effects, according to the perturbative treatment 
the Q 2 dependence of N{Q 2 ) is due to the hard scattering process. It was found 
that N(Q 2 ) decreases relative to s~ 10 , meaning that the hard scattering rate in the 
nuclear target is decreasing faster than the naive quark counting model prediction one 
would use at short distance after nuclear filtering. In perturbative QCD, however, 
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the quark-counting prediction is modified, due to the running coupling a s (Q 2 / A 2 cd ) 
and scaling of distribution amplitudes. The short distance QCD prediction goes like 
a] because there are five gluons in the amplitude: 

da/dtpQCD = (a s (Q 2 /A 2 qcd )) 10 s- w f(t/s) (3.7) 

It is interesting to compare the form including powers of a s with the hard scattering 
rate in the nuclear target. First, recall that the Q 2 dependence of a] causes serious 
disagreement with the data for isolated pp scattering in free space with A qcd « 100 
MeV for the quark-counting model. We would like to see whether the nuclear tar- 
get has filtered the events down to something like the shortest distance component. 
Although (3.7) may appear naive, it is adequate because of the usual ambiguity in 
the choice of scale. One can generate a range of reasonable theoretical predictions 
by choosing the Q 2 scale of a s in a typical range (—t/2) < Q 2 < (—1.5 t). To 
improve on this theoretically requires a calculation of next to leading logarithms. 

For comparison with the nuclear target, the ratio of the global s~~ 10 N(Q 2 ) to the 
pQCD predictions Eq. (3.7) was considered. The results are plotted as solid lines 
in Fig. (3.21). The asymptotic prediction |BT| is also illustrated as a dotted line. It 



falls within the region of scale ambiguity of the pQCD predictions. For comparison 
we also show s~ 9 ' 7 , a dashed line which is well outside the range of the perturbative 
predictions. The solid lines of s~ 10 N(Q 2 )/(da/dt pQ cD(Q 2 )) are rather flat, showing 
good agreement of the BNL nuclear hard scattering data with short distance QCD. 
For the Aluminum target we have several data points at several Q 2 which also fall 
fairly well within the band of perturbative predictions. These are rather spectacular 
results, but we emphasize that they should be viewed with caution because the BNL 
experiment was a pioneering one. If they are confirmed by upcoming higher precision 
data, it will lead a great deal of strength to the idea that QCD is cleaner after filtering 
in a nuclear target. 

A Scaling Law 

We next discuss a scaling law ||114|| predicted for the survival probability. Be- 



cause the scaling law predates the attenuation method, it was first applied to the 
transparency ratio, under the assumption that hard scattering rates might cancel 
out in electoproduction. The scaling law says that the survival probability is a func- 
tion of a dimensionless variable. The important dimensionless variable that it can 
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depend on is the effective number of nucleons encountered by the protons as they 
propagate through the nucleus (Fig. (3. 22)). This is proportional to the length of the 
nucleus (A l ^Fm) times the nuclear density n times the effective nucleon-nucleon 
cross section cr e ff- If the cross section goes like l/Q 2 , then the survival probability 
is a function of the dimensionless quantity nA 1 / z Fm/Q 2 . 

The advantages of a scaling law in testing dynamical assumptions are very great. 
Rather than testing values of ratios which can be fit with many parameters in some 
model, the scaling law method allows an almost model-independent check of the most 
basic assumptions. 

To experimentally determine whether this law is satisfied or not, consider the 
survival probability to be a function of Q 2 /A a and determine a from the experimental 
data. A basic complication is the fact that the experimentally measured transparency 
ratio T involves an unknown function f(Q 2 ), where f(Q 2 ) is the ratio of the hard 
scattering in the nuclear medium to the free space scattering. With a good data set 
over a range of the Q 2 and A plane one can nevertheless extract functional forms for 
both f(Q 2 ) and the survival probability P by fitting T = f(Q 2 )P(Q 2 /A a ). This is 
attractive because of its model independence: no theory needs to be used to model 
P, which is simply determined by the best fit. 

The data from the BNL experiment is limited, available only at two energies 
for several A. We therefore take a somewhat different approach here to check the 
scaling law with this data. We assume that the hard scattering inside the nuclear 
medium satisfies the short distance pQCD predictions. This is based on the idea 
of nuclear filtering; is supported by the global fit, and by the fact that the cross 
section in the nuclear target does not show any oscillations. We therefore set the 
average hard scattering cross section in Eq. (3.5) to be Z times the short distance 
perturbative QCD prediction given by Eq. (3.7). (In Eq. (3.7) we ignored the 
anomalous dimension which is found to make a very small difference in the results.) 
We can then calculate f(Q 2 ) and take it out of the transparency ratio to yield the 
survival probability. The survival probabilities for three different values of a are 
plotted in Fig. (3.22). The data points are circles (Q 2 = 4.8 GeV 2 ), squares (Q 2 = 
8.5 GeV 2 ) and crosses (Q 2 = 10.4 GeV 2 ). It is clear from the figure that a — 1/3 
is favored over other values of a. 

(There remains one subtlety. Although the scaling law procedure is quite explicit 
about the functional dependence on the two variables Q 2 and Q 2 /A a , the ansatz 
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(1) has a symmetry: the data and fit are unchanged under f(Q 2 ) — > f(Q 2 )/K, 
P(Q 2 /A a ) — > KP(Q 2 /A a ), where K is any constant. The absolute magnitudes of 
the survival probability (or the hard scattering) and the effective attenuation cross 
section are therefore not determined by this method.) 

Finally we note that the method of looking for scaling relations does not make any 
assumptions about the free space scattering. It only assumes a factorization in the 
nuclear medium. Therefore it is applicable irrespective of the source of long distance 
contamination in the free space scattering. 

The method has been criticized by several authors on the grounds that due to 
the time evolution of the hadron wavefunction as it propagate through the nuclear 
medium, the attenuation cross section changes with time. These authors claim it 
is not safe to treat it as a constant as was done in the original application of this 
method. We must point out that this is not a criticism of the method but of the 
model used to calculate the nuclear attenuation. Other models are possible, and 
might be used to predict other scaling laws to test these assertions. For a crude 
example, if a e ff increases linearly with time, then one might get an additional A 1 ^ 3 
dependence. Such models would be tested by scaling in Q 2 /A 2 ^ 3 . One can always 
extract an average attenuation cross section a e g; if there is color transparency in the 
data, this information will be directly available by a plot of a e g vs. Q 2 . It would seem 
logical to explore such broad simple regularities in theory and data before getting 
mired in the details of particular models. 

To summarize, the BNL data yielded an attenuation cross section which was 
consistent with the theoretical expectation that it should decrease like 1/Q 2 - More 
interestingly, the data yielded a hard scattering rate which is consistent with the 
theoretical expectations from short distance pQCD. 

Experimental Analysis of Electroproduction 

We next discuss the electroproduction experiments on color transparency. The 
NE18 group at SLAC has already published results|)4], |109|| of measurements of eA — > 
e'p(A — 1) on several nuclear targets with maximum momentum transfer Q 2 being 
6.8 GeV 2 . There also exists a proposal at CEBAF to perform a similar experiment 
with better statistics. These kind of experiments are beautifully complementary to 
the hadron beam experiments. The initial state of the large Q 2 virtual photon is very 
clean and well understood: one is measuring the photon's absorption and final state 
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interactions in the nuclear target. There is tremendous interest in this and other 
experiments which could measure meson knockout, transition form factors, etc. A 
clean theoretical matrix element is measured in the experiment: 

<p,A-l\J£JA> 

However, the experiment is somewhat insensitive to the attenuation cross section of 
the system crossing the nucleus, so that comparatively high precision data may be 
needed to extract an attenuation cross section. This is shown in Fig. (3.23), which 
shows some predicted survival probabilities for several attenuation cross sections a e g. 
Even with cross sections varying by a factor of about three, the curves on the log-plot 
are almost parallel. It follows that even for a substantially decreased a e g, the effects of 
color transparency (if any) might be masked by a simple change in the normalization. 
We have already seen that such a change in normalization is equivalent to a change 
in hard scattering rate, so it is tricky in electroproduction to separate the ambiguities 
of hard scattering inside the nucleus from the extraction of c e //- 

Qualitative Features of A-Dependence in Electroproduction 

To understand this in more detail one can make an analytic calculation. Assuming 
a uniform nucleon number density n — 1/6 Fm 3 and nuclear radius Ra = 1.2 FmA 1 / 3 , 
one can find the survival probability P p of one proton taken at random over the 
volume of a nucleus and crossing on a straight line with exponential attenuation. It 
is a function of y — Ra/1, where / = l/(ncr e ff) is the mean free path: 

P> cff , A) = P p (A 1/3 a cS ) = A + A e -2, _ JL (1 _ e -2y } ^ 

This formula was used to generate Fig. (3.24). Even though there is an integration 
over a sphere, this formula depends on only one variable y ~ A l l z na e ^. In the 
limit I » Ra we have 100% survival. In the limit of / << Ra the problem is also 
quite simple. Particles survive if they are knocked through no more than thin skin 
of depth I on the "back" side of the struck nucleus. The volume containing such 
survivors is nR A l; the area 7iR\ occurs because we want the projected area along 
the beam direction. The survival probability, in this limit, should be the ratio of the 
volume containing survivors to the total volume: 



P P {1 « R A ) 



nR A l 3 3 



4ttR a /3 Ana^A 1 / 3 4y 
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Now the ambiguity in examining data in this case is brought out by considering a 
transparency ratio 



where Na(Q 2 ) and Ni(Q 2 ) are nuclear and free-space hard scattering rates, and we 
indicate the possible Q 2 dependence of a e s(Q 2 ) explicitly. The left-hand side of the 
equation above is measured. But in the limit chosen there is no model-independent 
way to extract a c s{Q 2 )\ That is because there are two unknowns, and we can let 



where \(Q 2 ) is any function, and still fit the data. 

The way out is to use the A dependence to vary the target thickness over a 
reasonable range of I/Ra- With such a procedure the BNL data was fit using the 
shape of the A-dependence, resolving the ambiquity. The same procedure can be 
done in electroproduction, but it is more difficult. One can compare Fig. (3.23) to 
the survival P ppp of three particles in a nuclear target, all having the same exponential 
attenuation law, averaging over the volume, and with the 90° kinematics of the BNL 
experiment. This calculation was also done with realistic nuclear densities, which 
show considerable edge effects. To a surprisingly good accuracy, the survival of three 
particles is proportional to the cube of the survival of one, for A > 7. (The empirical 
proportionality constant is given by P ppp 1.9P p .) Because of this cube effect, the 
curvature in the A-dependence in the (p, 2p) experiment is easy to pick out. 

However, in electroproduction a high probability of survival occurs because the 
attenuated "particle" only has to cross a small part of the nucleus once. This means 
that the experiment does not necessarily have a high discriminating power in measur- 
ing the attenuation cross section, unless the data is of comparatively high precision. 
The reason for emphasizing this point is simply to balance the virtues and com- 
plementarity of each experimental approach. The virtues of electroproduction, in 
precision and understanding of the probe, are so well known that it might be helpful 
to mention the ironic limitations. Perhaps a CEBAF experiment could resolve the 
question. 

Q 2 Dependence 

The results of the SLAC NE-18 data as a function with Q 2 show only a mild 
increase with Q 2 compared to a plane wave impulse approximation. Based on the 




N A (Q 2 ) - KQ 2 )N A {Q 2 ) ■ a cS {Q 2 ) - A(Q> eff (Q 2 ) 
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scaling of the spectral densities of the NE18 |38[ experiment with Q , a contant 
transparency ratio was anticipated in Ref. |57|] well before the experimental results 
on the ratio were announced. This paper also preceeded several papers based on 
preliminary SLAC results. Plots [TUD| of the A dependence at fixed Q 2 do seem to rise 
systematically with increasing Q 2 , but it is not a very big effect. If the transparency 
ratio used is interpreted to measure the survival probability, then the experiment 
shows only a weak signal for color transparency. However, we have already seen that 
the transparency ratio cannot generally be interpreted as a survival probability. All 
that the NE18 data can determine now is some curve in the o~ e ff(Q 2 ) versus Na(Q 2 ) 
plane. If the electroproduction data becomes good enough to extract an attenuation 
cross section along with a hard scattering rate, then we suggest that would be the 
preferred procedure. It is possible that CEBAF (or future facilities such as ELFE) 
experiments will have sufficiently high statistics to make this possible. 

Interpretation of the SLAC Nuclear Hard Scattering Rate 

Theoretically, the description of the hard quasi-elastic scattering in electroproduc- 
tion begins with exclusive form factors: the elastic form factor, and transition form 
factors to whatever states end up converting to protons after travelling through the 
nucleus. There are plenty of reasons to believe that the electromagnetic form factors 
used to describe the process in a nuclear targets may not be the same ones as in free 
space. Moreover, the free space processes at laboratory values of Q 2 are not truly 
short distance dominated. Isgur and Lewellyn- Smith (ILLS) p5| and Radyushkin 



[ |118|| have made the case for this, claiming that none of the perturbative arguments 
apply at all. But if this is so, how can one explain the many regularities of power 
law behavior which are observed in many experiments? It would be unscientific to 
ignore this qualitative success! More recently, Li and Sterman and Li have shown 
that the ILLS arguments focus on an integration region that does not occur. In the 
new estimates, less than about one-half of the amplitudes in free space scattering 
could come from dangerously large quark separations.! The new ingredient is the 
appearance of Sudakov effects in the hard scattering form factor. Would these offen- 
sive regions be filtered away in a nuclear target? Nobody yet knows, but there are 
several indications that it should happen. We would then see a repeat of the BNL 
phenomenon: for experiments inside the nuclear target there would be a more pure, 

*This may sound like a negative conclusion but it an improvement on ILLS who claim that not 
even a tad of the amplitude is ever short distance. 
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short distance form factor, and a s should make its appearance in the naive way for 
the hard scattering. 

In Ref. ||68|| , the transparency ratio for electroproduction has been calculated by 
making the assumption that the elastic form factor of nucleon scales like a^/Q 4 . The 
free space form factor scales like 1/Q 4 . Taking the attenuation cross section, Fig. 
(3.21), extracted from the BNL data and calculating the survival factors P p using 
the formula given above, the authors obtained the curve for transparency ratio which 
are reproduced in Fig. (3.24) for the Fe target. The theoretical calculation is flat in 
approximate agreement with data. The transparency ratio calculated in Ref. [68 is 



defined slightly differently from the one reported by the NE18 group, but this should 
not effect the predictions significantly. In Ref. |68|] only the momentum dependence 
was predicted since the over all normalization of the elastic form factor in the nuclear 
medium is not known. 

A better procedure to interpret the NE18 data may be to use the Li and Sterman 
[ |92]j analysis for elastic form factor in free space and then model the filtering of the 
long distance part of the amplitude in the nuclear medium. It would be interesting to 
extract the effective attenuation cross sections for the short wavelength pieces which 
should be the only amplitudes that survive in a large nucleus. 



3.5 Ideas for future experiments 

In this section we review some ideas that have emerged in the literature for future 
experiments related to color transparency and nuclear filtering. There exists an 
enormous amount of work in this direction by many authors. We will also discuss 
the near-term experimental viability of some of these proposals. 

1) Exclusive Channels 

There exist a large number of different exclusive channels where we expect color 
transparency to occur. Examples include eA — > e'p(A — 1), eA — > eirA, eA — > 
eKA(A — 1) in electron accelerators, ttA — > p(A — 1)tt', pA —>■ p(A — l)p', pA — > 
7t + tt~(A — 1), pA — > K + K~(A — 1) in hadronic machines. 

Furthermore the corresponding hadron initiated processes where the nucleus is 
replaced by a proton are expected to show oscillations which will be absent in the 
nuclear case due to filtering of the soft components. This will lead to oscillating 
transparency in all these cases. A higher energy pA — > p'p"(A — 1) experiment will 
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also greatly help in distinguishing between the two proposed explanations [£3|, |119 
for these oscillations seen in the current BNL data. If the independent scattering 
interference mechanism is correct then the transparency ratio should continue to 
show oscillations at higher energy, 180° out of phase with the oscillations in the free 
space data. The charmed dibaryons explanation will require no further oscillations. 

As discussed in detail in Sections 4 and 5, the nuclear collision will also yield the 
short distance wavefunction of the participating hadrons, thereby giving important 
dynamical information about the hadron. It should be possible to use the wavefunc- 
tions measured in these experiments to predict the results of other experiments. 

A high statistics electroproduction color transparency will be performed at CE- 



BAF [39] which should help in greatly clarifying the relevance and precision of SLAC 
data. The higher energy machine ELFE might, however, be necessary to see a clear 
signal of color transparency in electron initiated experiments. 

2) Helicity Conservation 

Another set of very interesting nuclear medium effects arise because of the "hadron 
helicity conservation" rule in hard exclusive collisions. In exclusive large momentum 
transfer hadron hadron collision A+B — > C+D, a general prediction of short distance 
dominance is that 

X A + X B = A c + Xd (3.9) 

where A/ is the helicity of hadron I. However, this rule has been observed to be 
violated in several processes involving free space collisions. That hadron helicity non- 
conservation is so widespread proves conclusively that many processes must contain 
considerable contributions from long distances . A generic test of filtering to short 
distances in the nuclear medium is restoration of hadron helicity conservation. 

Many simple and direct tests of this idea can be done with polarized beams. 
These may be become available at BNL or future facilities, perhaps including the 
LISS facility now under discussion. It is also quite useful to use mesons and 
baryons which are "self-analyzing" in the final state. The A baryon and p meson are 
standard examples. Consider, for instance, the reaction 

rrA — > pp(A - 1) 

which can be measured without polarized target or beam. Hadron helicity violation 
has been measured here in the free space case [E3] at energies O(20GeV) and large 
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angles. A direct check of nuclear filtering is the restoration of hadron helicity conser- 
vation in large nuclei. By our counting this effect is due to an extra factor of b 2 , and 
will scale like A 1 ^ 3 /Q 4 . The helicity density matrix of the produced vector meson is 
analysed from the angular spectrum of its tot decay. 

The situation for electron beams is quite unsettled experimentally. A recent SLAC 
experiment had made the separation of the electric and magnetic proton form factors 
at high momentum transfer. The data for the proton's helicity flip form factor F 2 
is consistent with a power law decreasing like 1/Q 6 - Theoretically, the electric form 
factor F 2 (Eq. 2.1) is a higher twist object difficult to analyse within QCD. However, 
it has been pointed out in Ref. |)6| that the spin structure of this observable makes it 
possible to analyze it within pQCD. It is claimed that F 2 is determined by non-zero 
quark orbital angular momentum, which involved wave functions with a node at the 
origin. Because these wave functions are zero in the short distance limit, nuclear 
filtering should preferentially deplete them, providing a novel way to measure the 
wave function. Many interesting experiments can be done at moderate, fixed Q 2 . 
Thus one predicts that the analogue of F2/F1, when extracted from measurements 
in a nuclear target, should show a strong decreasing function of A. One also expects 
that the magnitude of the analogue of F 2 when measured in a large nucleus should 
be reduced compared to the free space value. This idea may be tested in the future 
as part of a proposal to CEBAF to measure form factors of nucleons in nuclei by 
polarization transfer 

3) Recoil Polarization 

It is well known that proton-proton elastic scattering at a few GeV of energy 
and small momentum transfers produces a spin polarization normal to the scattering 
plane. Since this violates the hadron helicity conservation rule, one can ask whether 
the effect might be reduced by nuclear filtering. An approved CEBAF proposal [|128 



to measure final state polarization from a hard-struck proton traveling through the 
nucleus will study eA — > e'p(A-l). A scattering plane is produced by comparing the 
momentum of the virtual photon and the final state proton. To measure the recoil 
polarization P n normal to this plane, one needs a polarimeter at the focal plane of the 
hadron spectrometer. The theoretical analysis of this experiment is tricky, because 



it is completely quantum mechanical; Greenberg and Miller |58| have a detailed 
treatment within the context of a compressed-state model. It is not clear whether 
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the energy at the first CEBAF machine will be high enough for color transparency 
to be observed in this way. A study at much higher momentum transfer will require 
a higher energy machine such as the European ELFE project @, |[. In that case, 
making a precise measurement of P n will also require substantial progress to be 
achieved in proton polarimetry at high energies. One might also look for correlations 
of polarizations of resonances such as the delta with the scattering plane to exploit 
the self-analyzing final state. 

4) Excited States 

An important element of color transparency experiments is that only the short 
distance part of the hadron wave function survives upon crossing the nuclear medium. 
If the short distance wave function for a hadron is close to zero, then the production 
of such a hadron will be considerably suppressed in a nuclear target. This leads to the 
very interesting prediction that the radial and orbital excited states of hadrons will be 
considerably suppressed when produced in a nuclear medium. This phenomenon can 
be observed in electroproduction experiments where the photon propagates e.g. as a 
p or p' meson through the nucleus. From time-scale considerations, one guesses that 
the formation time of the hadron has to be smaller than the longitudinal dimension 
of the nucleus to find this suppression. Otherwise the hadron is formed outside the 
nucleus and no significant difference should arise between different excitations. We 
believe that at ELFE and/or CEBAF energies this would be possible. Fermilab 
energies, however, turn out to be too large for this effect. A better picture there 
is that a quark pair rather than a physical meson propagates through the nucleus; 
the meson is typically formed long after the interaction of the quark pair with the 
nucleus. 

5) Scanning the Hadronic Wave Function 

There exist a large number of diffractive resonance production experiments which 
are expected to show color transparency. This effect may have already been confirmed 
by the experiment E665 on diffractive p production. Kopeliovich, Nemchik, Nikolaev 
and Zakharov |57 have discussed an interesting set of future experiments involving 
photoproduction of other resonances and their excited states. The authors argue 
that it is possible to scan the wavefunction from large to small distances by changing 
the momentum of the photon. Some of this is the same as nuclear filtering and the 
behavior of large Q 2 in selecting the short distance parts of the wave function, in 
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cases where the impulse approximation can be used (Section 4.5). Other aspects 
are unique features of the authors' theory of diffractive physics. Emphasising off- 
diagonal contributions and quantum mechanical mixing, the radially excited states 
are predicted not be suppressed to the same extent one might expect if only the short 
distance contributions were dominant. 

Quasi-exclusive meson production experiments can also be used to gain informa- 
tion about the distribution functions of the meson, as mentioned by Carlson and 
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Nikolaev has made numerous suggestions for CEBAF experiments, which can be 



found in extensive reviews ||105|| . A diffractive picture is assumed and often used quite 
imaginatively. A novel recent paper |J predicts color transparency for ss mesons 
produced under conditions of extremely low momentum transfer. This experiment 
seems to be ideally suited for CEBAF and it would be interesting to know the results. 

6) Diffractive Jets, and Explosive Destruction of the Proton 

It should be interesting to study diffractive scattering of pions on nuclear targets 
at very high energies in the process ttA — > jets + X. Frankfurt, Miller and Strikman 



4S] have suggested that color transparency might be observed in kinematic situations 
where two opposite isolated jets carry most of the beam energy. The idea is that the 
two-quark Fock component of the pion might contribute significantly, and that short- 
distance should occur due to the large momentum separating two observable jets at 
Fermilab energies. A shortcoming of this proposal is that it is quite vague, not 
specifying a definite cross section nor kinematic conditions to make the experiment 
very definite. For example, jet calculations in QCD are usually quite sensitive to 
the definition of the "jet". It follows that the normalization estimated in Ref. [f|8 



for the process, and claims that the x dependence of the distribition amplitude can 
be measured, cannot be considered reliable. The authors have suggested eliminating 
some of these problems by finding a signal which goes like A 2 , indicating coherent 
scattering with the nucleus, controlling the color flow this way. 

We suggest that the experiment should include a veto in the forward direction 
to eliminate co-moving higher Fock components. Such a veto, plus resolution of the 
total jet energies, might be sufficient to extract the 2-quark, short distance compo- 
nent. There seems to be no compelling reason to insist on the super- low momentum 
transfers of the A 2 dependence, because color is separated in the event: there is only 



68 



an implicit assumption that the formation of jets will not disturb the color trans- 
parency of the separating pieces. The region of diffractive scattering off separate 
nucleons therefore need not be excluded. Finally, there is no a-priori reason to use a 
pion beam. Using a proton beam, one could also look for 3 jets carrying the beam 
energy (within resolution), with a veto on other particles in the forward directions. 
In this way one might study "explosive destruction of the proton" into 3 independent 
valence quarks, and correlate this with ideas about the 3-quark Fock representation. 

7) Intrinsic Charm 

The idea of an intrinsic charm component of the proton wave function has emerged 
from the phenomenological analysis of charmonium production at large xp. It has 



been suggested [22] that the higher twist Fock component | uudcc) was responsible for 
leading charm production. When the proton scatters on the nucleus, the coherence of 
the Fock components may be broken and charmonium systems may be formed from 
the fluctuations. Now, the nuclear dependence of this process is expected to be typi- 
cal of soft interactions, i.e. a a ~ A 2 / 3 un and not linear in A as in incoherently adding 
physics. Thus, the nucleus again should act as a filter against intrinsic charm compo- 
nent contribution to charmonium production, in favor of a perturbatively generated 
cc pair from gluon emission. 

However, charmonium production and absorption in nuclear medium is a subject 
so full of controversies that it might be difficult to test color transparency ideas in 
these processes. A major barrier is the uncontrolled color flows, and the multitude 
of different processes, which contribute to a inclusive processes producing charm. 
Closely related is the observed suppression of J/ip and bottomonium production on 
nuclear targets. In that case, it is not settled yet whether the suppression is due to 
big absorption of the cc pair in the nucleus, or reduced production due to energy 
losses in the intitial state. 

8) Anomalous Color Transparency 

Consider the reaction 7 A — > 7r°A in the region of small momentum transfer —t < 
123j| . Experimentally a region of t can be isolated showing that the production 



m 



of the 7i° occurs by the "Primakoff effect", namely the 7771" vertex made by the 
incoming photon, the nuclear Coulomb field, and the 7r°. Most production by this 
mechanism occurs just before the tt° strikes the nucleus. 

This process may show a signal of "anomalous" color transparency ||123|| when 
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studied as a function of A and the 7r° energy. Although all momentum transfers are 
small, the 7771° vertex is known to be controlled by short distance: it is the original 
perturbative "chiral anomaly". If the qq pair contributing to tt° production are 
close together, then the nuclear attenuation measured just after the pair is formed 
should show small attenuation. Signals of this would occur using either the ratio 
or attenuation method; an example of predicted energy and A dependence of the 
transparency ration is shown in Fig. (3.25). The kinematics of the 7r° must be 
tuned to find a region where Primakoff is operating while a good fraction of the 
7r°'s produced penetrate the nucleus, but a window exists where this is possible. 
Because the anomaly is an object of great intrinsic interest, a "snapshot" of its short 
distance characteristics would be extremely interesting and a rather direct probe of 
fundamental microstructure of chiral symmetry. 

9) Virtual Photon Polarization Dependence 

It has been argued that transversely and longitudinally polarized photons may 
interact with the proton in different ways H, based on the possible occurence of 
protons in "compressed" fluctuations. For instance, if only transverse photons interact 
with compressed protons, the color transparency is claimed to depend on the degree 
of virtual photon polarization e through the equation 

R = R G + [R{Q 2 ) - T/i ' 



t/j% + e 



where r = -Q-^ , M p and [l v are the proton mass and magnetic moment, Rq the 
Glauber value for transparency and R(Q 2 ) a naively calculated transparency ratio. 
A different formula is found if only longitudinal photons interact with protons. The 
net effect of such a picture is an electron energy dependence of color transparency at 
fixed Q 2 . 

10) Double-scattering Events 

Studying the Q 2 dependence of the number of final state interactions of recoil 
protons in quasi-elastic scattering of electrons on light nuclei is a type of inverse way 



to look for color transparency. It is claimed that rescattering events such as 

e + A -> e + p' f + p T + (A - 2) 
have a cross-section which decreases much for increasing Q 2 in light nuclei. Model 
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calculations for the 4 if e case show that measurable effects can be obtained provided 
a full acceptance detector is used. 

11) Missing Transverse Momenta 

A way to access final state interactions (and its decrease through color trans- 
parency effects) is to study the broadening of the momentum distribution transverse 
to Q of scattered protons in (e,e',p) reactions. In the case of light nuclei, the de- 
scription of this broadening is dominated by single rescattering of the struck proton. 
An accurate theoretical description is however needed [ 104 | for rescattering without 



color transparency to separate CT effects from other Q 2 dependences. 
12) Transparency in semi-inclusive processes 



Rinat and Jennings [127] have extended the concept of transparency to also in- 
clude semi-inclusive processes. They argue that the transparency measured in these 
processes is not the same as the one measured in purely elastic scattering. Although 
this work is a commendable first attempt, we mention it also as an example of an 
area where much more conceptual work could be done. 

13) Color Opacity in the diffractive region of pp scattering 

The pp elastic cross section at very high energy (E rs 1000 GeV) with t of the 
order of 10 GeV 2 , Ref. |7D|, |oT|, shows a remarkable agreement |J3J with the t~ 8 power 
law behavior predicted by the independent scattering mechanism^]. In Fig. 3.26 
we show this for the reader's inspection. Based on this agreement, Donnachie and 
Landshoff have argued that the independent scattering amplitude dominates the 
scattering cross section in this region. This idea can be further tested by studing this 
kinematic region of pp scattering using nuclear targets. We predict strong attenuation 
of the protons participating in this process, and t~ 10 behavior if nuclear filtering 
succeeds in selecting the shortest distance components for A >> 1. 

In conclusion, there exist a large number of novel and interesting effects in exper- 
iments involving hard scattering or diffractive scattering which can test the ideas of 
color transparency and nuclear filtering. 
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3.6 Figure Captions 

3.1 Kinematics of the BNL pA —> p'p"(A — 1) experiment. 

3.2 The out of plane projections of the missing nucleon momentum distributions 
N(p y ) taken from Ref. [62j. Fig. (3.2a) through (3.2e) correspond to 6 GeV 



and Fig. (3. If) corresponds to 10 GeV incident momentum. 
3.3 The momentum dependence of the transparency ratio for the process pA 



p'p"(A — 1) reported in Ref. [28]. The dashed line at the bottom of Fig. (3.3b) 



corresponds to Glauber calculation. 

3.4 Kinematics of an eA — ► e'p(A — 1) experiment. 

3.5 The momentum dependence of the transparency ratio for the process eA — > 
e'p(A — 1) reported in Ref. ||109|| (filled data symbols). Open symbols are from 
Ref. IB3I. 



3.6 Comparison of the SLAC NE18 experimental results (Ref. |p3}] ) for the trans- 
parency ratio T' for eA — > e'p(A — 1) for Carbon with several theoretical pre- 
dictions. The dotted curve is the prediction of the naive parton model of Ref. 



ijj and the dashed curve is the correlated Glauber/quantum-diffusion model 



of Ref. [0 . The solid line is the Glauber prediction. 

3.7 The missing energy (a,b) and missing momentum (c) distributions for the SLAC 
eA — > e'p(A-l) experiment taken from Ref. [Q. Solid curves are Monte Carlo 
predicted with conventional nuclear spectral functions. No radiatve corrections 
have been done. 



3.8 The theoretical predictions of Farrar, Liu, Frankfurt and Strikman, Ref. j|T| for 
pA —>■ p'p"{A — l). The model predicts a monotonically increasing transparency 
ratio in clear conflict with the data, especially for the Al target (not shown by 
the authors; see Fig. 3.3). 

3.9 The nuclear cross section T • da pp /dt multiplied by by an overall power law 
factor s 10 . This plot, taken from Ref. p2] |, shows lack of oscillations in the 



nuclear cross section and provides evidence for nuclear filtering. 
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3.10 Kinematics of the Fermilab —>■ VA* experiment, where V is a vector meson. 
The final state A* is not measured. 

3.11 The t'{= t — t min ) distributions for Hydrogen (lower curve) and Calcium (upper 
curve), taken from Ref. The curve for Calcium clearly shows the separation 
into coherent and incoherent scattering processes. 

3.12 Transparency ratio T reported by the E665 experiment PJ for photoproduction 
of p mesons for the case of incoherent scattering from nuclear targets. 

3.13 Transparency ratio T reported by the E665 experiment [^, f4(| for photopro- 



duction of p and mesons for the case of coherent scattering from nuclear 
targets. 



3.14 Transparency ratio T reported by the E665 experiment [10] for photoproduction 
of p' mesons for the case of coherent scattering from nuclear targets. 

3.15 The differential cross section da/dt for proton-proton elastic scattering at fixed 



9 cm . The solid lines are the quark counting predictions. [18 



3.16 (a) The energy dependence of R\(s) ~ s 10 da/dt(pp)\ 90 o for high energy pp elas- 
tic scattering at 90° cm angle, (b) Prediction of oscillating transparency T(s) 



for A = 27 from [|119|| . There are no free parameters in the energy dependence. 
Upper and lower curves show possible effects of nuclear phase shifts on the 
prediction, which become negligible at large A. 



3.17 The A dependence of the transparency ratio T as reported in Ref. 



3.18 Fits to BNL data obtained in Ref. ||| for the A dependence of transparency 
ratio for a) (Q 2 = 4.8GeV 2 ) and b) (Q 2 = 8.5GeV 2 ) using the attenuation 
method. The solid line shows the best fit. For comparison the dashed lines 
show constrained fits using permutations of a e ff and best normalizations. 



3.19 Survival probability P pvv extracted in Ref. |)8| from BNL data by use of the 
attenuation method. The increase with Q 2 is evidence for observation of color 
transparency. 

3.20 The attenuation cross section a e ff extracted from the BNL data. Solid line: 
best fit including the intermediate Q 2 data points reported for Aluminum. 
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Dashed line: a 36 mb strong interaction cross section for comparison. The 
two data points are the global a e /f values extracted from Fig. (3.18). 

3.21 Comparison of the fit to the hard scattering rate N(Q 2 ) to pQCD prediction, as 
given in Eq. (5), including scale ambiguity. The ratio s~ 10 N(Q 2 ) / (da /dtpQco) 
is plotted as solid lines; a flat curve indicates agreement. The lower line uses the 
scale of the running coupling as (Q 2 = 0.5(-t)) in calculating da / dt P QCD- The 
upper line uses a s (Q 2 = 1.5(-t)). Data points are for the complete data set (box 
symbol) and the Aluminum data (no symbol). The long dashed curve represents 
s~ 9,7 / (da I dtpQCD) and the short dashed curve is the asymptotic prediction [|9l 



for da/dt P QCD using ot s (Q 2 = —t) and divided by the right hand side of Eq. 
(3.7). 

3.22 The survival probability/constant function of Q 2 /A a for the BNL data 
for three different choices of a. The hard scattering cross section has been set 
to be Z times the short distance pQCD prediction, given by Eq. [3.7]. The 
data points are circles (Q 2 = 4.8 GeV 2 ), squares (Q 2 = 8.5 GeV 2 ) and crosses 
(Q 2 = 10.4 GeV 2 ). The data shows scaling with the value of a = 1/3. 

3.23 Prediction of the survival probability P p in quasi-exclusive electroproduction 
versus A using a e ff(Q 2 ) extracted from the BNL data. Note that a change in 
normalization can nearly compensate a change in a e ff- 

3.24 Q 2 dependence of a "transparency ratio" consisting of a A s (Q 2 )P p for Fe under 
the assumption that the nuclear form factor goes like a 2 . Lower and upper 
lines show the range of theory predictions from varying the choice of scale Q 2 
of a s from (— 0.5Q 2 ,) < Q 2 < (— 1.5Q 2 ), where Q 2 is the photon momentum 
transfer. Vertical scale depends on overall normalization of hard scattering 
rate in nuclear medium which has not been specified. 

3.25 (a) Energy dependence of anomalous transparency in 7r° production for A = 



200, from Ref. ||123|| . T is the transparency ratio. Curves are shown for different 
values of parameter A, which is the amplitude for the qq pair to be small. The 
slope of the curves are quite different even for E n o as low as few GeV. b) A 
dependence of the tansparency ratio. 
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3.26 The differential cross section da/dt for pp elastic scattering at y/s = 400 GeV 
(hollow circles) and 494 GeV (filled circles) is shown on the left. On the right, 
y/s = 1482 GeV. Both figures are fit by da/dt = 0.09t" 8 . From ref. |3 . 
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4 Elements of pQCD 



Color transparency was proposed and came to be studied only after a long history of 
theoretical developments in QCD. Many differences in approach between pQCD and 
nuclear physics have caused confusion and perplexity when introduced along with 
color transparency. To save the reader some of the suffering of rediscovering known 
facts, we include here a brief history of QCD and the context of various approxima- 
tions. We hope to separate the solid foundations from the theoretical "opinion of the 
day" (reviewed in other sections). Our discussion will first be qualitative, providing 
background which might be helpful to newcomers. We then go into more quantitative 
and more detailed subtopics in the application of QCD. 

4.1 A Brief History 

In the earliest days of QCD ||102| , |136|1 , the theory was applied only in terms of the 
operator product expansion ( OPE) . The OPE is a method to express matrix elements 
of composite operators (known to be probing short distance) as an expansion using as 
a basis some fundamental operators of the theory. For example, the Fourier transform 
of the non-local operator product of two currents encountered in deeply inelastic 
scattering 

/ dye*" < p | TUMMo)) I P >= Ti(-9,u + ^f ) + ... (4.1) 
J q 

can be expanded in terms of local quark fields by writing 

T^^A^m^-y (4.2) 

where Q 2 = —q 2 , and A N is defined by 

AnP^p^.-p^ = t£<P,s| O w " 2 '-' w I p, s > . (4.3) 

^ s 

The cn{Q 2 ) are so called Wilson coefficients and On are matrix elements of known 
local operators. The cn{Q 2 ) are defined to satisfy renormalization group equations, 
so that their Q 2 dependence is calculable using perturbation theory. Note that the 
procedure assumes that there is only one "large scale", so that the renormalization 
group, which formally controls dependence on an ultraviolet regularization scale /x 2 , 
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can be used to deduce the dependence on the physically measured scale Q 2 . (Problems 
with more than one scale are treated in the parton model, described briefly below.) 
The matrix elements On are nonperturbative numbers coding the details about the 
particular hadronic state on which the current acts. In this, and most perturbative 
QCD approaches, the perturbative theory does not pretend to model the hadron, 
and makes no initial statement about it. Instead, the matrix elements are to be 
extracted from experiments or supplied from some non-perturbative model. This is 
a very important point. 

This separation of what could be calculated perturbatively from other unknowns 
made the first systematic tests of QCD possible. The simplest series of operators 
employ the lowest number possible of quark fields sandwiching a number of kinematic 
operator coefficients, namely (forgetting color matrices and indices) 

N = i>'y f , 1 D tJ ' 2 ...D flN ^ + permutations (4.4) 

where vp are quark fields and = (id** — gA^) is the gauge- covariant derivative 
containing the gluon field A M . The series is organized into a power series in 1/Q 2 . 
The leading terms in this series are called "leading twist". The word "twist" = 
(engineering dimension - spin) originated in the OPE notation for organizing the On 
into like combinations. The terminology "leading twist" and "higher twist" have 
subsequently been generalized away from its OPE origins to simply mean any case 
where a supposed leading order calculation (in terms of a 1/Q 2 expansion) is done, 
in contrast to a supposed power suppressed correction. 

The series can be summed in a suitable gauge where all of the matrix elements of 
the minus component of A vanish, to give 

q«p (x, Q 2 ) = J dy-e ixp+y ~ < ps | My , Or, + )#(0) | ps) > (4.5) 

which is the same definition as the parton corellation function, Eq (2.3). In general, 
hard inclusive experiments at leading order in 1/Q can be related in much the same 
way to corellation functions of quark or gluon fields, which are the parton distribu- 
tions. By using translational invariance and inserting a complete set of states, we 
can relate the matrix element above to a sum of probabilities: 

q a p(x) = £ |< N | MO) | ps >\ 2 5(xp + -p + - N + ) (4.6) 

N 
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This shows that we actually sum over all channels of all numbers of quarks, anti- 
quarks, and gluons in the target. Subsequent measurements show hadrons to be 
immensely complicated objects in the parton basis. The (leading twist) quark and 
gluon distributions turn out to be universal quantities; different experiments can be 
shown to measure the same matrix elements in combinations that can be unravelled 
through theory and experimental deduction. This is the strongest meaning of the 
word "factorization" , which more generally means that the theory can separate the 
process into certain generalized products of two or more subprocceses for study. 

The Parton Model 

The OPE is supplemented by the more intuitive and flexible parton model,^ orig- 
inally due to Feynman and Bjorken [fil ]. In the early days of the theory, the parton 
model was thought to be a weak sister and not up to the rigor offered by the OPE. 
This changed when it was realized that the parton model was actually more general 
and powerful than theOPE, which was formalizing a limited aspect of the parton 
picture. The key papers in engineering this breakthrough were the so-called DGLAP 
articles [Q. This established the physical origin of leading order scaling violations 
to be the substructure of a quark or gluon at a given resolution scale in terms of yet 
more quarks and gluons. It is important to realize that a theoretical state containing 
"three quarks and no gluons" refers to quarks defined at some scale of resolution. 
At a finer scale of resolution, that is at larger Q 2 , the same objects will be seen to 
contain more and more detailed quark and gluon field combinations. A gauge theory 
cannot be beaten and will not converge in the definition of "what is a quark." Gauss' 
law requires that the quarks be "dressed" with their Coulomb field at shorter and 
shorter distances, a procedure extending ad infinitum in resolution and to all orders 
of pertubation theory. 

The parton model is more flexible than the one-scale OPE because problems with 
several scales and kinematic regimes can be treated. For example, the Drell Yan 
process can be measured as a function of the invariant center of mass energy squared 
s, the lepton mass squared Q 2 , the lepton transverse momentum Q T , its rapidity 
y and an azimuthal angle, besides numerous spins. Using all these variables there 
may be regions of "easy" perturbative problems and others with almost untractable 
problems. In general, processes with more than one large scale can be considered, 



§ what we call here parton model is often referred to as the pQCD improved parton model in 
contrast with the original or naive parton model. 
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but problems involving big ratios of scales (such as the limit Q 2 jQ\ » 1) are often 
problematic. 

The scientific process of deduction and measurement of the numerous parton dis- 
tribution functions during these last twenty years | |136| ] is one of the most magnificent 
achievements of physics. Far more than static observables such as the hadron mass 
spectrum, these matrix elements contain secrets about how Nature makes hadrons in 
a confining theory. However, no theoretical model is yet able to reproduce these distri- 
butions from the QCD Lagrangian or a suitable approximation scheme. Conversely, 
these data have not been able to guide theorists toward a physical understanding of 
confinement mechanisms. The reason might be that these observables, coming from 
inclusive measurements, are quantities averaged over many unobserved states. The 
difficult study of exclusive processes, to which we now turn, is a way to make future 
progress in understanding hadron structure. 

Exclusive Processes 

The early study of hard scattering and exclusive processes |18|, |9^] predates QCD. 
By using a generic renormalizable vector theory with dimensionless coupling con- 
stant and Fermionic constituents, certain salient features of elastic scattering were 
extracted with limited model dependence. A key step in this analysis is the definition 
of the "fixed angle" kinematic limit, which is chosen for simplicity of the theoreti- 
cal analysis. In 2 — > 2 elastic scattering, the fixed angle limit means in terms of 
Mandlestam variables: 

s = const. \t\= const. | u | (4.7) 

with the ratios reasonably close to unity. It follows that the fixed angle limit is a 
"one-scale" problem. The resulting dependence on the single scale s (which in terms 
of the lab energy is approximately 2mEi a j 1 ) is relatively simple. This limit is not 
chosen to make the experimentalist's job easy - indeed, it defines a very difficult 
regime of measurements - but it is chosen to make the theory comparatively clean. 

Quark Counting 

The scaling properties of a perturbative theory in the fixed angle limit can be 
done "counting on one's fingers." Consider the diagram shown (Fig. 4.1) for the 
amplitude of qq — ► qq scattering. From the Feynman rules, this diagram has an 



79 



amplitude 

M(2 — ► 2) = g 2 u{ Pi )u{ Pj )u{ Pk )u{ Pl )/t — > g 2 s 1 / 2 s 1/2 /s = g 2 (4.8) 

where u(pi) is a spinor of momentum Pi . The Dirac traces for | M | 2 are readily 
evaluated. Even before doing that, one can anticipate some simple results simply by 
observing that in the fixed angle limit, uu scales like a momentum, which must scale 
like y^. The amplitude is a Lorentz scalar function of scalar kinematic variables s, t, 
and u, which are all proportional. Counting the 4 — u's and the vector propagators, 
then by elementary dimensional analysis, the amplitude is dimensionless and scales 
like "1", i.e a constant in the fixed angle limit. Note that mass terms, if calculated, 
would give a small contribution in the high energy limit: we can drop m 2 /s << 1. A 
mass m term is an example of a "subleading" or higher twist contribution in a power 
series in m^s. 

The power counting of amplitudes for scattering of more constituents follows from 
induction. (Figure 4.1) For example, add two external Fermion legs, and find two 
new propagators: the result by counting dimensions goes like 

M(3 — ► 3) = g 4 uuuuuu/ P / /t 1 /t 2 = g 4 /s (4.9) 

where t\ and t 2 are the gluon line transfers. The same pattern follows for n-constituents, 
which will give an amplitude-squared like s 4_n . Allowing for the kinematic flux fac- 
tors relating the amplitude squared to a cross section (for example, in the high energy 
limit da/dt =\ M \ 2 /s 2 ), then the overall power of s measured in a cross section 
might be used to "count" the constituents. The prediction that the PP — > PP fixed 
angle cross section goes like s~ 10 follows from this counting of three quarks in the 
proton. 

It is sometimes wrongly thought that this methodology assumes that the proton 
state has been modeled to consist of only three quarks. No; the spirit of the approach 
does not contain this assumption, or any particularly detailed model of the proton 
states. When detailed (and sometimes very thorough ||132|| ) analysis is done, it might 



allow more to be learned about the various contributions. If higher order Fock State 
projections are scattered, then their contribution will simply be small at higher s 
(goes the argument). The experimentally available values of momentum transfer 
are not always large enough to make such contributions small, however. Thus, the 
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argument is an asymptotic one, representing an idealized limit which may or may not 
apply to a realistic experimental situation. 

Is quark counting a rigorous result of QCD? Again, No. At finite Q 2 the entire 
theory of exclusive processes is controversial and extraordinarily interesting because it 
involves the details of how quarks are really arranged in hadrons. One reason that the 
quark counting predictions are not the same as QCD is that they are at the kinematic 
level of the old (sometimes called naive) parton model before scaling violations. At 
large Q 2 they are modified in QCD by renormalization group corrections, typically 
involving inverse powers of logarithms pi, 35fl. These represent the running coupling 
and effects of higher orders of perturbation theory. Even more than the case in 
inclusive processes, the logarithmic corrections are quite important because several 
quarks are involved in a typical exclusive process. As we show later in this section, 
the effects of the logarithms are sometimes equivalent to a non-integer power of 1/Q 2 , 
which can confuse the interpretation of experimental results. 

Independent Scattering 

A second reason that quark counting is not the same as QCD is due to the 



process known as "independent scattering" discovered by Landshoff [88]. To see 



how independent scattering works, consider a hard quark-quark scattering in its cm 
system. Suppose one of the final state quarks goes off in a direction given by p' . 
Consider a second pair of incoming quarks, unrelated to the first, which happen to 
collide independently. There is a chance that the second pair accidentally sends its 
final state quark in the same direction as the first. The final state quarks can then 
combine into a hadron. This is independent scattering. Independent scattering was 
anticipated in Section 2, with Fig. (2.6) showing "mirrors" tilted to arrange this. 
It is evidently an unlikely process, but all elastic hard scattering processes are rare. 
We will compare this process to the quark counting process by estimating the phase 
space integrals for two hard scatterings to coincide in direction and center of mass 
motion. 

The power counting of independent scattering goes as follows. (Figure 4.2). The 

outgoing beams of quarks must coincide in direction well enough to make hadrons 

in the final state; any discrepancy is set by the wave functions, which are defined^] 

^In perturbation theory it is necessary to separate bound state properties of the wave function 
from effects of gluon exchange. To avoid double counting the gluon exchange which produces large 
ItT, the bound state wave functions should have large kr tails subtracted. 
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to have small relative kr- The allowed kr values are much smaller than the beam 
energies, so we can approximate them as almost zero. Because each independent 
on-shell quark-quark scattering amplitude scales like 

g 2 uuuu/t (4.10) 

so that the independent scattering matrix element scales like 

[g 2 uuuu/t] n l A ~ g n/2 (4.11) 

up to logarithmic corrections. The scaling behaviour of the desired elastic scattering 
cross section comes then from the integration region constraint on 4-momenta set by: 
S 4 (k l +k 2 — k 3 — k 4 ). There are three large momenta for each scattering, and one out-of 
plane transverse momentum. This component of the transverse momentum is not as 
big as t/s but instead depends sensitively on what the hadronic wave function allows. 
It should be of order C < k\ > x l 2 in the state's wave function, which for purposes of 
counting is the same as CJ < b 2 > 1 / 2 , b being a transverse space separation. 
Each delta function of a big momentum counts as 1/ y/s, since 

S(P-P') ~ s- 1/2 S(x-x'), (4.12) 

where x and x' are dimensionless scaling variables. The overall probability amplitude 
for a pair of quarks to coincide in final state direction to make a hadron scales like 
the product of the delta functions of momentum, namely like C < b 2 > 1//2 (s)~ 3 / 2 . 
To calculate the cross section, recall that 

da/dt = const. | M \ 2 /s 2 , (4.13) 

giving < b 2 > s~ 5 for meson-meson scattering. As s — > oo, this beats the quark- 
counting process, which for meson-meson scattering goes like s -6 . 

Consider next proton-proton scattering. The argument goes the same way, but 
requires another quark-quark scattering to coincide with the first ones. This adds 
three more delta functions of big momenta, so the amplitude-squared is smaller by 
s -3 . Independent pp — ► pp scattering thus has 

da/dt = const.(< b 2 >) 2 s~ 8 . (4.14) 

This again beats the quark-counting process, which (recall) goes like s~ 10 . 
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How did this process manage to evade the power counting of the quark counting 
process? It is easy to show that the number of gluons and internal propagators is 
fewer than that the assumed in the quark-counting induction; the topologies of the 
low order diagrams are not the same. Because both quark counting and independent 
scattering were studied before QCD was established, early discussion focused on 
comparisons with data. At first it seemed as if pp scattering went like s~ 10 , creating 
a puzzle to explain the absence of the much bigger s~ 8 . One argument was made 
that quark counting diagrams might be more numerous and would dominate for that 
reason. However, when compared at the same order of perturbation theory, the 
independent scattering graphs are myriad and re-emerge inside the quark counting 
diagrams. This happens because internal gluons can become "soft": the effects of 
a soft gluon scales like the same power of s as if the same gluon was absent. The 
upshot is that many quark counting diagrams contain a region indistinguishable from 
independent scattering with a soft gluon. Independent scattering cannot in any sense 
be "absent" . Similarly, if "soft" gluons attached to an independent scattering diagram 
should receive enough momentum to be counted as "hard" , the diagram may merge 



into the quark-counting set. It was finally realized [17 that these physically distinct 
processes actually boil down to different integration regions found in the one theory 
of QCD. 

The independent scattering process had a confused history as these subtleties 
were only gradually appreciated. Closely related (and as much confused) is the issue 
of "Sudakov effects" (Section 4.5), at first thought to suppress the independent scat- 
tering regions, but which were subsequently shown actually to force the independent 
scattering to dominate in the limit of s — > oo. We believe that there is rather 
convincing evidence that independent scattering region of QCD has been observed 
and plays a major role in color transparency (see Section 3). However, the subject 
is unsettled, and the interplay of the independent scattering regions and the quark 
counting regions is currently a subject of active investigation. 



4.2 Setting up QCD Calculations: Stage One 

The process of setting up calculations in QCD can be divided into two stages. In 
stage one, generalized blobs representing Green functions are assembled to describe 
the process under study. This non-specific analysis can reveal integrations regions 
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where perturbative QCD may or may not be used, the most dangerous being particles 
with small momentum. A primary tool in setting up this study has been power 
counting methods which can often predict in advance which regions will become 
dominant at high energy. Sterman [ |130|| has given a thorough review which is barely 



summarized here. A rule of thumb is that a process can be "big" only if it develops a 
singularity. From the Landau rules, singularities in loop integrals occur for kinematic 
configurations that are allowed with classical, on-shell particles. (Off shell particles in 
diagrams are allowed to propagate only a short distance; they are "pinched" short.) 
The on-shell momentum space regions correspond to propagation over long distances 
and long times which allow the amplitude to get singular. 

Momentum Space Versus Coordinate Space 

Much of the intuition of the high energy limit comes from going back and forth 
between momentum space and coordinate space ones and trying to reconcile both 
pictures. Because the transverse coordinate is invariant under the boost, its inter- 
pretation is straightforward and essentially unchanged from the non-relativistic one 
134j|. Thus the transverse coordinate space valence meson wave function: 



il>{k,p) = Jdye ik - y <0\T(q(y)q(0))\p> (4.15) 

where Dirac indices have been suppressed, is related to the momentum space one by 
a simple Fourier transform: 

tffy = J d 2 k T ij{k T ,x)e liT - %T (4.16) 

The other coordinates of a particle depend on the approach. In old fashioned pertur- 
bation theory, particles are put on-shell but energy conservation is not imposed. In 
covariant perturbation theory, energy momentum is conserved, but the particles can 
go off-shell in intermediate states. The space time interpretation changes depending 
on the procedure. 

Null Plane Coordinates 

Most calculations are done in momentum space. One sets up a coordinate system 
oriented locally along each particle's momentum, which is large in the "+z" direction. 
The null plane momentum coordinates are p ± = (p° ±p 3 )/^/2. Similarly, the null 
plane space coordinates are x ± = (x° ± x 3 ) / y/2. With the particle moving fast along 
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the "+z" direction, the x + coordinate acts like the time coordinate in the limit of an 
infinite Lorentz boost. 

With these coordinates the invariant dot product is 

A ■ B = A + B- + A-B + - A T ■ B T . (4.17) 

Thus an ordinary plane wave is given by 

exp(—ip ■ x) = exp(—ip + x~ — ip~x + + ixt ■ fir) (4-18) 

Looking at this expression, one can surmise that p~ is the generator of null plane time 
x + , or the "null plane Hamiltonian." This is correct, as can be shown by looking at the 
algebra among generators.] 134 1 For an on-shell particle, the null plane Hamiltonian 



is just 

p- = (p T + m 2 )/2p + (4.19) 

It is not an accident that this has a non-relativistic appearance: a two dimensional 
Galilean group is one of the subgroups of the null plane algebra. The dependence on 
l/p + is an elementary result of relativity: when p + >> m, then the scale of the null 
plane time generator p~ goes to zero, causing time dilation as a kinematic effect. 

Not much more is needed to understand most arguments using null plane coordi- 
nates in QCD. There is an active industry in calculating things directly in null plane 
perturbation theory - it turns out to be more difficult and subtle than previously 
thought, because of gauge dependent singularities, and difficult interplay between 
the ultraviolet and infrared regions. Covariant calculations in Feynman gauge seem 
to have the least number of pathological singularities and unphysical surprises. 



4.3 Setting up QCD Calculations: Stage Two 

In stage two, blobs are related to wave functions and standard separate gauge invari- 
ant matrix elements. A convenient way to assure gauge invariance of a correlation 
function is to relate it to an on-shell Green function whenever this is possible. At 
this stage covariant integrals are 4 dimensional, something of a conceptual barrier to 
interpretation in non- relativistic quantum mechanical terms. The impulse approxi- 
mation, allows one to finesse some of the ignorance. 

The Impulse Approximation 
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The physical picture of the impulse approximation is a generic feature of Hamil- 
tonian dynamics, cutting across both the quantum and the classical realms. In clas- 
sical physics a particle's state is represented by a point in the classical phase space 
(q(t),p(t)). If at a time t — the system's Hamiltonian changes very suddenly, then 
the particle remains at the same point - namely (q(0),p(0)) - and just time evolves 
according to the new Hamiltonian. It does not have time enough to move across 
phase space, but after the impulse, simply evolves according to the new rules. 

In quantum mechanics the mechanism is much the same. The most convenient 
way to see this is to use the Schroedinger picture. In this picture the dynamical 
phase space coordinates of the system are given by the wave function. If the old 
wave function for time t < evolves to | s(0) >, and if at t > we suddenly turn on 
a perturbation, then the time evolution of the new state | s'(t) > is simply given by 
the old state evolved with the new Hamiltonian H'\ 

| s '{t) >= e- im | s(0) >= e- lE '^ | E' n >< E' n | s(0) > . (4.20) 

where | E' k > are the energy eigenstates of the new Hamiltonian H' . The probability 
that a particular transition occurs in the impulse approximation is simply given by 
wave function overlaps from the initial to the suddenly perturbed system. 

Note that we do not have to follow the details of time evolution during the instant 
the sudden perturbation is switched on, so long as this change is much faster than 
the time scales already present. Sometimes this causes confusion, and opportunities 
to misuse the energy-time uncertainty principle. It would be a mistake, for example, 
to assert that if the perturbation is turned on in a short time At then the system 
after the perturbation will have energies ranging up to 1/At by the uncertainty 
principle. Looking at a generic low-energy wave packet we would see almost no high 
energy components, revealing the error. The facts do not change when the system is 
relativistic: in that case, the underlying Hamiltonian structure still exists, although 
the form of the Hamiltonian must be consistent. 

The Impulse Approximation in QCD: Factorization 

Let us now see how the impulse approximation separates the hard scattering from 
the rest of the dynamics. In Feynman diagrams the configurations of the system to 
be summed over are represented by momentum space integrals. The step where the 
impulse approximation is used is in the integration over each particle's respective k~ 
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when it occurs in loops. Since the light-cone "time" coordinate is x + , then integrating 
over k~ sets the internal times of the system to a value, namely zero. 

With this in mind, consider the diagrams for a hard exclusive reaction shown in 
Fig. (4.3). Let the scattering amplitude be called M. Following the Feynman rules, 
we have a d 4 k = dk + dk~d 2 kT integral for every loop. The momentum k flows through 
the wave functions, and into the perturbative part of the diagram (H) where hard 
gluons are exchanged. Schematically the diagram can be written 

M = j d 4 k l d 4 k 2 d 4 k 3 d 4 k^{k l )^{k 2 )^*{k 3 )^*{k 4 )H{k i ,Q) (4.21) 

where we omit for simplicity the dependence on external variables. Now, if the 
impulse approximation can be used we can make a power series expansion of the 
hard scattering kernel H in powers of k~ and integrate each wave function over its 
k~ . This will be a good approximation if the typical k~ in the wave functions is small 
compared to momenta occuring in H, which is what was meant by hard scattering. 
(Recall that k~ scales like m 2 /2k + for an object of invariant mass m 2 .) The Lorentz 
boost properties make the decoupling in k~ a generically good approximation. 

This leaves the transverse momentum integrals. The independent scattering re- 
gion is tricky and will be discussed separately, //the hard scattering kernel allows a 
power series expansion to be made in kr, as assumed in the quark-counting region, 
then we can again de-couple most of the integrals, obtaining 

M~ J dx 1 dx 2 dx 3 dx4®(x 1 )®(x 2 )®(x 3 )$(x 4 )H(xi,Q) (4.22) 

We then have relations for the scattering of massless, on shell partons, which are 
gauge invariant, parton-model calculations. 

These wonderful tricks are foundations of perturbative QCD. Ignorance about a 
hadron's quark wave function is permissible. Just as in deeply inelastic scattering, it 
is enough to make a list of what the wave functions might be, and arrange to have 
the experiment measure the wave functions. Certain general features of the wave 
functions will be discussed below. 

Power Behavior of Exclusive Processes 

The H part of the quark counting diagram is now evaluated at the on-shell point 
k~ = kx = 0. This is gauge invariant in perturbation theory. If we count the 
dimensions of the diagram and the number of far off shell gluons and fermions, we 



87 



obtain the quark-counting results for the scaling of the diagrams with Q 2 . The 
power counting was designed to be simple and crude and anticipated the QCD- 
based analysis. The overall normalization is not determined; it is an important, but 
extremely complicated quantity. It depends not only on the sum over diagrams, but 
also on the distribution amplitudes, defined next. 

The Distribution Amplitude 

The object 



is called the distribution amplitude, and was invented by Brodsky and Lepage |JT], [35|| . 
As Q 2 —>■ oo, the distribution amplitude tends to evaluate the wave function at the 
point of zero transverse space separation, by Fourier analysis. One can see this 
explicitly by transforming the integrals to obtain 



One sees that the Q 2 dependence of <p(Q 2 ,x) "maps out" the 6-dependence of the 
wave function in the region b 2 > 1/Q 2 . The wave function in question has had its 
hard perturbative "tail" at large removed, and put into the hard scattering. The 
working assumption is that the wave function is unknown but smooth, characterized 
by hadronic physics scales of a Fermi or so. There is also information available from 
the renormalization group, which can organize the limit of extremely short distances 
in powers of £n(b 2 ). For laboratory values of Q 2 it is unlikely that this information 
is enough and it is necessary to model the wave function. Finally one is left with the 
^-integrals, which flow through the diagrams and cannot be eliminated. 

We have re-obtained a main result: the quarks participating have been shown 
to be separated by a transverse distance b 2 > 1/Q 2 . The only essential change 
from the non-realistic mirror analysis of Section 2 is the trading of two Lorentz 
invarient "b" variables for three non-realistic "r's". As in the non-relativistic case 
short-distance is a statement about a dominant integration region, not about any 
compresive preparation of the initial state. In free space scattering, for example, the 
cartoon picture is shown in Figure (4.3.5). 

Wave Functions 

A priori, one knows little about quark wave functions except for symmetry prop- 
erties. It seems unjustified to assume on the basis of non-relativistic quark models 




(4.23) 




(4.24) 
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that the pion would be dominated by the "s-wave" (m = 0) component, since non- 
relativistic quarks are some kind of quasi-particles and not the same as light-cone 
quarks. Let us make a list of the quark wave functions for the simplest case, a pion. 
The wave function to find a quark- ant iquark pair in a pion with momentum p is 
actually a Green function, the Bethe Salpeter two point amplitude |5B |: 



<0|T^ Q (0)^(x,6 T )) I Ps>= {^7s + B\p,p] + ^75 + ^75} ^ (4.25) 

The coefficients A — D depend on position while the Dirac matrices represent the 
spin projections. (The expressions are not gauge-invariant, because they are building 
blocks to be used in a particular procedure.) Now assume that the pion is moving 
fast, namely its momentum p » 1 GeV. The first two wave functions (A, B) scale 
linearly with p and are the big ones by "power counting" - namely, they come from 
the parts of the quark Dirac operators which get big in a boost. The other wave 
functions are kinematically sub-leading. The list is also organized by powers of b; it 
makes it easier to pick out the short distance (A, D) components. There is a similar 
list for the Lorentz covariant vector meson and baryon wave functions. 

Angular Momentum 

One can also make an expansion in L z , the SO (2) orbital angular momentum of 
the quarks about the particle's momentum direction, according to 

i,(b,x) = Y / e imip ^ m (\b\,x) (4.26) 

m 

By continuity and differentiability a wave function with orbital angular momentum 
m has to go like b m as b —>■ 0; this will be important below. The normalization of each 
component of the wave function is a dynamical question. One thing from experiment 
that is known is the normalization of the pion wave function from pion decay: 

^(0) = J -J^y-^T, x) = U- const. (4.27) 

This is a solid number because to a very good approximation the H^-boson only 
couples to a pair of quarks. 

Logarithmic Corrections 

As mentioned before, the overall power behavior of the QCD process is kine- 
matically the same as the QED models used before QCD was invented. There is a 
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substantial difference due to the logarithmic corrections. The typical result can be 
understood rather easily. Postponing the discussion of Sudakov double-logarithms 
to Section 4.5, let us see how leading logarithms are resummed into the distribution 
amplitude, in a way which is much reminiscent of the now standard DGLAP evolu- 
tion equations |53| for structure functions in deep inelastic reactions. For simplicity, 
we restrict to the simple meson case. In axial gauge, one may isolate big logarithmic 
factors in ladder-type diagrams, leading to an expansion: 

$LLO, Q) = +K [ dlL Vqq^qq(u, x) <f (u) (4.28) 

J 



K 2 rl 



where 



— / duV qi j^ q q(u, x) \ du'V q q^ q q(u' , u) <p (u') + ... 

21 Jo Jo 



2 a s {p?) 
k = — In 



P « 5 (g 2 ) 



with p = (11 - ^)/4, rif being the number of quark flavors and V qq ^ qq is the 
following kernel 

Vqq^gq(u, x) = - < - (l 



, , ( 1 + ) 6{u -x) + -(l + ) 9(x -u)\, (4.29) 

J I V u — X J + X \ x — uj + I 



where the ()+ distribution comes from the colour neutrality of the meson. 
The equation on $ may be rewritten as 

( d§\ r 1 

\K) = Jo duV ( u > x )®( x >Q)> ( 4 - 3 °) 

the solution of which has been known for a long time as: 

$(x, Q) = x{\ -x)J2 ®n(Q)C^ 2 \2x - 1); (4.31) 

n 

where the are Gegenbauer polynomials which satisify: 

" du u(l - u) V(u, x) C^ /2 \2u - 1) = A n x(l - x)C ( * /2) (2x - 1), (4.32) 
with eigenvalues A n . One thus gets 

$„(Q) = <P^)e A ^ = $(/i) [^^j , (4-33) 
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where exponents monotonically decrease, beginning with: 

-j- = 0, -j- = -0.62, . . . (4.34) 

Using the previously derived normalization: 

f 1 dx$(x, Q) = $ (Q) t dxx(l - x) = ^ = A, (4.35) 
Jo Jo 6 

the expansion may be written as: 

$(z, Q) = 6f w x(l -x) + $ 2 (lnQ 2 )- a62 a;(l - x)(2x - l) 2 + . . . (4.36) 

We thus know the asymptotic pion distribution amplitude: 

$(x,Q) ~q-k>o 6/^(1 — x). (4.37) 

We do not know the realistic distribution amplitude at accessible energies, since the 
constants <$>2,---*&n cannot be derived from perturbative QCD. This is where the 



QCD sum rule approach enters [29| to yield values for moments of the distribution 



[ dx(2x- l) 2n $(x,/i), ... (4.38) 
Jo 



and then reconstruct the constants $2,.-. ( 3?n. In practice, the reconstruction of 
//) is fairly model dependent, as research in this area continues. 
A similar result is obtained in the baryon case. The proton valence wave-function 
can be written as a series derived from the leading logarithmic analysis, similar to 
Eq. (4.35) for the pion case, but here in terms of Appel polynomials |29| : 



Q) = 120-rj ,-o {1 + ^ (^4t£) 1 + I f^S) 2 + - 

2 \(ys{Qo)J 2 \®s{Qo)J 

(4.39)' 

Here the slow Q 2 evolution entirely comes from renormalization group factors as(Q 2 ) x , 
the Aj being calculated to be increasing numbers: 

20 , 24 

Xl= 9P ' A2= 9^ 6tC -' 

and PA tabulated Appell polynomials 

P 1 (x i )=x 1 -x 3 , P 2 (xi) = 1 - 3x 2 
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Ai are unknown constants which measure the projection of the wave-function onto 
the Appell polynomials: 

Ai = dxidx 2 dx 3 8{x\ + x 2 + x 3 — 1) (j>(xi) Piixj) (4.40) 
Jo 

Independent Scattering, Again 

Recall that the quark counting analysis contains an assumption which is not true 
in general: the assertion that the transverse momentum integrals can be decoupled 
from the hard scattering. Because of this problem, the resulting power counting, etc., 
are doubtful. Instead, quark-counting should be understood as a particular model for 
QCD, a model defined by the region it emphasizes. The complementary region where 
the transverse momentum integrals do not obey the quark-counting assummptions 
is just the independent scattering region. 

In this region a new singularity occurs near = as one might anticipate, be- 
cause the process goes with all particles on-shell. This "Landshoff-pinch" invalidates 
the assumption that H can be usefully expanded in a power series in kx, and conse- 
quently invalidates the factorization a la Brodsky - LePage. For some time it seemed 
that no factorization existed and that the problem might be beyond the "divide and 
conquer" strategy of standard pQCD. This was finally solved by Botts and Sterman 



[ 17] , who realized that a different, new form of factorization had to be created to 
accomodate the different region. 

The basic trick is quite simple, and intuitively appealing. Given that the ampli- 
tude contains unavoidable convolutions over transverse momentum, it is a good idea 
to express this in terms of the transverse spatial coordinate, where one will obtain 
products. Thus, instead of decoupling the k T integrals, one just Fourier transforms 
from kxi to 6j and gets the fixed angle amplitude: 

r i 4 
M(s,t)~/ dXidbtHdxiYjJlipiixubi), (4.41) 

A cartoon going with this expression shows that one has a simple overlap between 
hadrons with different impact parameters. (Figure 4.3). The terms inside the expres- 
sion just have to be evaluated, much the same way as the ^-integrals simply have to 
be calculated. 

The essential feature of the impulse approximation remains. Although the wave 
function integrals are more complicated, their slow bound state time-evolution is 
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successfully separated from the fast time evolution of the hard scattering. There 
is no reason why pQCD should not apply to the resulting amplitude, albeit with 
more detailed dependence on the transverse overlap integrals. As it turns out these 
are extremely interesting, because they contain information we want to know about 
where the quarks might be found. We will return to discuss this in more detail below. 



4.4 Hadron Helicity Conservation as a Test of Models 

The almost chiral nature of QCD without heavy quarks leads to interesting simple 
properties at the perturbative level, namely the fact that a quark or antiquark con- 
serves its helicity when emitting or absorbing a vector particle, a photon or a gluon. 
Long time confinement physics allow however orbital angular momentum to mix with 
quark spins to construct the hadron polarized states, so that polarization observables 
at the hadron level remain an interesting physics question. The fact that hard ex- 
clusive reactions probe the short distance structure of hadrons leads to interesting 
consequences. We will first present this in the context of the assumption of the quark 
counting model, and then examine the differences caused by independent scattering. 

Hadron Helicity Conservation in the Quark Counting Model 

Within the quark counting model, note that a particular orbital angular momen- 
tum component of the wave function is selected by Eq.4.21. Expanding the wave 
function according to Eq. 4.22, then only the m = term contributes in Eq.4.23. 

The wave functions used in this approach therefore are selected by the hard scat- 
tering to give the hadron helicities A n as the sum of the quark helicities 

A„ = E A ^ ( 4 - 42 ) 

i 

In the high energy limit the quark helicities do not flip when gluons are exchanged: 
this is because pQCD is almost perfectly chirally symmetric. It is quite easy to see 
this. The rule for Dirac quark spinors is 

7 5 w(p, A) = Xu(p, A) ; fZ(pA) 75 = -Xu(p, A) . (4.43) 

One can insert 1 = A75 on an initial state spinor, then permute 75 through all the 
Dirac matrices until it reaches the final state. One immediately finds that the "ft 
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term in propagator preserves helicity, while the mass term flips it. Terms of order 
m/p can be dropped when m is only a few MeV. Putting these results together gives 



the famous hadron helicity conservation rule |l"9fl. In a reaction between hadrons 



A + B — > C + D, then the sum of the helicities going in to the reaction is the same 
as the sum going out: 

+ Ab = Ac + As 

This rule is evidently as general as the factorization; it is an exact dynamical 
symmetry of the quark-counting model. 

The hadron helicity conservation rule is a remarkable achievement: it allows us 
to test the general model independent of detailed assumptions about the unknown 
wave functions. Unfortunately, its success when confronted by experiments is very 
uneven. There currently seems to be no evidence of strong violations of the rule when 
applied to electromagnetic form factors. Yet in experiments involving scattering of 
hadrons, it is violated in almost every case tested! Despite this checkered record, we 
will later discuss reasons to believe that the hadron helicity conservation rule will be 
important in the study of color transparency. 

Hard Scattering Hadronic Helicity Flip 

Recall that the hadron helicity conservation rule represents an exact symmetry 
of the quark - counting factorization, because the hard scattering is "small" and 
"round". But perhaps it is not true in free space, if independent scattering is really 
contributing there. The crucial question whether (or not) the symmetry of the model 
is a property of the entire perturbative theory. For this we follow the treatment of 
Ref. 70. 

First, note that the non-perturbative Hamiltonian of QCD does not conserve spin 
and orbital angular momentum separately, but instead generates mixing between 
them. This Hamiltonian operates over the long times in which wave functions are 
formed. Its effects are largely unknown; is why a list of all the orbital and spin 
projections in a meson wave function was given earlier. Thus if a non-zero orbital 
angular momentum component somehow enters the hard scattering - and this is a 
crucial point - then the long-time evolution before or after the scattering can convert 
this angular momentum into the observed hadron spin. It is not necessary to flip a 
quark spin in the hard interaction, because the asymptotic hadron spin fails to equal 
the sum of the quark spins. Such a mechanism is totally consistent with the impulse 
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approximation. 

It turns out that the independent scattering processes are not "round" but instead 
are "flat" (Fig (4.2)). The origin of the asymmetry is kinematic. As before let x be the 
direction perpendicular to the scattering plane and y be a vector in the scattering 
plane. In the independent scattering mechanism, the two (or more) uncorrelated 
scattering planes are separated at the collision point by a transverse out-of-plane 
distance b x . The in-plane transverse position separations Ab y that contribute are as 
small as possible, namely of order the uncertainty principle estimate 1/Q, because 
this is the direction of large momentum flow carried by the gluons. The kinematic 
"short-distance" in the problem is, however, only the in-plane distance. The relevant 
out-of plane transverse distance is set by the hadronic wave functions. 

Returning to the kinematic discussion , the "out-of plane" transverse momentum 
is the main focus in studying independent scattering. To control this variable, it 
is convenient to make a Fourier transform to conjugate space variables b. We let 
the x-direction be out of plane; this is the same for all the hadrons. We let the z/i 
variables be chosen relative to the hadron momenta directions, so b y is in the (vector) 
b x x x pi direction for each momentum Botts and Sterman[T7] then found a useful 
formula for the independent scattering amplitude, which we write in the meson case 
for simplicity: 



Q 

(4.44) 



where 



V a p(x,b;Qv) 



dy 
2tt 



<0\T(q a (y)q j3 (0))\n(Qv)> 



y=y v'+bri' 1 



with Dirac indices aft. 

For the discussion here, note|SB[ that the scattering plane breaks rotational sym- 
metry with the out-of-plane direction x. Immediately one is struck by the absence of 
any selection rule favoring m = 0; instead, all orbital angular in the wave functions 
are allowed. It is as if hadrons "flatten" under impact in the in-plane direction y, 
forming a cigar-shaped hard scattering region. 

We therefore have a new rule of "hadron helicity nonconservation" when indepen- 
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dent scattering occurs: 

Aa + A B ^ A c + Xd (4.45) 

There are a number of tests of this idea. Basically, we want to have spin information 
on every possible reaction in free space, and every similar reaction in the nuclear 
target. When one observes hadronic helicity violation one is observing the orbital 
angular momentum makeup of the hadron! This will be discussed further below. 



4.5 Sudakov Effects 

QED History 

In 1956, V. Sudakov [|131|| studied the large Q 2 dependence of the electromagnetic 



form factor of an electron in perturbative quantum electrodynamics. At one-loop 
order, he observed that the form factor received a "double-log" correction of order 
a£n 2 (Q 2 / A 2 ) . Here A << Q is an infrared cutoff of some kind. These terms come 
from the integration region where the internal photon loop momenta are "soft" , with 
momenta much less than Q. At two loop order the same contribution was squared, 
and multiplied by 1/2!. Because £n 2 (Q 2 ) is a large factor, standard perturbation 
theory in a is inadequate when a£n 2 (Q 2 ) becomes large. 

For this reason the soft photon region is calculated to all orders in perturbation 
theory. Conversely, it can be shown that this region is the only one giving terms 
of order [a£n 2 (Q 2 / A 2 )} 3 . Sudakov showed that the vertex function from summing 
terms of this kind to all orders is the exponential of the one-loop graph, times the 
bare vertex. The "Sudakov form factor" goes like exp(— a const £n 2 (Q 2 /A 2 ). Note 
the minus sign in the exponential. For large Q 2 the form factor goes to zero faster 
than any power of Q 2 \ 

Sceptics have pointed out that the procedure is awkward, which is true enough. 
The value of the sum to all orders is smaller than some of the low-order terms ne- 
glected, causing worries about self-consistency. Although asymptotic trends can be 
identified, quantitative values of leading-log expressions must be interpreted with 
care. A helpful rule is that at a given order of approximation, £n(zQ 2 ) = £n(Q 2 ) 
for any fixed number z as Q 2 goes to infinity, up to subleading corrections. These 
technical issues have dominated much work, but the physical picture that emerged is 
beautiful enough to justify the procedure. 
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Physical Picture 

Sudakov correctly interpreted the rapid decrease of the form factor with Q 2 as 
follows. At large Q 2 , a struck electron has a chance to radiate an enormous number 
of photons into a very large phase space. Because the phase space is so large, these 
inelastic events dominate the probability. There remains very little probability for an 
elastic event. This explains qualitatively why the elastic form factor must be small 
when Q 2 is large. (As a check on the calculation, the calculation of probability "one" 
minus the inelastic events does indeed recover the double-log exponentiated. 1 14]. 

The physical picture indicates that the phenomenon of a strongly damped elastic 
process is very general. It is important to realize that each (on-shell) charged leg 
emerging from any Green function can radiate vector particles copiously, so that 
Sudakov effects are not generally described by simple form factors. Sudakov's form 
factor was simply the first example of far-reaching phenomena, now called "Sudakov 
effects" that must occur in many processes in a gauge theory with massless particles. 

The result seems to be deeper than its modest perturbative origin.. The exponen- 
tiation is not an accidental pattern, but closely related to the gauge transformation 
properties of the matter fields. This follows from the observation of Grammer and 
Yennie [5% that the polarizations of soft gauge photons are primarily longitudinal, 
and can systematically be summed using Ward identities. These become phases on 
the matter fields, because this is how the fields transform. In low orders of perturba- 
tion theory QCD is much like QED and the classic procedures of infrared QED can 
be used almost without change [[JT], |ST|. Higher orders in QCD are tricky, but theory 
can organize the calculations into ordered sums, which add up all terms of the form 
[a J £n K (Q 2 / A 2 )} J , where J and K are related by some pattern. The case of J < K 
is called subleading logarithms. Implementing this systematically in QCD took quite 
some work to accomplish, but the underlying gauge symmetry makes it possible. The 
techniques were eventually brought to a high level of perfection by Collins and Soper 
[ j30|| . Further progress has been made by Sterman ||130|| . 

Running coupling and the chromo- Coulomb phase shift 

Here we will present a fast, heuristic way to understand the exponential functional 
form physically. Consider an e + e~ pair approaching each other before they produce 
a time-like photon. The electron and positron are not free particles but are accom- 
panied by clouds of soft photons. These can be exchanged before the event, and 
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the low-momentum region will produce a phase-shift on the particles' wave function. 
That phase should be given by an eikonal form, in a gauge theory the line integral 
being ie J A - dx where A is the vector potential. For A one might semi-classically in- 
sert the the Coulomb potential, but more generally by dimensional analysis there will 
be an integral of the form ie 2 J dx/x between two endpoints. The closest approach 
in space between the particles is of order 1/Q, and the furthest relevant point set by 
the infrared cut-off, of order 1/A. Thus there is a phase proportional to a £n(Q 2 / A 2 ) . 
This can be shown to be the "Coulomb phase shift" [ 135|| . In perturbation theory it 
is recovered from summing diagrams to all orders. 

It is convenient to employ analyticity to reconstruct the function from its analytic 
properties in the complex plane. The analytic function of Q 2 having a branch cut 
on the real axis and going like exp(i7r£n(Q 2 / A 2 ) for Q 2 > is exp(— £n 2 (— [Q 2 — 
iO]/A 2 ). This is the Sudakov form factor already quoted. This motivates the intimate 
relation between the eikonal phase, which is closely tied to gauge symmetry, and the 
exponential form with double-logs. Note that in reconstructing a contour integral 
knowledge of the region of Q 2 < A 2 is needed. This is no problem for QED, but 
indicates infrared sensitivity of quantity in QCD. 

In QCD much the same eikonalization occurs for quarks and gluons |3C, 133]. The 
quark's Coulomb field is modified by the running coupling, which in coordinate space 
can be written a s (x) = const/ £n{l/ x 2 A 2 q CD ). Repeating the calculation above, the 
QCD "chromo- Coulomb" phase effects are 

1/A 

eXP (- K / T faUfcJ ~ 6XP "^"MQVAU)- (4-46) 



i/Q 2 



While this way of calculating requires justification, the result suggests that sums in 
perturbation theory should recover a phase shift going like £n£n(Q 2 / Aq CD ) , which 
is indeed true. A bonus is that bothersome infrared cutoff effects produce a simple 
constant phase that might be ignored. The appearance of Aqcd as opposed to the 
infrared cutoff is an indication of a perturbatively calculable result. The region of 
momentum from the infrared cutoff to a renormalization point can be held fixed, 
while the region from the renormalization point to the big scale evolves with Q 2 by 



the renormalization group. |TT2| Thus, the QCD "Chromo-Coulomb" phases can be 
calculated as purely imaginary anomalous dimensions ||1 1 2|| , on the same perturba- 
tive footing as other anomalous dimensions. These physical arguments have been 
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confirmed by careful and sophisticated work [ 129 , and now there is no doubt of the 
calculability of these phases. 

Experimental Evidence for Sudakov Effects 

Many years ago it was realized that inclusive processes involving large ratios of 
scales would probe Sudakov effects as a novel form of radiative correction in QCD. The 
classic experimental application was Drell Yan lepton pair production at measured 
transverse momentum. Subsequent study has confirmed the qualitative features of 
the theory, creating a convincing example supporting the experimental relevance of 
the effect. Several other processes involving small transverse momentum seem to be 
equally well treated by now. 



Even earlier, Polkinghorne ||116|| had pointed out that the independent scattering 
process occured via nearly on-shell scattering, and its description might require Su- 
dakov resummation in a theory with vector exchange. This statement lead to claims 
that the rapid decrease with Q 2 of Sudakov form factors would eliminate the inde- 



pendent scattering contributions altogether. Mueller ||101|| showed that this argument 
was incorrect, and that independent scattering actually dominates at asymptotic en- 
ergies, as we review below. 

Another approach is to test experimentally for the presence of Sudakov effects 
by looking for the Coulomb phases in processes with sufficient coherence, such as 
elastic scattering. Evidence for the QCD coulomb phase in pp fixed angle scattering 
was presented in Ref. 137. It was proposed that interference between independent 
scattering regions and quark counting regions would explain the energy-dependent 
oscillations, of about 50% magnitude of da/dt, as experimentally observed. Similar 
oscillations have also been observed in the energy dependence of rr — p hard elastic 
scattering. The disappearance of these oscillations observed in the BNL experiment 
is consistent with filtering of the process to short distance. This has been discussed 
in Section 3. 

QCD Analytic Reconstruction, Bin-by-Bin in Transverse Space 

Again one can find the analytic function which has the running coupling-modifed 
phase reported above. To leading log order the function in the exponent has to be 
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meaning that the complete Sudakov factor is 

This expression contains both the infrared cutoff and Aqcd, showing that its calcu- 
lation is sensitive both to unphysically soft gluons and the calculable perturbative 
region. This is the QCD "Sudakov form factor", which shows the effects of the 
running coupling. If the dependence on an infrared cutoff were taken literally, the 
expression would be doubtful, because infrared-divergent gluons do not physically 
exist. Curiously, the rate of change of phase of the expression with Q 2 does not 
depend on any infared cutoff, so at most the scale of infrared cutoff can be in doubt, 
reflecting the leading- log character of the result. 

The role of the infrared cutoff in QCD is actually provided by bound state and 
confinement effects. Given a color neutral "atom" of quarks and antiquarks, then 
there is a tendency toward destructive interference of long wavelength radiation from 
the various color charges inside, just as in QED. This effect can be incorporated by 
considering the Sudakov sums from the whole bound state coherently. By transla- 
tional symmetry, the radiation of a gluon of momentum q from an emitter at the 
origin and an opposite emitter at a transverse position b are added like 1 — e~ tbq . 
Expanding around q = 0, the "monopole" term vanishes, leaving perturbative dipole 
emission into the momentum region q > 1/6. The upshot is that coherence gives a 
cutoff on radiation exactly like an infrared cutoff. 

Up to subleading proportionality constants, for color singlet hadrons we can trade 
the infrared cutoff for the transverse size with A = 1/6, giving a Sudakov factor of 
the form [lTj 

The calculable coefficient c depends on the color algebra, and the number of particles 
entering and leaving the reaction. The chromo-coulomb phase is not shown, but 
takes the form of the exponential of some calculable color matrices. The form of the 
expression above can be considered a transverse "bin-by-bin" model for the results of 
bound state coherence, acting much like a wave function. Presumably it applies well 
to the short distance perturbative region of b — > 0, while application to large distance 
regions of order the hadron size is problematic. Figure (4.4) shows the damping 
effects of this function on large b values of the wave function. 
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In an axial gauge, such leading logs are contained in diagram topologies which 
are wave function- like, allowing one to count their properties without engaging the 
details of the hard scattering. In any gauge the next-to leading logarithms are more 
difficult and will not be found to "factorize" so simply. We refer the reader to system- 
atic procedures in the literature which re-organize the pertubation theory without 
overcounting diagrams and momentum integration regions. 

Sudakov-Damped Regions Versus Sudakov Suppression 

a. Asymptotics 

There is no doubt that the Sudakov effects are real, but there is controversy about 
their numerical importance. Early arguments claimed that the Q 2 dependence would 
rapidly damp out processes, such as independent scattering, associated with Sudakov 
type corrections. These arguments have not been supported by subsequent studies. 
The Q 2 dependence is simply not sufficiently rapid to make strong suppression for 
experimentally accessible values of Q 2 . The term "Sudakov suppression" is obsolete 
when applied in this way. 

A different approach is to find which regions of b are contributing at a given value 
of Q 2 . The interplay of these two scales causes one to depend on the other. First, let 
us note that the regions of large logarithms of b are the regions of 1/Q < b < 1/Aqcd] 
these regions are of interest and experimentally accessible. Second, inspecting the 
Sudakov formula, one observes that the region of b scaling like Q~ A , where A is a 
positive (but not necessarily integer) power, gives a constant, unsuppressed exponent 
as Q 2 is increased. This is a candidate for a dominant region. With this observa- 
tion and a saddle-point approximation, Mueller [101] came to the conclusion that 



independent scattering would become a short-distance process in the limit of asymp- 
totically large Q 2 , namely Q 2 going to infinity in a mathematical limit. In this limit 
Mueller estimated an effective power behavior for the amplitude, showing that the 
contributiuon was larger than the quark-counting one. 

In the same region and with a similar saddlepoint approximation, Botts and 
Sterman JlT] showed that the Chromo-coulomb phase effects disappear. This can be 



anticipated, because a power-behaved amplitude does not have an energy dependent 
phase. 

However, in the absence of truly enormous Q 2 , are the saddle point arguments 
quantitatively reliable? In numerical studies, Botts [|16j subsequently showed that 
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the asymptotic region argument sets in at about Q 2 = 1 TeV 2 , which is far too high 
to be relevant to experiments. 

b. Realistic Values of Q 2 

There is considerable current interest on the regions of b relevant to current ener- 
gies. At accessible energies and momenta many studies [[TJ], 71 1 have shown that com- 
paratively large values of b ranging up to 1/Aqcd make substantial contributions. 
This can be seen by consulting Fig. (4.4). In practice, many authors investigate 
the effects of the Sudaokov factor multiplying models of the "soft" hadronic wave 
functions which extend out to about a Fermi. Li and Sterman were able to use 
the Sudakov "wave function" in a calculation of the electromagnetic form factor, il- 



lustrating this effect [72]. While the regions of very large b are completely eliminated 
by the Sudakov effect, some worrisome regions which cannot easily be called "short 
distance" remain to give sizable contributions. 

If finite values of b do give substantial contributions, then the Chromo-coulomb 
phase should continue to produce oscillations. One can learn much from spin, because 
spin dependence is another piece of evidence that very short distance is simply not 
achieved in real reactions. There is considerable consistency in the picture between 
electromagnetic form factor studies and hadron hadron studies. Recently, the effects 
of all wave functions in titt — > pp scattering have been explored jSBJ showing that 
hadronic helicity violation might also be attributed to non-short distance effects even 
in the presence of Sudakov corrections. 

Color Transparency and Nuclear Filtering as Sudakov Effects 

Recall that in color transparency, big 6-zones of the wave function interact with 
the nuclear medium, while small 6-zones survive. In Sudakov effects, big 6-zones of 
the wave function radiate inelastically into free space, while small b zones survive to 
make exclusive channels. The mechanisms of color transparency and nuclear filtering 
have the same basic feature of "survival of the smallest" . In comparing the two, the 
"filtering" aspects depend on the medium, which can be vacuum or nuclear matter. 
At any particular value of Q 2 , the nuclear medium sould be more effective in selecting 



short distance than the vacuum because of its stronger interactions. [120, 115 



This suggests that repeating the calculation of Sudakov effects in a nucleus, and 
in particular summing pattern in perturbation theory, might lead to a quantita- 
tive, perturbative theory of color transparency. Unfortunately there is no theoretical 
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agreement on how to model nuclear matter in QCD. 

Let us consider this idea in a medium with the average properties of nuclear 
matter, but lacking some of the detailed correllations. A perturbative gluon can go 
into free space, or be absorbed and re-emitted by spectator nucleons. The result 
is renormalization of the gluon propagator in the nuclear medium. There are some 
constraints to be respected in modeling this. The effects of the nucleus will break 
Lorentz symmetry while remaining gauge invariant. In a derivative expansion, one 
could find terms in an effective action £ e ff of the form 

C eff = eE 2 /2 + B 2 /2/2 + - 

where E and B are the gauge-covariant electric and magnetic fields, and e and /i are 
contants representing the medium. In principle e and fi can depend on position. From 
causality the theory would be non-local, a further complication. Let us investigate 
the local form above anyway. In the temporal Aq = gauge defined by the vector 
n M = (1, 0, 0, 0), the renormalized gluon propagator takes the form 

d, u = (l/e)l/(w 2 - k 2 + (1/M - l)k 2 )/V (4.50) 

where R^ u is the usual tensor 

Rav — ~9uv H T\ n uk u + nykn) — - -—kj\, v . 

n ■ k (n ■ ky 

(It is well known that other tensor terms can also contribute in a medium.) The point 
of this example is the renormalization of the pole value by 1/e. Screening in QED 
makes e > 1. Antiscreening in QCD would be expected to make e < 1, increasing the 
value of the gluon propagator near the pole. 

Neglecting [Xjie^i) — l)k 2 compared to k 2 in the denominator, we can quickly 
anticipate the results. They are equivalent to g 2 — * g 2 /e, which is coupling constant 
renomalization to a stronger coupling. The sum of Feynman diagrams in leading-log 
region can also be found using the new propagator. It is equivalent to known Sudakov 
corrections replacing g 2 — > g 2 /e in the final answer: a nuclear Sudakov factor of the 
form 
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In Figure (4.4) we show comparisons of the relevant 6-space profiles of this expression 
as e is varied. For any given value of Q 2 , as e is varied to smaller values, the effects of 
filtering become more pronounced, exactly as anticipated in the coherence arguments. 

Although theory cannot yet demonstrate such a formula rigorously, we believe 
that future progress in defining the nuclear medium may lead to results of this kind. 
It would be extremely interesting to know whether color transparency and nuclear 
filtering might be described within such a "leading log" framework, because then the 
theory would be tightly constrained. Leading log sums are probably more reliable 
than the present situation embedded in higher twist concepts. One can hope that 
the perturbative calculability of aspects of color transparency and nuclear filtering 
might be put on a secure footing by more work along such lines. 
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4.6 Figure Captions 

4.1 Fermion - fermion scattering with vector exchange is dimensionless and scale 
invariant in the fixed-angle limit (left). Adding a pair of Fermions connected 
by a hard scattering, the amplitude is smaller by a factor of 1/Q. 

4.2 Coordinate space picture of meson-meson independent scattering. Scattering 
planes (dashed) and directions of quarks (arrows) are independent, but to make 
a final state hadron the important region is the case where they are parallel. 
The intersection of the Lorentz contracted hadrons (initial state pancakes are 
shown) and the region of integration over the transverse separation between 
quarks (b) inside the hadrons is a cigar shaped overlap region. Pancake wave 
functions of outgoing hadrons have been omitted. 

4.3 Factorization in QCD. Small internal momenta (k) circulating in wave functions 
(blobs) are separated from the hard scattering kernel H by integrating over 
them. 

4.3.5 (a) Hard scattering in the impluse approximation selects a region where the 
quarks are close together in a full-sized hadron. (b) Compressed state hypoth- 
esis assumes that the state is quantum-mechancially prepared to exist in a 
small spatal region. 

4.4 (a) The Sudakov factor for meson-meson scattering plotted as a function of the 
transverse separation b, in units of Fermi, using Xqcd — 100 MeV. The region 
which survives after coherent emissions of gluons are taken into account is re- 
stricted to b < 1/A; large regions of b associated with large color dipole moments 
are strongly damped (b) Same as (a) but using g 2 — » g 2 /e as a renormalized 
coupling, with e = 1, 1/2, 1/3, 1/4 (top to bottom). Smaller e leads to stronger 
suppression of large b - regions. Both (a) and (b) have been "flattened" in the 
short distance region b < 1/Q where standard Sudakov expressions should be 
patched into short distance expansions. 
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5 Models of Propagation 



In the preceeding section we have explained and developed the concept of color trans- 
parency from the pQCD point of view. These ideas are certainly valid once we go to 
high enough energies. However it is not clear exactly at what energies this formalism 
will be applicable. It is therefore useful to approach the problem also from lower 
energy points of view which, however, necessarily require modeling. This problem 
has been addressed extensively in the literature. We review in this section some of 
the ideas and models that have emerged. 



5.1 The Glauber Theory Framework 

Color transparency effects represent the deviation of the experimentally measured nu- 
clear cross sections in comparison to Glauber theory |54|, which provides a framework 
for calculating the attenuation of fast moving hadrons through a nuclear medium us- 
ing standard hadron-nucleon interactions. It is based on the assumption that the nu- 
clear potential changes very slowly in comparison to the center of mass motion of the 
hadron. We do not cover the Glauber formalism here since it is reviewed extensively 
in the literature, see for example Ref. [?]. Let us, however, make a few comments. 
Glauber theory has been used as a "bottom line" against which color transparency ef- 
fects might be measured. There exist many different calculations [ITl pO, p>§, M 



of this type for both hadro-production and electro-production experiments. These 
calculations typically use pQCD predictions or free space data for the single hard 
scattering, and Glauber formulation or some variation for the soft scatterings with 
the nucleus. Although minor differences in assumptions might be expected, it is 
somewhat disturbing that the results of these calculations vary considerably, of the 
order of 10-30%, from one another. Nevertheless almost all of these authors agree 
that they do not find much Q 2 dependence of the transparency ratio purely by in- 
clusion of nuclear medium effects, and the results of these "traditional" calculations 
disagree with the BNL hadro-production experimental data. Furthermore, most of 
these calculations remain below the BNL results for the transparency ratio. The 
only exception is the work of Frankel and Frati j|5) who surprisingly reproduce the 
BNL data purely in terms of Glauber calculation (claiming that an important effect 
comes when the finite size of the proton is subtracted from the distance z that the 
nucleon traverses in the nucleon, which drastically modifies their calculation of the 
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suppression factor P(r) = exp(— apZi),in the survival probability). In contrast to the 
hadro-production experiments, the SLAC electro-production data is reasonably close 
to Glauber predictions. 

Color transparency effects can be included in the Glauber formalism, if it is mod- 
ified in several ways. At the simplest level it is sufficient to use a reduced hadron- 
hadron interaction strength. The reduction in hadron-hadron cross-section is pre- 
scribed by QCD and goes like 1/Q 2 where Q is the momentum transfer. However 
many authors [£y], [77], |85| have argued that because of the expansion of a proposed 



small wave packet, the interaction strength will also depend on the position of the 
propagating hadron. At the laboratory energies of the order of about 10 GeV, this is 
claimed to be a large effect not calculable from fundamental theory, thereby making 
model calculations necessary. Whether this is true or not is controversial. Many such 
model calculations will be reviewed. Another important effect, which is crucial for 
hadro-production experiments, is the possibility of sizeable contributions due to long 



distance components in free space scattering ||119|| . These contributions are filtered 
in the nucleus and introduce an additional parameter, which for large nuclei might 
be incorporated as a shift in the hard interaction strength. 

5.2 A Simple Model Calculation 

The basic idea of color transparency is that the process of hard scattering involves a 
region of the wave function of very small spatial size (r « 1/Q) which then evolves 
through the nucleus. Most authors have interpreted this to mean that the entire 
wave packet is small, a "compressed state." Since this wavepacket is not an eigen- 
state of the Hamiltonian, it expands and evolves as it travels through the nucleus. 
Modeling the wavepacket and its subsequent evolution can be treated by introducing 
a complete set of states. The various matrix elements necessary for this calculation 
can be obtained by using some models like the non-relativistic quark model or the 
Skyrme model. Partial information about these matrix elements can be also obtained 
phenomenologically. 

Let us first discuss this approach using a toy model of Blaizot, Venugopalan and 



Prakash |15| which is quite illustrative and highly pedagogical. The main simplifying 
idea is to use a nonrelativistic harmonic oscillator model to discuss the evolution 
of a cc pair produced in a diffract ive photoproduction or electroproduction process, 
assuming factorization of the nonperturbative soft scattering process from the hard 
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scattering part. In this case the straightforward impulse approximation (Section 
2) does not apply to diffractive kinematics. However, Feynman diagrams contain 
momentum transfers of order m c or greater, indicating short distance. It is natural 
to describe the cc pair as a very small wave packet at the time of the hard scattering. 
The interaction of the system with the nuclear medium is then assumed to be of the 
form 

H 1 = eSz (5.1) 

where z is the relative coordinate for the charm pair and £ is a model color "electric 
field." The wave packet at the hard scattering is parametrized by a gaussian with 
size r ~ 1/Q, Q being the momentum transfer in the hard scattering, which is much 
smaller than the size of the ground state wave function. 

The evolution of this wavepacket after the hard scattering is then determined by 
the total hamiltonian, that is the sum of the harmonic oscillator Hamiltonian plus 
the perturbation Hi given in equation (5.1). One then defines the absorption factor 
S such that, 

m =e s 

P(0) 

where P{T) is the probability of finding the system in the ground state of the unper- 
turbed harmonic oscillator hamiltonian at time T, where T is larger than the time 
that the charm pair spends inside the nucleus. In the case when the radius of the 
small wave packet is much smaller than the radius of the ground state wave function, 
S is found to be 

S - : : 2Po z - 9 z , (5-2) 



where po = eSr, r is the time over which the charm pair interacts with the nucleus 
and < p\ >=< <p\pl\4> >• in the limit when the size of the initial wave packet goes 
to zero, < p 2 z > goes to infinity and therefore S — > 0. The model therefore leads 
to color transparency and gives reduced attenuation as the size of the initial wave 
packet becomes smaller. The above result for S (Eqn. (5.2)) is valid in the case when 
the charm pair undergoes no free expansion prior to its interaction with the nucleus. 

The calculation with the compressed state initial conditions can be compared 
with one consistent with the impulse approximation, a factor of exp(— iQ ■ r) [Eq. 
(2.5)]. Markovoz and Miller |96[] have analyzed such a model to test the assumptions 
behind color transparency. They model the "nucleon" as being the ground state of 
an electrically neutral system of two quarks interacting via the Coulomb potential. 
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They find that the basic three assumptions namely, i) selection of a small region in 
a high momentum transfer reaction, ii) reduced interactions between a small "color" 
neutral wave packet and the nucleus and iii) long expansion time of the wave packet 
compared with the time of propagation in nucleus, are satisfied in their model. They 
find a rate of expansion of the wave packet that is inversely proportional to the 
momentum transfer Q and is therefore extremely small. 

5.3 Classical Geometrical Expansion Model 



Farrar et al |^TJ studied both hadroproduction and electroproduction color trans- 
parency experiments within the framework of a classical physics geometric model for 
the expansion of the wavepacket. The authors argue that a small size configuration is 
formed at the hard scattering and estimate its transverse expansion during its travel 
through the nucleus: they propose two models, the first one, called naive, being 
driven by the equation: 

b ~ t ~ {m/E)z ; a abs ~ z 2 

where E/m is a time dilation factor and z the distance traveled during the time t. 
The second model, said to be inspired by perturbative QCD, is given by 



Z 1/2 (T abs ~ Z . 



Their ansatz for an effective cross-section of the expanding particle is 



a eff ~ a tat 



[{z/lf + ^-(1 - (z/lf)}e(l -z) + 6(z - I) 



(5.3) 



where k is 1 or 2 depending on the model, I is an estimated effective length of the 
Lorentz-contracted nucleus,n is the number of valence quarks in the expanding object 
and kt is an average internal transverse momentum. The authors advocate different 
expressions for I in their two models. While they favor k — 1, one notes that a \ft 
rise of b(t) corresponds to a transverse velocity which is superluminal and violates 
causality fTOH , casting into doubt its usage as a classical cross section. 



The merit of this approach is its simplicity and physical motivations. It does not, 
however, pretend to be more than an educated guess, a quasi-classical picture. The 
resulting predictions of this paper are given in Fig (3.8). As discussed in Section III, 
it is no surprise that the BNL data cannot be understood within such a model. As 



109 



for (e,e',p) SLAC data, the Q 2 region accessed is too narrow and the data's precision 
is too limited to allow definite conclusions to be drawn. One can infer that the most 
optimistic models in this framework are disfavoured. 

5.4 Hadronic Expansion Models 

The completeness of the set of hadronic states remains an assumption of strong 
interaction-physics dating from ^-matrix days. It is not clear how this is to work in 
a confining theory. In most works in the hadronic basis, it is simply assumed a priori 
that one may expand a small configuration present at the hard scattering as a series 
of states selected from a discrete and a continuum spectrum. One very interesting 
goal is to understand how coherence among the internal coordinates (destructive 
interference of gluon exchange) might be understood in terms of coherence between 
the external (hadronic) coordinates. To realistically implement this idea requires a 
judicious choice of approximations which may be difficult to justify. This lead various 
groups to propose diverse descriptions of the expansion of the mini-states. 

In the specific example of eA — > e'p(A — 1), Jennings and Miller write the ampli- 
tude Ai a to knock a hadron out of the nuclear shell model orbit a as a sum of two 
terms, 

M a = B a + ST a 

where B a refers to the Born term and ST a refers to the scattering terms shown 
diagramatically in fig. (5.1). The Born term corresponds to no interaction with the 
nucleus and is given by, 

B a =< p\T H {Q)\a >= F(Q 2 ) < p - q\a > , 

where T# refers to the photon-proton hard scattering, q is the momentum carried by 
the photon, p is the momentum of the outgoing proton and F{Q 2 ) is a generic proton 
form factor (taken as the free space magnetic form factor Gm)- 
The scattering term is written as 

ST a =< p\UGT H (Q)\a > , 

where G is the Green's operator for the wave packet emerging from the hard scattering 
and U represents the nuclear interactions of the ejected hadron. This matrix element 
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may be evaluated by inserting a complete set of hadronic states. The resulting 
expression for the scattering amplitude is then written as: 

r e iqZ 

ST a = Y, d 2 BdZdZ'——U Nm (B, Z)G m (Z, Z')e iqZ < B, Z'\a > F mN (Q 2 ) , 

(5.4) 

where the sum over m extends to all the hadronic states. The matrix elements of the 
operator U that appear in this expression are written as, 

U N>m (B, Z)= J d 2 b < N\b > U(B, Z; b) < b\m > , 

where \m > are the eigenstates of the internal Hamiltonian and \N > is the nucleon 
state. The variables B and Z are the transverse and longitudinal distances between 
the ejected hadron and the nuclear center and the variable b refers to the internal 
coordinates of the hadron, the transverse component of the interquark separation. 
The inelastic form factor is given by 

FmN (Q 2 ) =< m\T H (Q)\N> 

Finally the eikonal Green's function G is taken to be equal to 

G m (Z, Z') = 6{Z-Z')- 



e i Pm (z-z r ) 



2ipm 

where p 2 m = p 2 + M 2 — and M is the nucleon mass. 

To compute the scattering amplitude ST a , the interaction matrix C/jy.m is modeled 

aS b 2 

U{B,Z;b) = -i2E?-p{B,Z) w . 

where a is the nucleon-nucleon cross section, p is the nuclear density, bn is the size of 
the hadron and b is the position operator measuring the transverse separation of the 
quarks inside the hadron. The nuclear interaction amplitude, given in equation (5.4) 
can then be computed within some specific model for the hadron wave functions. 

Jennings and Miller [[F4| (JM) illustrate their ideas by using a simple harmonic 
oscillator model in two dimensions and restricting the sums to only the ground state 
iV and the first excited state iV* contributions to the scattering amplitude. The 
contribution of both these states turns out to be important to calculate the scattering 
amplitude. Within this model and using different values for the excited state mass the 
authors concluded that the results of the calculation do not agree with the BNL data 
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on pp scattering. In particular they do not obtain the decrease in the transparency 
ratio experimentally observed above 10 GeV. 

In Ref. |T5| , JM repeated their calculation by including in the effects due to the 
independent scattering contribution yi 19 ] and also the charm threshold proposal . 
They concluded that both of these effects go in the right direction for explaining the 
data, but that their calculation still did not yield a sufficiently large oscillation. 

The authors also performed their calculation using a more realistic treatment 
of baryonic components of the wave packet created at the hard scattering. This 
treatment [7£| involves using as much information as possible directly from the ex- 
perimental data to compute the hadronic matrix elements thereby reducing the model 
dependence. Assuming first that at the hard scattering a point like configuration is 
formed which has no interaction with the nucleus, 



UT H (Q 2 )\N>=0 



and inserting a complete set of hadronic states, they obtained 



9 , 9 < a, Ml\T H (Q 2 )\N > 
dM x < N(q)\a p \a,M 2 x > ' 1 ' " 



F(Q 2 







(5.5) 



where a p is the proton-nucleon total cross section. As the wave packet propagates 
through the nucleus its interaction with the nucleus no longer vanishes. Defining an 
effective mini hadron-micleon cross section, a e g(l), 



a eS (l)=a p + J2 [ dM 2 x <N(q)\a p \a,M 2 x >e 

a J{M+m^) 2 



i('Px -p)l 



<a } M 2 x \T H (Q 2 )\N> 



(5.6) 



where each of the components has picked up a phase factor e lPxl with 



v\ = P 2 + M 2 N - M\ , 



one is then lead to compute <r e ff (I). One thus needs the two matrix elements, for which 
only partial experimental information is available. In particular the relative phases 
of these matrix elements are unknown. These matrix elements are then written in 
terms of cross sections, 



< N(q)\a p \a,M 2 x > DD ~- 



2„DD, 



d 2 a 



a 



dtdM\ 



1/2 
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DIS <a,M 2 x \T H (Q 2 )\N>-- 



1 d?a DIS ( 



a) 



1/2 



(Tm dVtdE 

where DD and DIS refer to diffractive dissociation and deeply inelastic scattering. 
(A subtlety in Eq. (5.6) is that the experimentally measured final states in diffractive 
dissociation and deeply inelastic scattering are not the same, either theoretically or 
experimentally, so going from Eq. (5.5) to Eq. (5.6) requires dynamical assumptions.) 
Further modelling is still needed: one may parametrize the sum of matrix elements 
as a simple function g(M x ) = (M/Mx) ■ The resulting calculated results, however, 
still do not agree well with the data. 

Finally in Ref. fTFfl , JM further improve their calculation by including several 
kinematic and dynamical effects which had previously been ignored. The new effects 
included were a proper treatment of the longitudinal (parallel to the virtual photon 
or incident proton momentum) component of the momentum of the detected protons 
and careful inclusion of Fermi motion. The authors then find that their results do 
agree with the experimental results on pA quasielastic scattering as published by the 
BNL group. 

Anisovich et al. have also performed a dynamical calculation of color trans- 
parency in pA — ► p'p"(A — l) process with the inclusion of the mechanism ||119|| for in- 
terference of soft and hard hadronic components. The authors use the diquark model 
of the proton, which allows for a considerable simplification of the calculation. This 
paper can be recommended as a reasonably complete calculation including proper 
kinematics and attention to detail. Their results, with the inclusion of color trans- 



parency effects, agree with BNL data p5 |. Similar results have also been obtained 



by Kohama et al. [83 



5.5 Nonrelativistic Quark Model Calculations 



Kopeliovich, Nikolaev and collaborators |84| , |85| , [87] , |106| , |107|| have applied the non- 
relativistic quark model to color transparency calculations. The model is expected 
to be closer to reality in the case of heavy quark systems like charmonium than for 
light quark hadrons. Let us take the concrete example of diffractive charmonium 
production |84|, in the process — > ccN to describe this type of calculation. 
Using a harmonic oscillator model to describe the interaction between quarks, the 
Lagrangian of the quark antiquark system in interaction with the nucleus may be 
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written as 



L eS (r,f,t) = L{r,f) + -^a(T T )p A (r(t)) 



(5.7) 



where r is the interquark separation, a is the total cc nucleon cross section, is the 
nuclear density, 7 and v are the Lorentz factor and the velocity of the cc pair in the 
laboratory frame and L(r, f) is the Lagrangian of the cc system in free space. 



Here \i is the reduced mass p = m c /2 of the cc system and the frequency uj = 
(M^> — M^)/2 is adjusted to fit the mass difference of the low lying charmonium 
states. The evolution operator representing the propagation of the cc pair through 
the nucleus is then given by 



This path integral is estimated by breaking it up into several smaller steps such that 
for each step the nuclear density can be taken to be constant. This facilitates the 
calculation since then an analytic expression is available for each step of constant 
density for the case of harmonic oscillator. The full evolution can be calculated by 
convoluting the different steps. 

The remaining unknown is the initial wave function for the cc pair. Neglecting the 
interquark interaction, the authors assume the transverse part of the wave function 
as 

# 7 (a:, k T ) oc (m 2 c + k^y 1 

where a is the light cone momentum fraction carried by one of the quarks. In position 
space b this model wave function is 



where A = 0.536 fm and B = 0.11 fm. In the low energy region, the authors model 
the wave function differently as, 



where the parameter a is allowed to vary in a reasonable range, a = l/m c , 2/m c , 3/m c . 




U 




w const x [exp(-& 2 /^ 2 ) - exp(-6 2 /^ 2 )] 



oc exp(-6 2 /a 2 ) 
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The transparency is then calculated by taking the ratio of nuclear to free space 
cross sections, 

nr, fTP , Jd 3 rp A (r)\ < V f \U\V- m > | 2 

Ni ] ~ A\ < > p 

The final results of Kopeliovich et al's calculation for charmonium are shown in fig. 
(5.2). The authors have also used their formalism to calculate transparency ratio 
for photo-production of p ]8?| and <fi mesons and their radial excitations. The 
results for p and p' are shown in fig. (5.3). In this case the use of a nonrelativistic 
oscillator model is, of course, not very credible. Nevertheless their results for p 
meson were found to be in good agreement with the E665 experimental results. 
However they also predicted substantial difference between the transparency ratio for 
p and p', whereas experimentally not much difference was found between these two 
particles [[|(J. Their predictions for and 0' mesons || for CEBAF energies remain 
to be tested. Finally this procedure has also been applied |85| to the calculation 
of transparency in pA —>■ p'p"(A — 1) in quasielastic collision. Within their model, 
the authors show that their theory does not reproduce the BNL experimental results. 
This result is obtained both with and without the inclusion of independent scattering 
contributions. 



Gardner [51] has also done a detailed study of color transparency in photopro- 
duced charmonium using a non-relativistic harmonic oscillator model. The initial 
state is assumed to be a small wave packet formed out of superposition of many 
eigenstates. The nucleons in the nucleus are modelled as strings which carry uni- 
form U(l) color field. The calculation shows that the quantum mechanical coherence 
effects in the final state interactions are very important. 

In most of these calculations the true initial state is guessed in a resonable way, 
but actually unknown. Because of the sensitivity of the calculation to coherence 
of the wave packet, it might be interesting to try different variations of the initial 
state. One idea which has not been explored is to tune the initial state in various 
ways while observing which choices propagate best. In this way, one might explore 
Nature's quantum mechanically selected states that actually will survive. 

The Skyrme model may be used |37| to test for the existence of color trans- 
parency. Within this model the baryons emerge as solitons within an effective theory 
of mesons. The model is known to describe quite well the low energy properties of 
strong interactions. Although color transparency is a high energy phenomenon such 
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model calculations may qualitatively describe the onset of color transparency at low 
energies. Both the studies find that Skyrme model supports the appearance of color 
transparency. 



5.6 Fermi Motion Effects 



Jennings and Kopeliovich ( JK) have argued that color transparency of the nom- 
inal elastic point is possible only due to Fermi motion. The authors proceed within 
the hadronic basis by taking the specific example of electroproduction. They argue 
that in the limit Bjorken XBj = 1 only the proton can be produced on shell, if the 
Fermi momentum is neglected. Therefore no excited states are produced, and within 
their model it is not possible to produce an ejectile with small size which requires 
interference of several amplitudes. 



This argument appears to be wrong for elementary kinematic reasons ||104|| . The 
excited states needed in the model are immediately created diffractively when the nu- 
cleus makes the "measurement" by being probed by the fast ejectile system. However 
the authors make interesting experimental sugggestions, which focus attention on the 
observability of elastic events at xbj 7^ 1 due to Fermi motion. The authors then 
conclude that the transparency effects become considerably smaller with increase in 
XBj- By using a two component model, JK conclude that color transparency effects 
will be negligible in electroproduction for momentum Q 2 < 7 GeV 2 , whereas consid- 
erable effects are predicted at Q 2 > 10 GeV 2 . JK further argue that the unexpected 
behavior of the BNL (p, 2p) scattering experiment is understandable once the effects 
of Fermi motion (as measured by a momentum fraction x) are carefully included. 
Their argument is that for a fixed beam energy, higher center of mass energy s/s 
corresponds to higher Fermi momenta in the direction opposite to the incoming pro- 
ton. This means that for fixed beam energy, x is larger for larger \/s. This explains 
the decrease in transparency observed at 12 GeV since according to their arguments 
transparency decreases with increasing x. The authors argue that this does not hap- 
pen at the lower energy of 6 GeV due to "strong mixing of eigenstates." The authors 
also concede that the explanation of the 10 GeV points is problematic within this 
framework. 



Bianconi, Boffi and Kharzeev [H], [12|, |13j have also studied the phenomenon of 



color transparency within the hadronic basis. They agree on a crucial role of Fermi 
motion in exciting higher mass states which can interfere to form a small sized wave 
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packet that can exhibit color transparency. 

Others do not agree with these conclusions. Anisovich, Dakhno and Giannini 
[||] calculate transparency in electroproduction experiment using the deuteron as 
the target nucleus. This is one of the simplest systems that can be studied which 
is expected to show some color transparency effects. The calculation involves ap- 
proximating the struck proton and the spectator neutron interaction by means of 
one-pomeron exchange. Anisovich et al also allow for the possibility that the photon 
produces resonances which subsequently emerge as protons after diffractive scatter- 
ing. No energy dependence of transparency due to Fermi motion effects emerge from 
their calculation. They further find that the onset of color transparency is very slow 
in all cases considered, and suggest that for the investigation of color transparency in 
electroproduction experiments both intermediate-energy as well as next generation 
higher energy accelerators will be necessary. 



Nikolaev et al |104| have also developed a detailed theory of Fermi motion and 



have studied the possibility of the strong XBj dependence proposed in Ref. [[72 
Working within the framework of the diffraction operator technique, which treats the 
evolution and the final state interactions of the wave packet consistently, the authors 
find that the XBj dependence is much weaker than what was claimed in Ref. [72]. 

We note that the conclusion of Refs. [ 72 , 12 , namely that color transparency is 
not possible without Fermi motion, appears to be result of the picture that a small 
sized wave packet is formed because of the large momentum transfer process. The 
phenomenon of color transparency, however, does not require the formation of such a 
wave packet. It only requires (Section 2) that short distance components of hadronic 
wave functions dominate the exclusive process under consideration. This is discussed 
in more detail below. Although it is very important to properly account for Fermi 
motion effects in order to extract signals of color transparency, there is no agreement 
in the literature with the conclusions of Refs. [^, [12|] concerning any crucial link 
between Fermi motion and color transparency. 



5.7 Quantum Chromotransparency 

We have so far outlined a variety of approaches which aim to describe color trans- 
parency experiments using hadronic basis or non-relativistic quark models for hadron 
structure. In this section we show that color transparency can be obtained directly 
from QCD in the high energy limit. We further show how the factorization of hard 
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scattering and the nonperturbative wave function is achieved in this limit. A consid- 
erable advantage of the pQCD approach is that it represents a unified theory which 
can be held accountable. The effects of coherence are what the theory predicts, and 
not imposed from the outside by a model. The initial state is also not modeled as an 
outside input, but is obtained from a definite calculation. 

We follow Ref. | |124|| by assuming factorization for the ampitude M between the 
hard and soft components. The amplitude M for a hard scattering with 2 legs (Fig. 
5.4) can be written 

M(Q\A) = {J] / [dxid 2 k T A UjP'fa, k Tf )H(xi, x f , k Ti , k Tf , Q 2 )^{x h k Ti )\ } 
I i j J L J J 

(5.8) 

where H (x, kr, Q 2 ) represents the hard scattering, i/)a(%, &t) represents the hadronic 
wave function inside the nuclear medium, and Q 2 is the characteristic momentum 
scale in the collision. In this equation we are emphasising the internal quark coor- 
dinates and suppressing integration over center of mass hadron variables, indicated 
by curly brackets. The wave function represents the overlap of the short distance 
hadronic wave function with its propagation through the nuclear medium and can be 
expressed as 

iP A (x',k' 2 T ) = |y dxd 2 k T (s(x - x')5 2 (k T - k' T ) - F A s, (^(k T -k' T )] 2 Wz,4)} 

(5.9) 

where ipo( the free space hadronic wave function and Fa is the nuclear scat- 

tering amplitude. This is nothing more than the light-cone region of 1 — Ga used in 
Section 2. For simplicity here we suppress the i, f indices. The nuclear filtering of 
soft components of the wave function will imply that only the short distance part of 
the wave function survives. Therefore the function i/joa(x, kj) may be approximated 
very well by its short distance form, in contrast to the analogous quantity in free 
space. Going to the quark transverse separation b space we get: 

ij) A (x, b) = f A (s, x, b)ipo(x, b), (5.10) 

where the tilde's denote the Fourier transforms and Ja = \—F A is the nuclear survival 
amplitude. In Eq. (5.9) we have suppressed a kernel for transitions from x — >• x' and 
left the x-dependence as diagonal. The distribution amplitude can then be written 

as 

A (x,Q 2 ) = ( Q d 2 k T [ d 2 b T e ih ^~fA{s 1 x 1 b 2 )^{x,b) 
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J o 



The above integral will get its dominant contribution from the b < l/Q region of 
the wave function. This is exactly analogous to the selection of short distance by 
exp(— iQ ■ r) in the non-relativistic case (Section 2). 

We assume that the wave function tp is a slowly varying function of b for small b. 
This assumption is supported by bound state and renormalization group studies. We 
can then approximate the above integral by replacing the wave function by its value 
at b pa l/Q and cutting off the integral at the same value of b. By substituting this 
expression into Eq. (5.8) we see that we have a factorized form for the amplitude M: 



M(Q 2 ,A) = |[dz]{/ A (6«l/Q) 



ijV( Xl ,l/Q)H(x,0,Q 2 )^(x 2 ,l/Q) 



< » l/Q) . . . f A (b-f » l/Q)}H(x, 0, Q 2 ) >, (5.12) 



where the pointed and curly brackets indicate restoration of the integrals over x and 
hadron center of mass coordinates, respectively, which code information on nuclear 
size and density. Note that H(x, 0, Q 2 ) is independent of A and hadron center of 
mass coordinates for A » 1. That allows us to factor the hard scattering out 
of the integral over center of mass coordinates. Generally one cannot factor the x 
dependence of the cross section. One might hope that the nuclear medium effects do 
not much change the x integrals over the wave functions. 
The transparency ratio T can then be written as 

T(Q 2 ,A) = da/dtU 



A d<r/dt\f ree space 

I < {Ui ~ l/Q) . . . h{V ~ l/Q)}H(x, 0, Q 2 ) > | 2 

A(L(J J dt free space 

P N (Q 2 ,A)R(Q 2 ), (5.13) 



where iV is the number of participating particles that cross the nucleus. The first 
line is an approximation based on the expectation that the ^-convolutions more or less 
factor. If so, then the transparency ratio can be factored into two pieces: the survival 
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probability P(Q 2 , A) and the ratios of the hard scatterings R(Q 2 ). Furthermore, the 
leading Q 2 dependence of the nuclear hard scattering can be obtained by considering 
the leading order perturbation theory diagrams convoluted with the short distance 
hadronic wave function. 

Let us now demonstrate that perturbative QCD predicts color transparency at 
asymptotic Q 2 [|124| . For this purpose, let us consider equation (5.11) at large Q 2 , 
and show that the nuclear-filtered wave function reduces to the short distance free 
space wave function as Q 2 — > oo. This may be done systematically by using a Mellin 
transform, defining 

r%Q- N u*,Q) 

Jo Q 



v ; r(i + iv/2y ' v ; 

using the Mellin transform of the Bessel function. Here iPa,n{x) is the impact pa- 
rameter moment ^ 

^Pa.n(x) = —b N tjj A (x,b). 
Jo 

Note that the b moments are defined with the negative — N of the Q moments; this 
is because powers of b will turn out to map into powers of 1/Q (modulo logarithms). 

Inverting <pN, complex iV-plane singularities at iV = are leading twist; those at 
N = —2 give 1/Q 2 higher twist, and so on. The wave function ipA{b) determines the 
N = singularity: 

/•OO rjh 

Mx) = Jo T bNtfj ° m ~ b2Al/3n <s + ■■••] 

Note that the iV — > behavior is fixed by the short distance free space wave function 
ipo(b). Here the effect of interaction with the nucleus is a factor "1". The y4 1//3 <7g ff 
dependence affects the iV — > —2 behavior, and is thus suppressed by 1/Q 2 - This is 
transparency. It is remarkably general in the asymptotic limit Q 2 — > oo. 

To illustrate the idea that at large Q 2 nuclear wave functions reduce to the short 
distance free space wave functions, let us assume a typical example of the hadronic 
wave function ipo{b) = exp(-mb), and model the nuclear filtering amplitude Ja by 
/a(&) ~ exp[— {b 2 A l ^a' cS ) 1 / 2 ]. Inserting this into equation (5.11), we get 



i(Q 2 ) = MQ 2 



y/[m + (AV3^ ff )l/2]2 + Q2 



m 



m 2 + Q 2 
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which goes to the free space distribution (fi (Q 2 ) as Q 2 — ► oo. By adjusting parameters 
m and o' rme ^ many models can be made. More sophisticated models are of course 
possible. The advantage (and purpose) of the factorization/distribution amplitude 
formalism is that certain universal functions are identified. These are definite matrix 
elements, acccessible to short distance renormalization group techniques. If one also 
accepts that Pomeron exchange is calculable for short distance dominated processes, 
(and certainly there is a great deal of work in the literature claiming this is so), 
then the entire effects of interactions, and all of color transparency, are perturbative 
problems. Part of the underlying protection from infrared effects is due to large 
Q 2 , the original motivation |20|, |100|| for color transparency. The other part is large 



A, which also acts like an infrared cut-off in the nuclear filter. 1 119 , 124] A totally 
perturbative treatment is evidently realistic, and a very exciting prospect. However, 
its exploration would take us into new research and remains for future work. 

5.8 Resume 

We have established that pQCD predicts color transparency in the limit Q 2 — > oo. In 
order to compare with experiments, however, one needs numbers for the magnitude of 
this effect at finite energy. The usual approach requires modeling the nuclear medium 
which introduces considerable uncertainty into the calculation. 



As emphasized in Ref. [ 114}j , model uncertainty can be considerably reduced by 
noting that the nuclear effects should follow a scaling law. The scaling law arises 
because the survival probability, being dimensionless, can only depend on a dimen- 
sionless variable. In the high energy limit the only dimensionless variable that is 
relevant is the effective number of nucleons encountered by the hadron, which is 
proportional to naeffA 1 ^ 3 , where n is the nuclear density, and a e ff is the effective 
attenuation cross section. In the high energy limit a e ff is set by the dominant b re- 
gion proportional to l/Q 2 , and therefore the scaling law follows. We have discussed 
applications of this in Section 3. 

In closing this section on Models of Propagation, a comment might be suitable. If 
the various models can be understood well enough to extract other broad regularities, 
such as modified scaling laws or results that are independent of adjustable parameters, 
then clearly their scientific discriminating value might be enhanced. Science would 
be well served by first establishing definite trends, and afterwards fitting models. 
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5.9 Figure Captions 

5.1 The Born B a and first scattering ST a terms for ( r y*p,p) scattering at high Q 2 , 
Ref. 



5.2 The predictions of Ref. [|87| for photoproduction of charmonium. 

5.3 The predictions of Ref. |S7] for photoproduction of p and p'. The experimental 
results 0, H, for these mesons, shown in Figs. (3.12-14), agree nicely with 
the predictions of Ref. [37] for the case of the p meson but fail to agree for the 
P'- 

5.4 Factorization of the hard scattering H(k, k', ...Q) from the soft evolution of the 
wave functions ip and i/j*. 
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